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Chapter  1  Introduction 


The  primary  goal  of  this  research  is  to  fully  measure  and  understand  the  effects  of  various 
shaped  roughness  elements  in  two-dimensional  and  three-dimensional  high  Reynolds 
number  turbulent  boundary  layers.  Roughness  elements  with  heights  less  than  7%  of  the 
boundary  layer  thickness  are  analyzed  in  considerable  detail.  Mean  flow  and  turbulence 
characteristics  associated  with  the  various  protuberances  are  discussed  extensively. 
Another  significant  attribute  of  this  study  is  to  determine  the  effects  of  spatial  smoothing 
of  roughness  elements  on  the  decay  of  turbulence  quantities  in  a  turbulent  boundary 
layer.  This  is  the  first  time  that  a  detailed  fundamental  study  has  been  done  on  roughness 
elements  of  these  shapes  and  sizes.  Moreover,  it  is  essential  to  note  that  detailed  analysis 
has  not  been  done  previously  concerning  the  spatial  smoothing  of  related  roughness 
elements.  The  end  result  of  this  research,  combined  with  the  previous  efforts  of  George 
and  Simpson  and  continued  future  efforts,  is  to  provide  detailed  measurements  that  will 
enable  one  to  obtain  a  deterministic  conceptual  scheme  in  order  to  better  define  how 
roughness  affects  the  physics  of  the  flow  in  turbulent  boundary  layers.  This  study  is  also 
very  significant  in  that  it  will  provide  detailed  modeling  test  cases  which  will  allow 
comparison  with  and  adjustment  to  turbulence  models  that  are  used  to  calculate  such 
flows. 

The  rest  of  this  chapter  is  broken  up  into  five  different  sections.  Section  1 . 1  will  discuss 
laminar  and  turbulent  flows  over  isolated  protuberances.  Section  1.2  will  go  over 
previous  efforts  related  to  turbulent  flow  past  surface  mounted  obstacles.  This  section 
will  look  at  obstacles  that  have  a  height  more  comparable  to  that  of  the  boundary  layer 
thickness.  Previous  literature  on  the  distribution  of  roughness  elements  and  rough-wall 
turbulent  boundary  layers  will  be  looked  at  in  section  1.3.  A  brief  look  into  some  related 
modeling  and  computational  efforts  used  to  calculate  roughness  flows  can  be  found  in 
section  1 .4.  Finally,  the  organization  of  the  research  and  thesis  is  presented  in  section 
1.5. 
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1.1  Laminar  and  Turbulent  Flow  Over  Isolated  Protuberances 

The  earlier  experiments  involving  flows  over  single  protuberances  were  conducted  in 
laminar  boundary  layers.  The  emphasis  of  these  early  laminar  boundary  layer  studies 
was  on  defining  a  critical  Reynolds  number  based  on  the  element  height  which  would 
cause  transition  to  turbulent  flow.  Due  to  the  fact  that  studies  conducted  in  laminar 
boundary  layers  are  not  directly  related  to  the  current  research,  only  a  few  examples  will 
be  discussed  on  this  topic.  Tani  (1961)  studied  the  effect  of  two-dimensional  and  isolated 
roughness  in  laminar  boundary  layers.  Through  this  research  it  was  discovered  that  if  the 
fluid  speed  is  low  and  the  height  of  the  roughness  element  is  small  compared  to  the 
boundary  layer  thickness,  then  transition  will  not  occur  and  the  effect  of  the  roughness 
element  will  only  propagate  downstream  a  short  distance.  Flows  involving  laminar 
boundary  layers  can  have  relevance  to  the  current  research  in  regards  to  the  formation 
and  propagation  of  horseshoe  vortices.  These  vortices  will  only  be  present  in  laminar 
boundary  layers  if  there  is  the  existence  of  a  pressure  gradient  large  enough  to  cause 
shear  layer  roll-up  at  the  upstream  junction  of  the  roughness  element  base  and  the  wall. 
Experiments  that  involve  the  development  of  the  horseshoe  vortex  structure  can  be  found 
in  Sedney  (1973)  and  Gregory  and  Walker  (1955).  Gregory  and  Walker  performed  a 
qualitative  analysis  on  isolated  protuberances  using  the  china-clay  and  smoke 
visualization  techniques  at  wind  tunnel  speeds  of  37  m/s  and  5  m/s  respectively. 
Elements  tested  were  on  the  order  of  the  same  size  as  the  present  research  while  having  a 
roughness  height  to  boundary  layer  thickness  ratio  (k/8)  less  than  1  in  all  cases.  Varying 
heights,  5.1  mm  and  smaller,  of  a  60°  cone  and  cylinders  were  tested.  The  horseshoe 
vortex  structure  related  to  each  element  was  observed.  Sedney  also  investigated  flow 
past  cones,  hemispheres  and  cylinders  using  a  smoke  visualization  technique.  Again,  the 
presence  of  horseshoe  vortices  in  these  flows  was  detected  and  the  qualitative  flow 
features  analyzed.  Other  experiments  that  have  been  performed  on  isolated  roughness 
elements  in  laminar  boundary  layers  with  the  intent  to  study  the  effects  of  flow  transition 
or  the  development  of  vortex  structures  due  to  the  presence  of  roughness  element  are  as 
follows.  Ichimaya  (1999)  performed  experiments  on  a  2  mm  high  cylinder  in  a  2.2  mm 
thick  boundary  layer  in  order  to  determine  the  effects  of  a  single  roughness  on  boundary 
layer  transition.  Barrett  et  al.  (1993)  took  measurements  at  various  plane  locations  in  the 
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downstream  wake  region  of  a  0.8  mm  diameter  sphere  in  a  1.4  mm  thick  boundary  layer 
and  has  shown  the  development  of  vortical  structures.  Finally,  Klebanoff  et  al.  (1992) 
investigated  the  development  of  a  turbulent  boundary  layer  induced  by  a  hemisphere 
placed  on  a  flat  plate. 

There  is  a  very  limited  amount  of  literature  related  to  isolated  roughness  elements  in 
already  developed  turbulent  boundary  layers,  hence  the  necessity  for  the  present  research. 
Previous  studies  have  mainly  been  restricted  to  the  evaluation  of  the  drag  induced  by 
these  small  protuberances.  Wieghardt  (1946)  conducted  experiments  in  order  to  quantify 
the  drag  of  various  sizes  of  cylinders.  These  experiments  were  performed  in  a  67  mm 
thick  boundary  layer  at  a  Rex=7.2(106),  where  x  is  the  streamwise  distance  over  which  the 
flow  develops.  Hoerner  (1958)  presents  results  for  common  protrusions  such  as  bolt  and 
rivet  heads,  brackets  and  various  other  joints.  Relations  of  protuberance  sizes  and  shapes 
to  their  associated  drag  characteristics  were  investigated  and  presented  as  drag 
coefficients  normalized  by  an  ‘effective’  dynamic  pressure.  These  results  were  garnered 
from  Hoemer’s  work  as  well  as  various  other  people’s  work  such  as  Nikuradse, 
Schlichting,  and  Moody.  Hoerner  also  discusses  the  characterization  of  roughness  by 
relating  it  to  the  concept  of  sand  grain  roughness.  Young  and  Patterson  (1981)  reported 
on  aircraft  excrescence  drag  which  also  involved  Wieghardt’ s  data.  These  data  were 
presented  as  a  drag  coefficient  normalized  by  the  undisturbed  friction  coefficient  versus 
the  undisturbed  roughness  Reynolds  number  with  respect  to  the  element  height  and 
undisturbed  friction  velocity,  Ux.  Gaudet  (1987)  examined  the  drag  on  various  forms  of 
three-dimensional  excrescences  including  cylinders  and  mushroom  shaped  rivet  heads. 
These  results  as  well  as  Wieghardt’ s  data  are  presented  as  a  functional  form  related  to 
wall  variables. 

Previous  work  that  does  investigate  the  turbulence  structure  behind  wall  mounted 
protrusions  is  that  of  Fontaine  and  Deutsch  (1996).  The  flow  field  behind  a  wall 
mounted  Gaussian  bump  with  a  height  of  16.4  wall  units  (k+=16.4)  and  a  base  diameter 
of  13  wall  units  was  analyzed  in  a  turbulent  boundary  layer  with  a  momentum  thickness 
Reynolds  number,  Ree,  equal  to  730.  This  Ree  is  based  on  a  skin  friction  velocity,  UT, 
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equal  to  0.4  m/s  and  an  equivalent  momentum  thickness  based  on  centerline  velocity.  It 
was  found  that  the  protrusion  produced  a  pair  of  counter-rotating  vortices,  each  1 5  wall 
units  in  diameter.  Reynolds-stress-producing  events  such  as  sweeps  and  ejections  were 
discovered  to  be  retarded  within  the  convergence  region  of  the  vortices.  The  work  of 
George  and  Simpson  (2001)  is  directly  applicable  due  to  the  fact  that  the  current  research 
is  an  extension  of  these  previous  experiments.  Previous  work  completed  by  George  and 
Simpson  was  performed  on  three  cylinders  and  one  Gaussian  shaped  element  in  the  same 
boundary  layer  conditions  as  are  described  in  Chapter  2.  The  three  cylinders  tested  had  a 
base  diameter  of  1 .98  mm  and  were  0.38  mm,  0.76  mm,  and  1.52  mm  tall.  The  Gaussian 
element  also  had  a  height  equal  to  1.52  mm.  These  experiments  will  be  discussed  in 
more  detail  in  the  following  chapters.  It  is  important  to  note  that  higher  Reynolds  shear 
stresses  were  found  downstream  of  the  elements  due  to  the  horseshoe  vortices.  An 
increase  in  skin  friction  drag  was  also  found  behind  the  element.  This  is  due  to  the 
significant  amount  of  fluid  being  drawn  toward  the  wall. 


1.2  Turbulent  Flow  Past  Surface-Mounted  Obstacles 

Studies  discussed  in  this  section  are  ones  related  to  obstacles  having  dimensions  more 
comparable  to  the  boundary  layer  thickness,  or  larger,  than  those  of  the  protuberances 
discussed  previously.  As  compared  to  section  1.1,  flows  in  this  regime  have  flow 
features  that  are  more  prevalent,  easier  to  distinguish,  and  easier  to  make  qualitative 
conclusions  about.  There  have  been  many  two-dimensional  shapes  examined  including 
square  bars,  triangular  cylinders  and  forward  and  backward  facing  steps  to  name  a  few. 
There  are  a  few  important  differences  between  two-dimensional  and  three-dimensional 
flows  as  is  described  by  Hunt  et  al.  (1977)  and  Martinuzzi  and  Tropea  (1993).  These  key 
differences,  respectively,  are  that  two-dimensional  flows  have  a  closed  mean  flow 
separation  bubble  behind  the  element  whereas  three-dimensional  flows  do  not,  and  flows 
over  two-dimensional  obstacles  take  longer  to  reattach  as  compared  to  flows  over  three- 
dimensional  obstacles.  To  this  end,  only  three-dimensional  obstacles  will  be  discussed  in 
the  rest  of  the  section. 
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Wind  effects  over  three-dimensional  bluff  bodies  has  been  the  center  of  many  studies 
concerned  with  determining  drag  and  pressure  distributions  on  man-made  obstacles  as 
well  as  natural  obstacles.  The  study  of  Okamoto  et  al.  (1977)  was  motivated  by  the  flow 
around  cones  and  how  they  could  be  related  to  the  flow  around  a  mountain.  Cones  with 
vertex  angles  equal  to  60°,  90°,  120°,  and  150°  were  used.  All  cones  had  base  diameters 
equal  to  100  mm  and  heights  equal  to  86  mm,  50  mm,  29  mm,  and  13  mm  respectively. 
Experiments  were  conducted  in  a  turbulent  boundary  layer  with  a  thickness  of  7.6  mm 
and  a  ffee-stream  velocity  of  10  m/s.  Surface  pressure  measurements  on  the  cone  and  flat 
plate,  and  wake  velocity  measurements  were  taken  to  analyze  the  flow  field.  Drag  and 
lift  coefficients  calculated  from  the  surface-pressure  distributions  were  presented  as  well. 
Surface  oil  flows  were  also  performed  to  visualize  the  formation  of  horseshoe  vortices  as 
well  as  to  define  the  separation  points  on  the  cone  surface.  Another  study  by  Savory  et 
al.  (1988)  focused  on  flow  over  a  hemisphere  and  its  relation  to  flow  over  a  domed 
building.  This  particular  investigation  involved  a  hemisphere  with  a  base  diameter  equal 
to  190  mm  in  three  separate  turbulent  boundary  layers  with  thicknesses  equal  to  367,  258, 
and  86  mm.  All  flows  had  a  free-stream  speed  of  10.7  m/s.  Drag  measurements, 
pressure  measurements  and  wake  velocity  profiles  were  taken  in  order  to  quantify  the 
hemisphere’s  presence  in  the  boundary  layer.  Cases  in  which  the  obstacle  was  more  or 
less  contained  within  the  boundary  layer  are  Arie  et  al.  (1975a  and  b),  Taniguchi  et  al. 
(1981),  Sakamato  et  al.  (1982),  and  Sakamato  (1985).  Different  aspects  of  circular 
cylinders,  rectangular  cylinders,  and  cube-shaped  obstacles  were  examined.  Obstacle 
characteristics  such  as  aspect  ratio,  ratio  of  obstacle  height  to  boundary  layer  thickness, 
and  the  ratio  of  friction  velocity  to  free  stream  velocity  were  related  to  the  form  drag  of 
each  obstacle.  Conclusions  stemming  from  this  work  are  that  the  obstacle  shape  and  size 
and  the  thickness  of  the  boundary  layer  are  key  factors  in  describing  the  flow 
characteristics. 

A  recent  study  was  performed  involving  the  flow  around  surface  mounted  pyramids  in 
thick  and  thin  turbulent  boundary  layers,  Martinuzzi  and  AbuOmar  (2003).  The  effects 
of  apex  angle  and  angle  of  attack  were  related  to  wake  periodicity  and  vortex  shedding 
based  on  surface  pressure  measurements  taken  on  the  wall  and  pyramid  sides,  velocity 


Chapter  1:  Introduction 


-6- 


measurements,  and  mean  flow  surface  patterns.  For  broad  pyramids  it  was  found  that 
periodicity  in  the  wake  does  exist  but  cannot  be  related  directly  to  vortex  shedding. 
Periodic  fluctuations  are  seen  in  the  surface  pressure  measurements  on  the  sides  of  the 
pyramids  as  well  as  in  the  velocity  field  in  the  wake.  Also,  for  pyramids  placed  in  thick 
boundary  layers  no  wake  periodicity  is  observed  and  the  wake  characteristics  are 
different  as  compared  to  those  of  the  pyramid  in  a  thin  boundary  layer. 

The  flow  around  a  surface-mounted  cube  has  been  the  focus  of  numerous  studies;  a  few 
of  these  are  Castro  and  Robins  (1977),  Hunt  et  al.  (1978),  Scholjield  and  Logan  (1990), 
Martinuzzi  and  Tropea  (1993),  and  Sousa  (2002).  The  investigation  by  Castro  and 
Robins  is  one  of  the  few  previous  studies  that  has  analyzed  the  effect  that  orientation  has 
on  the  flow  field.  Mean  velocity  and  turbulence  intensities  were  measured  in  the  wake  of 
each  case.  These  measurements  show  the  effect  that  the  vortices  have  on  the  near  wake 
flow  regions.  The  velocity  deficit  in  each  case  was  shown  to  have  fully  decayed  at  a 
location  equal  to  six  lengths  downstream.  It  was  also  discovered  that  the  cube  oriented 
45°  relative  to  the  flow  produced  a  stronger  downstream  vortex  pair  than  the  cube 
oriented  90°  relative  to  the  flow.  The  same  result  has  also  been  found  to  be  the  case  in 
the  present  research.  Schofield  and  Logan’s  analysis  concentrates  on  how  the  major  flow 
features  are  influenced  by  the  model  geometry  as  well  as  the  incident  shear  flow.  Sousa 
utilized  digital  particle  image  velocimetry  to  identify  the  location  of  the  large-scale 
vortical  structures  present  within  the  flow  field  surrounding  the  cube. 

Another  study  by  Logan  and  Lin  (1982)  evaluates  the  wall  shear  recovery  behind  prisms 
of  different  aspect  ratios  at  downstream  centerline  locations.  This  study  concluded  that  a 
quicker  recovery  was  seen  by  three-dimensional  obstacles  as  opposed  to  two-dimensional 
obstacles  due  to  the  presence  of  the  streamwise  vortex  structure.  One  of  the  most 
complete  flow  field  pictures  of  the  flow  around  a  cube  is  given  by  Martinuzzi  and  Tropea 
(1993).  Two-component  laser-Doppler  anemometry  measurements  along  with  surface 
pressure  measurements  and  surface  oil  flow  visualizations  were  utilized  to  describe  the 
flow  field.  The  aim  of  this  study  was  to  determine  the  separation  and  reattachment 
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pattems  on  the  prism  as  well  as  to  obtain  quantitative  data  that  describes  the  mean  and 
fluctuating  velocity  field. 


1.3  Rough-Wall  Turbulent  Boundary  Layers 

Due  to  the  limited  amount  of  directly  applicable  data  obtained  in  the  present  research, 
this  section  will  be  a  very  restricted  look  into  the  previous  research  done  on  rough-wall 
turbulent  boundary  layers.  For  a  more  in-depth  look  at  previous  studies  done  to  analyze 
the  flow  structure  in  sparsely  distributed  rough-wall  turbulent  boundary  layers  please 
refer  to  George  (2004).  The  current  subset  of  data  will  be  compared  with  the  data  and 
analysis  done  by  George  and  Simpson  (2000-04)  in  later  chapters. 

Past  studies  of  rough-wall  turbulent  boundary  layers  have  mainly  involved  the 
characterization  of  roughness  elements  as  equivalent  sand  grain  roughness.  This  idea 
was  stemmed  from  Nikuradse  (1932)  and  has  been  in  wide  use  ever  since.  A  significant 
amount  of  attention  has  also  been  devoted  to  the  effect  of  two-dimensional  roughness 
elements  (ribs,  grooves,  etc.).  Perry  et  al.  (1969),  Perry  et  al.  (1987),  and  Krogstad  and 
Antonia  (1994)  are  a  few  examples  found  in  previous  literature  related  to  two- 
dimensional  roughness  effects.  George  and  Simpson  have  taken  a  more  systematic 
approach  to  describe  rough-wall  turbulent  boundary  layers.  They  have  performed 
experiments  on  only  the  isolated  roughness  element  and  then  used  that  very  same  element 
to  construct  an  array  of  roughness  elements.  Information  gathered  from  these  efforts  has 
been  extensive  in  terms  of  the  mean  velocities,  turbulence  quantities,  and  skin  friction,  as 
well  as  the  related  fluid  dynamic  quantities  that  can  be  derived  from  the  above 
parameters.  The  current  research  is  directly  related  to  that  of  George  and  Simpson  in  that 
seven  isolated  roughness  elements  are  tested  as  well  as  a  fetch  of  one  of  those  isolated 
elements  (Gaussian  element).  Consequently,  with  the  addition  of  the  current  research, 
the  combined  data  base  will  encompass  a  wide  range  of  the  effects  that  each  roughness 
element  has  on  the  flow  field.  This  knowledge  in  turn  can  be  assembled  into  a  deeper 
understanding  of  how  the  physics  of  the  different  flows  are  connected. 
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1.4  Related  Modeling  and  Computational  Efforts 

Modeling  the  effects  of  surface  roughness  is  an  area  of  concern  in  many  practical 
engineering  applications.  Consequently,  there  has  been  a  lot  of  effort  in  recent  years  to 
calculate  flows  over  different  rough  surfaces  using  various  computational  fluid  dynamic 
(CFD)  models.  Many  current  roughness  models  to  this  point  have  involved  the  use  of 
empirical  ‘constants’  and  equivalent  sand  grain  roughness.  These  underdeveloped 
concepts  have  little  direct  relationship  to  realistic  roughness  and  cannot  predict  accurately 
and  consistently  the  flow  characteristics  for  different  roughness  shapes.  The  calculations 
of  various  roughness  flows  are  greatly  influence  by  the  choice  of  turbulence  model.  Most 
models  have  a  difficult  time  with  the  calculation  of  mean  flow  characteristics.  This 
difficulty  only  becomes  more  pronounced  when  turbulence  quantities  are  calculated.  The 
shortcomings  of  these  models  are  due  to  effects  not  taken  into  account  in  the  calculation 
methods.  Patel  (1998)  also  discusses  the  downfalls  of  current  CFD  models  and 
comments  that  high  Reynolds  number  flows  over  rough  surfaces  is  the  ‘Achilles  Heel’  of 
CFD.  Patel  suggests  the  need  for  further  experimental  research  concerning  the  effects  of 
roughness  in  high  Reynolds  number  turbulent  boundary  layers.  This  section  gives  some 
general  insight  as  to  what  the  objectives  of  some  models  are,  how  they  have  been  used  to 
calculate  different  roughness  flows,  and  the  results  of  such  efforts. 

Patel  (1998)  shows  the  experimental  data  of  Parthasarathy  and  Muste  (1994)  compared  to 
the  calculations  of  Yoon  et  al.  (1997)  using  the  standard  k-co  model  in  an  asymmetric 
channel  flow.  A  2  mm  thick  wire  mesh  and  a  series  of  rib-roughness,  8.9  mm  wide  x  13 
mm  tall  (50.2  mm  spacing),  were  tested  using  an  equivalent  sand  grain  roughness 
determined  from  empirical  correlations.  Mean  flow  characteristics  were  calculated  to 
within  no  less  than  a  7%  difference.  In  contrast,  the  Reynolds  shear  stress  was  not 
calculated  very  close  at  all.  Conclusions  from  this  experiment  were  that  the  standard  k-co 
model  does  not  capture  the  physics  of  the  flow  and  deteriorates  for  increasing  roughness 
size. 

A  force  field  approach  using  large  eddy  simulation  (LES)  is  discussed  in  Cui  et  al. 
(2003).  A  body-fitted  grid  is  avoided  by  applying  an  orthogonal  Cartesian  grid  to 
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complex  geometry  and  additional  stability  limits  are  not  introduced  by  the  computational 
scheme.  The  body  force  magnitude  is  determined  internally  during  numerical  solution, 
thus  making  it  a  more  general  model  than  those  requiring  empirical  inputs.  This  force 
field  LES  approach  was  compared  to  a  two-dimensional  wavy  rough  boundary  with  a 
mean  pressure  gradient  and  separation  (data  of  Buckles  et  al.).  In  this  case,  the  model 
calculates  premature  separation,  overestimates  reverse  flow,  and  calculates  late 
reattachment.  Although  this  model  does  have  its  downfalls,  it  does  show  improvements 
from  previous  models. 

A  modified  van  Driest  damping  function  model  was  introduced  by  Krogstad  (1991). 
This  model  simulates  the  viscous  stresses  near  the  wall  by  manipulating  the  amount  of 
viscous  damping  applied.  Through  this  ‘manipulation’  process  various  parameters  are 
adjusted  to  match  experimental  results.  This  model  was  compared  with  data  from  a 
three-dimensional  wing  body  junction  flow  using  24-grit  sandpaper  to  roughen  the 
surface.  The  model  was  tested  on  two  different  streamlines  within  the  flow:  line  1  is  a 
line  of  symmetry  that  leads  up  to  an  ordinary  separation  point,  and  line  3  passes  outside 
the  separated  region.  Line  2  leads  up  to  the  horseshoe  vortex  where  three-dimensional 
separation  exists  and  was  not  calculated.  The  model  predicts  early  separation  on  line  1 
due  to  the  fact  that  no  viscous-inviscid  interaction  is  accounted  for.  The  model  does 
calculate  line  3  better  due  to  the  smaller  pressure  gradients  and  less  interaction  seen  along 
this  line. 

Shim  et  al.  (2000  and  2001)  compare  four  different  turbulent  models  used  to  calculate  the 
flow  field,  namely  the  lift  and  drag  coefficients,  involving  iced  airfoils  at  various  angles 
of  attack  (data  from  Addy  et  al.).  Two  different  elemental  flow  fields  are  examined.  A 
somewhat  smooth  Gaussian-shaped  element  (rime  ice)  and  a  similarly  shaped  element 
with  prominent  horns  near  the  peak  of  the  element  (glaze  ice).  In  these  flow  regimes  two 
1 -equation  models  are  utilized  and  two  2-equation  models  are  used.  The  1-eqn.  models 
are  the  Baldwin-Barth  (BB)  model  and  the  Spalart  and  Allmaras  (SA)  model.  The  BB 
model  is  derived  from  the  k-e  model  and  its  objective  is  to  accurately  calculate  the 
turbulence  phenomena  in  separated  regions  and  shear  layers  which  cannot  be  handled  by 
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algebraic  models.  The  SA  model  has  the  same  objective  as  the  BB  model  but  is  not 
derived  from  a  2-eqn.  model.  The  SA  model  includes  a  destruction  term  that  depends  on 
the  distance  to  the  wall  and  is  forgiving  in  terms  of  near-wall  resolution  and  stiffness. 
This  model  also  provides  a  smooth  laminar  to  turbulent  transition  but  only  at  points 
specified  by  the  user.  The  SA  model  is  less  grid  sensitive  than  the  BB  model  and  more 
accurate  in  calculating  boundary  layer  profiles,  skin  friction  and  other  wake  properties. 

The  2-eqn.  models  are  the  Chien  k  -  e  (k  -  e)  model  and  the  shear  stress  transport  (SST) 
model.  The  k  -  e  model  includes  terms  describing  the  behavior  of  the  turbulent  shear 
stress,  turbulent  kinetic  energy  (TKE),  and  rate  of  dissipation  near  a  solid  wall  using  a 
Taylor  Series  expansion.  This  model  retains  a  proper  physical  behavior  of  the  balance 
between  dissipation  and  molecular  diffusion  of  the  TKE.  The  SST  model  is  developed 
from  the  original  k-co  model  of  Wilcox.  The  SST  is  identical  to  the  k-w  model  in  the 
inner  region  of  flow  and  changes  gradually  to  the  standard  k-e  model  toward  the  outer 
edge  of  the  boundary  layer.  As  compared  with  the  original  k-co  model,  the  SST  model 
has  the  ability  to  account  for  the  transport  of  the  turbulent  shear  stress  in  boundary  layers 
with  adverse  pressure  gradients. 

All  models  were  found  to  have  good  results  at  low  angles  of  attack  but  the  k-e  model 
did  not  converge  at  high  angles  of  attack.  The  SST  model  was  determined  to  give  the 
best  calculation  for  the  rime  ice  case.  The  maximum  percent  difference  for  the  lift 
coefficient  was  =13%  whereas  the  maximum  percent  difference  for  the  drag  coefficient 
was  =9%.  However,  the  SA  model  gave  the  best  calculations  for  the  glaze  ice  flow 
regime.  The  maximum  percent  difference  for  the  lift  coefficient  was  —7%  whereas  the 
maximum  percent  difference  for  the  drag  coefficient  was  =8%.  The  deviations  of  the 
models  at  higher  angles  of  attack  are  attributed  to  the  formation  and  shedding  of  large 
separation  bubbles  and  the  resultant  unsteadiness  within  the  flow. 
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1.5  Organization  of  Thesis 

This  thesis  is  comprised  of  five  main  chapters  and  six  appendices.  Each  chapter 
discusses  the  most  significant  aspects  of  the  research  whereas  the  appendices  contain 
supporting  material  related  to  the  main  discussion.  Chapter  2  discusses  the  apparatus  and 
instrumentation  employed  to  conduct  the  research.  Information  contained  in  this  chapter 
is  related  to  the  wind  tunnel  setup,  LDV  measurement  system,  data  post-processing, 
quality  of  flow  conditions,  description  of  the  tested  roughness  elements,  and  LDV  probe 
alignment  techniques.  Chapter  3  contains  the  results  of  the  research  conducted.  All 
mean  velocities,  turbulence  quantities,  skin  friction,  and  other  quantities  that  are  derived 
from  the  previous  parameters  can  be  found  in  this  chapter.  Oil  flow  visualizations  are 
also  described  within  Chapter  3.  Chapter  4  provides  the  related  discussion  to  the 
experimental  results  as  well  as  data  correlations  and  flow  field  schematics.  These 
schematics  give  a  description  of  the  flow  field  using  the  data  gathered  from  Chapter  3. 
Chapter  5  provides  the  conclusions  and  results  gathered  from  the  previous  four  chapters. 
The  references  used  in  this  study  can  be  found  immediately  following  Chapter  5. 

Appendices  A  through  F  can  be  found  at  the  end  of  the  thesis.  Appendix  A  contains 
information  related  to  various  trip  arrangements  and  the  generation  of  the  desired 
boundary  layer  characteristics.  Appendix  B  is  the  derivation  of  the  drag  equation  used  in 
the  single  roughness  element  analysis.  Appendix  C  contains  the  uncertainties  that  are 
pertinent  to  the  conducted  research.  Appendix  D  includes  a  brief  analysis  of  the  data 
related  to  the  Gaussian  fetch  of  roughness.  A  discussion  related  to  the  flow  angle  data 
can  be  found  in  Appendix  E.  Finally,  Appendix  F  gives  the  y-z  coordinates  that  define 
the  shape  of  both  the  large  and  fine  grooved  elements. 


Chapter  2  Apparatus  and  Instrumentation 


The  purpose  of  this  chapter  is  to  familiarize  the  reader  with  the  experimental  apparatus 
and  techniques  used  in  conducting  the  present  research  in  the  Department  of  Aerospace 
and  Ocean  Engineering  at  Virginia  Polytechnic  Institute  and  State  University.  Virginia 
Tech’s  small  boundary  layer  wind  tunnel  was  resurrected  and  many  modifications  were 
performed  in  order  to  conduct  the  present  research.  Pitot  tube  and  three-velocity 
component  laser-Doppler  velocimetry  (LDV)  measurements  were  taken  and  used  to 
determine  the  flow  characteristics  of  the  wind  tunnel  following  setup.  All  LDV  data 
obtained  for  the  various  roughness  cases  was  obtained  in  the  same  fashion  and  are 
discussed  below.  Other  details  that  will  be  discussed  include  the  following:  wind  tunnel, 
LDV  system,  flow  conditions,  tested  roughness  elements  and  their  respective 
manufacturing  techniques,  and  data  acquisition  and  post  processing. 

2.1  Small  Boundary  Layer  Wind  Tunnel 

The  function  of  the  wind  tunnel  is  to  deliver  low  turbulence  air  at  constant  temperature 
and  velocity  to  the  test  section.  The  wind  tunnel  is  a  closed-circuit  design  and  can 
accommodate  the  necessary  seeding  of  the  flow  for  LDV  measurements.  The  small 
boundary  layer  wind  tunnel  previously  located  in  Lab  #7  in  the  basement  of  Randolph 
Hall  is  now  located  in  Lab  #6.  In  order  to  accommodate  the  ever  increasing  need  for 
space,  many  modifications  were  made  to  the  previous  wind  tunnel  setup.  For  information 
concerning  the  previous  setup  see  E.  J.  Smith  et  al.  (1990),  as  well  as  K.  R.  Saripalli  and 
R.  L.  Simpson  (1980).  The  result  is  the  most  space  efficient  and  versatile  design  for  the 
small  boundary  layer  wind  tunnel.  The  dimensions  of  the  wind  tunnel  setup  are 
approximately;  6.25  meters  long,  1.73  meters  wide  and  2.20  meters  high.  There  are 
seven  main  parts  to  the  wind  tunnel  system:  air  conditioning  system,  speed  control  valve 
and  filter  box,  blower  system,  plenum  chamber,  contraction  section,  test  section,  and  the 
return  ducting.  Figures  of  the  following  descriptions  can  be  found  at  the  end  of  the 
chapter.  Figure  2.1  gives  a  schematic  of  the  wind  tunnel  arrangement  while  Figures  2.2 
and  2.3  show  different  views  of  the  actual  setup  of  the  wind  tunnel.  The  modifications 
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will  be  discussed  below  along  with  a  description  of  the  wind  tunnel  setup  in  the  order  the 
flow  travels  through  the  system. 

2.1.1  Air  Conditioning  System 

The  first  component  that  the  flow  encounters  in  the  experimental  setup  is  the  air 
conditioner  (AC).  The  flow  travels  through  30.48  cm  diameter  ducting  that  is  attached  to 
the  AC.  In  order  to  accommodate  the  duct  work,  the  AC  rests  on  top  of  a  7.62  cm  by 
3.81  cm  aluminum  C-channel  support  bracket  bolted  to  the  top  of  the  plenum  chamber, 
see  Section  2.2.4.  The  function  of  the  AC  was  to  maintain  a  constant  temperature  in  the 
wind  tunnel  test  section.  The  AC  utilized  is  a  model  ACM244XL0  manufactured  by 
Whirlpool  with  a  1  phase  60  Hz  motor  and  has  a  capacity  of  25.3  MJ  (24,000  BTU’s). 
By  controlling  only  the  fan  and  cooling  adjustments  on  the  AC,  it  was  possible  to  obtain 
an  equilibrium  temperature  condition  in  the  test  section  of  ±0.56°C  at  25°C  (±1°F  at 
77°F).  It  is  important  to  note  that  the  AC  condenser  never  ceased  running  in  order  to 
maintain  a  constant  temperature.  If  such  a  condition  did  happen,  and  only  the  fan  was 
left  running  for  a  period  of  time,  there  would  be  an  unwanted  temperature  fluctuation 
within  the  system. 

A  series  of  ductwork  was  constructed  and  added  to  the  AC  in  order  to  expel  unwanted  hot 
air  out  of  the  laboratory.  Sheet  metal  was  cut  and  subsequently  fitted  and  sealed  to  the 
AC  exit  and  then  connected  to  a  rectangular  duct  system  that  measured  41.91  cm  by 
21.59  cm  and  was  185.42  cm  tall.  This  ducting  then  connected  to  the  wall  which  led  to 
an  outlet  for  the  unwanted  hot  air.  The  normal  room  AC  and  heating  system  was  utilized 
in  order  to  keep  the  temperature  of  the  laboratory  between  22.22°C  and  24.44°C  (72°F  to 
76°F).  With  this  accomplished  it  was  possible  to  maintain  the  ±0.56°C  at  25°C  (±1°F  at 
77°F)  condition  within  the  test  section. 

2.1.2  Speed  Control  Valve  and  Filter  Box 

The  next  component  in  the  flow  circuit  is  the  speed  control  valve.  This  valve  is  a 
butterfly  valve  contained  inside  the  30.48  cm  diameter  ducting.  It  is  possible  to  limit  the 
amount  of  air  ingested  by  the  blower  by  rotating  the  valve  in  order  to  obtain  a  constant 
speed  in  the  test  section.  With  the  valve  completely  open,  the  wind  tunnel  can  achieve  a 
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test  section  velocity  of  approximately  44  m/s.  Following  the  butterfly  valve  the  flow 
passes  into  the  filter  box.  The  filter  box  is  constructed  out  of  1.27  cm  thick  plywood  on 
all  sides  except  for  the  side  that  it  connects  to  the  blower  intake.  This  side  is  constructed 
of  0.64  cm  thick  aluminum  in  order  to  bolt  the  filter  box  to  the  blower  intake.  A  0.32  cm 
thick  rubber  gasket  lies  between  the  filter  box  and  the  intake  flange  of  the  blower  in  order 
to  reduce  the  transmission  of  vibration  throughout  the  system.  The  dimensions  of  the  box 
are  as  follows;  68.90  cm  tall,  64.77  cm  wide,  and  43.82  cm  deep.  The  ducting  connects 
to  the  filter  box  via  a  sheet  metal  connection.  A  fine  wire  mesh  screen  is  fastened  to  the 
inside  top  of  the  filter  box  in  order  to  ensure  that  there  will  be  no  debris  taken  into  the 
blower. 

2.1.3  Blower  System 

Once  the  air  passes  through  the  filter  box  it  enters  the  suction  side  of  the  blower.  A  new 
blower  was  selected  in  order  to  provide  the  most  versatility  to  the  small  wind  tunnel 
setup.  The  final  blower  selection  made  was  a  New  York  Blower  Size  194  Gl-DH  Series 
20  Fan,  Arrangement  10  with  a  DH  wheel  and  clockwise  rotation.  This  general 
industrial,  centrifugal  blower  is  driven  by  a  constant  speed  v-belt  drive,  5.59  kW  (7.5 
horsepower),  3500  rpm,  230  Volt  motor  which  is  covered  for  safety  precautions.  This 
blower  is  capable  of  delivering  62.30  m3/min  (2200  cfm)  of  air  at  30.48  cm  H2O  static 
pressure  at  2529  rpm,  5.16kW  (6.92  Bhp),  and  21.1 1°C  (70°F). 

2.1.4  Plenum  Chamber 

As  the  air  exits  the  blower  it  enters  the  plenum  chamber.  The  purpose  of  the  plenum  is 
twofold.  First,  the  air  ‘settles’  and  large  scale  fluctuations  produced  by  the  blower  are 
eliminated.  Second,  the  acoustic  noise  is  absorbed  and  the  flow  is  quieted.  The  plenum 
is  a  rectangular  box  that  is  constructed  out  of  1.91  cm  thick  plywood  with  inside 
dimensions  of  1 16.91x1 16.91x157.48  cm.  The  inlet  and  outlet  dimensions  of  the  plenum 
are  50.80x50.80  cm  and  58.42x58.42  cm  respectively.  The  plenum  is  symmetric  about 
the  centerline  that  extends  from  the  inlet  to  the  exit  of  the  plenum.  The  plywood  is 
supported  via  0.318  cm  thick  3.81x3.81  cm  steel  angle  which  is  bolted  to  the  plywood 
sides.  The  plenum  is  connected  to  the  blower  outlet  flange  and  the  contraction  via  0.3 1 8 
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cm  thick  rubber  in  order  to  eliminate  any  mechanical  vibrations  to  the  test  section.  The 
entire  plenum  chamber  is  sealed  with  silicone  rubber  in  order  to  prevent  any  leaks. 

To  facilitate  the  ‘mixing’  of  the  flow,  there  is  a  baffle  plate  placed  inside  the  plenum 
chamber,  see  Figure  2.4.  The  baffle  plate  is  63.50x63.50  cm  constructed  from  1.91  cm 
thick  plywood  and  secured  to  the  top  and  bottom  walls  of  the  plenum  by  means  of  the 
previously  discussed  steel  angle.  It  is  positioned  63.50  cm  aft  of  the  plenum  exit  and 
centered  between  the  plenum  side  walls.  Directly  behind  the  downstream  side  of  the 
baffle  plate,  a  1 .27  cm  diameter  copper  tube  extends  to  the  center  of  the  baffle  plate,  see 
Figure  2.5.  This  tube  introduces  the  ‘seed’  particles  to  the  flow  to  ensure  the  air  will  get 
uniformly  mixed  with  ‘seed’  before  the  flow  reaches  the  contraction  and  test  section. 
The  seeding  particles  used  in  the  current  research  are  di-octyl  phthalate  (DOP).  The 
liquid  DOP  is  converted  to  nearly  constant  diameter  particles,  having  a  Gaussian 
distribution  centered  around  1  pm  in  size,  via  an  aerosol  generator  system.  Once  the 
DOP  is  in  its  aerosol  form  it  is  introduced  into  the  plenum  chamber.  In  order  to  catch 
excess  liquid  DOP,  a  20.32x20.32x6.35  cm  stainless  steel  drip  pan  is  located  directly 
under  the  copper  tubing  on  the  bottom  wall  of  the  plenum.  A  ball  valve  is  attached  to  the 
drip  pan  on  the  outside  of  the  plenum  that  enables  one  to  release  the  contents  of  the  pan. 

A  7.62  cm  thick  layer  of  polyester  urethane  foam  with  a  standard  anechoic  wedge  surface 
geometry  (Smith,  et  al.  1990)  covers  the  inlet,  top,  bottom,  and  side  walls  of  the  plenum 
as  well  as  the  upstream  side  of  the  baffle  plate.  This  foam  helps  to  absorb  the  acoustic 
noise  generated  by  the  blower.  For  more  information  on  the  foam  selection  see  Smith  et 
al.  1990.  The  final  effective  inside  dimensions  of  the  plenum  chamber  are 
101.60x101.60x149.86  cm. 

2.1.5  Contraction  Section 

Following  the  exit  of  the  plenum  chamber,  the  flow  enters  the  contraction  section.  The 
contraction  section  consists  of  two  different  contractions  as  well  as  a  rectangular  section 
that  houses  three  screens  and  a  honeycomb  ( Saripalli  and  Simpson,  1980).  The  purpose 
of  the  contraction  is  to  connect  the  plenum  to  the  test  section  as  well  as  increase  the  flow 
velocity  and  reduce  the  turbulence  levels  of  the  flow.  The  contraction  section  is 
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constructed  from  0.71  mm  thick  sheet  metal.  The  first  contraction  that  leads  to  the 
rectangular  section  is  45.72  cm  in  length  and  has  inlet  and  outlet  dimensions  as  follows, 
53.34x53.34  cm  and  38.10x22.86  cm.  Once  the  flow  has  entered  the  rectangular  section 
it  encounters  the  honeycomb  which  serves  the  purpose  of  eliminating  large  scale 
turbulence  within  the  flow.  After  passing  through  the  honeycomb,  the  flow  reaches  the 
three  screens.  These  screens  are  located  13.97  cm  apart  with  the  purpose  of  letting  the 
wakes  from  each  of  the  previous  screens  die  out  before  reaching  the  following  screen. 
The  main  function  of  the  screens  is  to  reduce  the  turbulence  level  of  the  flow.  The 
second  and  final  contraction  has  an  inlet  dimension  of  40.64x24.13  cm  and  an  outlet 
dimension  of  10.16x24.13  cm. 

2.1.6  Test  Section 

Once  the  air  exits  the  final  contraction,  it  enters  the  test  section.  The  test  section  and 
contraction  are  lined  up  to  facilitate  a  smooth  flow  transition.  Diagrams  of  the  test 
section  within  the  wind  tunnel  setup  can  be  seen  in  Figures  2.1  and  2.2.  The  test  section 
has  an  overall  length  of  199.39  cm  and  is  24.13  cm  wide.  The  test  section  floor  is 
constructed  of  hard  ‘fin-form’  plywood  with  a  polished  surface.  There  are  two  openings 
in  the  test  section  floor,  in  which  inserts  can  be  placed,  that  make  it  possible  for  data 
acquisition.  The  opening  nearest  the  contraction  measures  40.64  cm  by  23.50  cm  and  has 
a  0.48  cm  thick  smooth  aluminum  plate  inserted  into  the  opening.  The  aluminum  plate  is 
supported  via  three  evenly  spaced  sections  of  wood  laid  perpendicular  to  the  tunnel  floor. 
Each  piece  of  wood  is  routered  on  their  respective  ends  to  provide  the  necessary  support 
and  a  smooth  continuation  of  the  test  section  floor.  The  second  opening  in  the  floor 
measures  71.12  cm  by  23.50  cm  and  has  0.79  cm  thick  Plexiglas  inserted  into  the 
opening.  This  Plexiglas  insert  is  supported  in  a  similar  fashion  as  discussed  previously. 
The  only  difference  being  that  four  wooden  sections  were  utilized  for  support,  one  on 
each  edge,  in  order  to  allow  the  necessary  LDV  access  needed  for  data  acquisition.  There 
is  a  20.32  cm  diameter  hole  cut  into  the  Plexiglas  in  which  a  20.32  cm  diameter,  0.64  cm 
thick  piece  of  float  glass  can  be  placed.  It  is  through  this  piece  of  float  glass  where  the 
LDV  measurements  are  taken.  Actual  locations  of  measurements  will  be  discussed  in 
further  detail  in  the  following  chapters.  It  is  important  to  note  however,  that  the  center  of 
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the  glass  disc  is  located  1 16.72  cm  from  the  contraction  exit.  The  side  walls  are  20.96  cm 
high  and  are  constructed  out  of  0.64  cm  thick  glass. 


The  test  section  ceiling  has  three  main  sections.  There  are  two  removable/adjustable 
sections  and  a  permanent  section,  all  of  which  are  constructed  from  0.95  cm  thick 
Plexiglas.  The  permanent  section  is  located  the  last  13.97  cm  of  the  test  section  and  has  a 
height  of  10.48  cm.  The  rest  of  the  ceiling  is  split  into  two  parts  in  order  to  facilitate  the 
removal  and  replacement  of  the  trip  arrangement  discussed  previously.  It  is  important  to 
note  that  the  two  removable  ceiling  sections  are  reinforced  with  2.54  cm  by  1.27  cm 
aluminum  c-channel  and  are  secured  together  in  a  manner  which  permits  a  smooth 
transition  from  one  section  to  the  other.  The  forward  most  section  has  a  length  of  5 1 .44 
cm  and  is  the  section  that  can  be  removed  for  the  replacement  of  the  sandpaper,  etc.  The 
rest  of  the  removable  ceiling  has  a  length  of  133.98  cm  and  transitions  smoothly  to  the 
permanent  ceiling.  There  are  twenty-four  0.95  cm  holes  drilled  along  the  centerline  of 
the  removable  ceiling  in  order  to  provide  the  necessary  pitot  probe  access  to  obtain  a 
zero-pressure  gradient  flow.  The  first  hole  is  22.86  cm  from  the  contraction  exit.  All 
subsequent  holes  are  spaced  6.35  cm  apart,  except  between  holes  ten  and  eleven.  The 
spacing  between  these  two  holes  is  12.70  cm.  This  is  due  to  the  fact  that  a  support  beam 
interferes  with  the  regular  hole  placement.  At  hole  number  twenty,  there  are  two  more 
holes  drilled  3.81  cm  on  either  side  of  the  centerline  hole.  These  holes  give  access  to  the 
tunnel  in  order  to  probe  the  boundary  layer  and  determine  the  two-dimensionality 
characteristics  of  the  flow.  See  Table  2.1  below  for  port  #’s  and  their  related  locations 
relative  to  the  contraction  exit. 

Table  2.1 :  Port  #’s  and  related  distance  from  contraction  exit 


Dist.  from  Contraction 
Exit 

Port# 

Inches 

Centimeters 

1 

9 

22.86 

2 

11.5 

29.21 

3 

14.0 

35.56 

4 

16.5 

41.91 

5 

19.0 

48.26 

6 

21.5 

54.61 

7 

24.0 

60.96 

8 

26.5 

67.31 
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9 

29.0 

73.66 

10 

31.5 

80.01 

11 

36.5 

92.71 

12 

39.0 

99.06 

13 

41.5 

105.41 

14 

44.0 

111.76 

Center  of 
Glass  Disc 

45.95 

116.72 

15 

46.5 

118.11 

16 

49.0 

124.46 

17 

51.5 

130.81 

18 

54.0 

137.16 

19 

56.5 

143.51 

20 

59.0 

149.86 

21 

61.5 

156.21 

22 

64.0 

162.56 

23 

66.5 

168.91 

24 

69.0 

175.26 

There  is  a  single  static  pressure  port  148.59  cm  from  the  contraction  exit.  This  port  has  a 
diameter  of  0.05  cm  and  was  utilized  in  probing  the  boundary  layer.  The  ceiling  has  five 
aluminum  crossbars  that  allow  it  to  be  bolted  to  an  aluminum  c-channel,  attached  to  the 
glass  side  walls,  to  alleviate  unwanted  wandering  of  the  setup.  Two  bolts  located  at  each 
ceiling  crossbar  position  enable  one  to  adjust  the  ceiling  height.  The  ceiling  height  was 
adjusted  in  order  to  produce  a  zero  pressure  gradient  flow.  Actual  ceiling  height 
measurements  are  discussed  for  the  final  configurations  (single  element  and  Gaussian 
fetch)  in  Section  2.4.2,  Table  2.2. 


It  is  important  to  note  that  the  blower,  contraction  section,  and  test  section  are  all  bolted 
to  the  floor  to  alleviate  any  wandering  of  the  wind  tunnel  system. 

2.1.7  Return  Ducting 

Once  the  flow  exits  the  test  section  it  enters  a  rectangular  to  circular  expansion  having 
inlet  dimension  of  24.13  cm  by  10.48  cm  and  outlet  dimensions  of  30.48  cm  in  diameter. 
The  expansion  takes  place  over  a  length  of  26.67  cm.  This  expansion  then  connects  with 
the  rest  of  the  circular  ducting  which  is  also  30.48  cm  in  diameter.  The  ducting  utilized  is 
a  general  commercial  type  of  ducting.  In  order  to  assist  the  space  efficiency  of  the  tunnel 
setup,  the  ducting  was  fashioned  in  a  vertical  return  loop.  For  safety  precautions,  as  well 
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as  ease  of  use,  the  bottom  of  the  ducting  is  182.88  cm  off  of  the  floor  and  is  located 
directly  above  the  test  section.  A  0.32  cm  thick  section  of  rubber  acts  as  a  joint  between 
two  pieces  of  ducting  within  the  return  loop  help  to  alleviate  any  mechanical  vibrations  to 
the  test  section.  The  return  system  is  supported  via  three  1.91  cm  by  1.91  cm  steel  angle 
brackets.  Two  brackets  are  directly  on  either  side  of  the  rubber  connection  and  the  other 
bracket  can  be  found  attached  to  a  steel  angle  located  on  the  front  end  of  the  plenum 
chamber.  The  return  ducting  is  attached  to  the  AC  and  then  continues  into  the  filter  box. 
All  connections  between  ductwork  and  other  various  parts  of  the  wind  tunnel  system  are 
sealed  with  silicon  rubber  in  order  to  prevent  leaks. 


2.2  LDV  System 

The  LDV  system  utilized  consists  of  a  miniature  LDV  head,  optical  table,  a  signal 
conditioning  arrangement,  3  digital  signal  processors,  traverse  system,  and  a  particle 
seeding  system.  This  setup  allows  the  coincident  and  instantaneous  measurement  of  the 
U,  V,  and  W  components  of  velocity.  Consequently,  all  mean  velocities,  Reynolds 
stresses,  and  triple  products  can  be  gathered  from  this  measurement  scheme.  This 
measurement  system  is  also  described  by  Chesnakas  and  Simpson  (1994)  as  well  as  Tian 
(2003)  and  Tang  (2004).  The  LDV  system  has  a  nearly  spherical  measurement  volume 
of  ~50  pm  in  diameter  and  fringe  spacing  equal  to  ~5  pm.  These  parameters  were 
calculated  using  equations  from  Durst  et.  al  (1995)  and  Durst  et.  al  (1981)  respectively. 
Very  near-wall  measurements  can  be  obtained  with  this  LDV  system.  Measurements 
with  this  system  closely  agree  with  direct  numerical  simulation  (DNS)  data  for  a  two- 
dimensional  turbulent  boundary  layer  which  gives  added  confidence  in  the  obtained  data. 
An  uncertainty  analysis  for  the  LDV  system  and  results  can  be  found  in  Appendix  C. 
Kuhl  and  Simpson  (2001)  have  investigated  the  velocity  bias  effects  and  transient  time 
broadening  effects  of  the  measurement  system  and  determined  them  to  be  negligible. 
Likewise,  Olcmen  et.  al  (1998)  investigated  the  instrument  broadening  effects  on  the 
measured  frequency  by  the  digital  signal  processors  (Macrodyne  FDP3 1 00)  and  they  too 
were  determined  to  be  negligible. 
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2.2.1  Optical  Table 

A  Coherent  argon-ion  laser  powers  the  LDV  measurement  system.  The  laser  is  encased 
in  a  Plexiglas  box  that  has  a  small  air  conditioning  unit  placed  inside.  Both  the  box  and 
the  air  conditioner  serve  to  keep  the  laser  in  a  clean  dust-free  environment.  A  view  of  the 
optical  table  setup  near  the  laser  exit  can  be  seen  in  Figure  2.6,  whereas  a  full  view  of  the 
optical  table  is  shown  in  Figure  2.7.  Following  the  exit  of  the  main  beam  out  of  the  laser, 
the  beam  is  directed  by  a  series  of  mirrors  into  a  prism.  The  main  beam  is  a 
conglomeration  of  various  wavelengths.  Thus,  the  prism  separates  the  colors  and 
associated  wavelengths  in  order  for  the  desired  green  (514.5  run)  and  blue  (488  nra) 
beams  to  be  utilized.  Once  the  colors  are  separated  the  beams  are  allowed  to  travel  a 
sufficient  distance  in  order  for  them  to  spread  out  from  one  another.  The  green  beam  is 
then  directed  by  a  mirror  into  a  polarization  rotator.  This  rotator  enables  the  manual 
shifting  of  beam  power  (vertically  or  horizontally  polarized  laser  light)  in  order  for 
optimum  signal  quality.  From  the  polarization  rotator  the  green  beam  passes  through  a 
beam  splitting  cube.  This  cube  splits  the  vertically  and  horizontally  polarized  beams  and 
directs  them  90°  away  from  one  another.  After  the  green  beams  are  split  they  travel  into 
their  respective  Bragg  Cells  which  shift  the  beams  by  a  known  frequency  amount,  which 
in  this  case  is  -27  MHz,  0  MHz,  and  50  MHz.  The  purpose  of  the  Bragg  Cell  is  to  enable 
one  know  the  direction  in  which  the  fringes  are  moving.  This  in  turn  will  eliminate  any 
directional  ambiguity  that  may  arise  within  the  measured  Doppler  signal.  It  is 
worthwhile  to  note  that,  for  this  system,  the  maximum  data  rate  was  obtained  when  more 
power  was  in  the  unshifted  0  MHz  green  beam  as  opposed  to  an  equal  distribution  of 
power.  The  blue  beam  requires  no  beam  splitting  cube  and  is  split  into  two  equal 
intensity  beams  via  the  40  MHz  Bragg  Cell  driver.  An  unshifted  beam  of  0  MHz  and  a 
shifted  beam  of  40  MHz  will  come  from  this  Bragg  Cell.  In  order  to  have  a  clear  signal 
from  each  of  the  3  beam  pairs,  it  is  necessary  to  shift  the  beams  by  a  frequency  amount 
large  enough  that  will  eliminate  cross-talk  between  the  signals.  All  five  beams,  -27  MHz, 
0  MHz,  50  MHz  (green)  and  0  MHz,  40  MHz  (blue),  are  then  directed  into  their 
respective  beam  launchers  where  a  microscope  objective  focuses  the  beam  into  the  fiber 
optics.  From  there  the  beams  travel  through  polarization  preserving  optical  fibers  and 
into  the  LDV  head  itself. 
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2.2.2  LDV  Head 

The  LDV  head  is  comprised  of  two  separate  aluminum  beam  assemblies  that  house  the 
necessary  optics  to  transmit  the  3  green  beams  (514.5  nm)  and  2  blue  beams  (488  nm). 
The  beam  assemblies  are  at  a  right  angle  with  respect  to  one  another.  Within  these 
assemblies  are  beam  collimating  lenses  as  well  as  focusing  lenses  which  make  it  possible 
to  obtain  a  constant  diameter  laser  beam  that  can  be  aligned  at  the  focal  point  of  the  lens. 
Figure  2.8  shows  a  close-up  view  of  the  head  where  Figure  2.9  shows  the  LDV  head 
while  acquiring  data.  The  off-axis  backscatter  signal  is  received  through  an  array  of 
receiving  optics  that  is  also  located  in  the  center  of  the  LDV  head.  Both  the  green  beam 
assembly  as  well  as  the  receiving  optics  assembly  are  adjustable  which  aids  in  the  beam 
alignment  process. 

The  LDV  head  is  attached  to  a  base  which  in  turn  is  attached  to  a  three-dimensional 
traverse  system.  The  current  system  is  on  the  same  traverse  base  as  a  comprehensive 
LDV  system  being  developed  by  Lowe  (2004).  This  enables  the  full  utilization  of  both 
systems  without  constant  realignment  issues.  This  traverse  system  allows  full  mobility  of 
the  LDV  head  in  the  flow  field.  The  vertical  (y)  traversing  is  accomplished  via  a 
National  Aperture  Inc.  MC-4B  Series  (Servo  3000)  controller,  where  the  spanwise  (z) 
and  streamwise  (x)  directions  are  controlled  with  a  Velmex  VP9000  Series  controller. 
For  additional  information  on  the  x-z  traverse  system  see  Lowe  (2004).  The  entire 
traverse  system  and  the  LDV  head  were  situated  underneath  the  tunnel  floor  so  as  to  not 
interfere  with  the  flow  field  in  any  way.  Measurements  could  be  made  through  a  20.32 
cm  diameter,  0.64  cm  thick  piece  of  float  glass  inserted  into  the  wind  tunnel  floor.  The 
center  of  this  glass  disc  is  1 16.72  cm  from  the  contraction  exit. 

2.2.3  LDV  Seeding 

An  aerosol  generator  designed  by  Echols  and  Young  (1963)  was  used  to  seed  the  wind 
tunnel.  This  generator  used  dioctal  phthalate  as  the  seeding  material.  The  mean  particle 
size  is  ~1  pm  in  diameter.  The  air  pressure  going  into  the  generator  was  kept  near  1 1  psi 
in  order  to  maintain  optimum  signal  quality.  The  seeding  material  was  introduced  into 
the  wind  tunnel  via  a  0.5  inch  (1.27  cm)  diameter  copper  pipe.  The  air/seed  mixture 
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exhausted  from  the  air  conditioner  was  sent  through  an  Aprilaire  Electronic  Air  Cleaner 
before  being  exhausted  outside. 

2.2.4  Signal  Processing 

With  all  five  laser  beams  crossing  at  the  receiving  lens  focal  point,  instantaneous 
coincident  three-velocity  component  data  is  sent  to  the  signal  setup.  The  signal  receiving 
setup,  stationed  on  the  optical  table,  gathered  the  Doppler  bursts  from  the  multimode 
receiving  optical  fiber.  Two  photo-multiplier  (PM)  tubes  from  Electron  Tubes  Limited, 
models  9125B  and  9124B,  received  the  separated  green  and  blue  light  respectively.  The 
signal  from  the  photo-multiplier  tubes  was  then  transmitted  into  an  electrical  voltage 
signal  and  sent  to  the  signal  conditioning  electronics. 

All  of  the  electronics  are  located  on  a  single  circuit  board.  The  green  signal  coming  out 
of  its  respective  PM  tube  is  divided  into  two  signals  related  to  the  50  MHZ  and  -27  MHz 
shifted  beams.  These  two  signals  are  then  taken  along  with  the  blue  40  MHz  signal  and 
amplified  as  well  as  mixed  with  their  respective  RF  signals.  These  RF  signals  are 
produced  by  three  variable  RF  generators  that  allow  the  extraction  of  the  meaningful 
Doppler  signals.  These  mixed  signals  are  then  transferred  to  their  respective  Macrodyne 
FDP  3100  frequency  processors.  The  three  Macrodynes  were  set  with  a  coincidence 
window  equal  to  10  ps  in  order  to  obtain  coincident  data  from  the  three  signals.  A 
validation  ratio  from  each  of  the  three  Macrodynes  was  always  within  a  range  of  97%  to 
99%.  The  acquired  number  of  samples  per  point  varied  between  15,000  near  the  wall  to 
30,000  away  from  the  wall.  The  Macrodyne  signals  were  sent  to  an  IBM  386/33C  PC 
with  a  Dostek  1400A  Laser  Velocimeter  Interface,  and  a  TCEM  daughter  board  option, 
via  three  25-pin  parallel  cables. 


2.3  Post  Processing 

Due  to  the  significant  amount  of  electronics  used  to  measure  the  LDV  signals  there  is  an 
inherent  amount  of  noise  that  must  be  taken  into  account  when  processing  the  LDV  data. 
The  method  utilized  to  remove  the  irrelevant  noise  from  the  data  was  that  same  as  that 
used  by  Olcmen  and  Simpson  (1995).  A  parabola  was  fit  to  each  side  of  the  logarithm  of 
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the  velocity  component  histogram  ordinate  over  a  range  of  1%  to  80%  of  the  peak 
histogram  value.  Histogram  values  above  and  on  the  outskirts  of  the  parabolic  fit  to  the 
histogram  were  considered  noise  and  discarded.  If  one  velocity  component  was  deemed 
unacceptable  then  all  three  velocity  components  were  discarded.  Following  the  initial 
noise  reduction  the  data  were  rotated  to  match  the  tunnel  coordinate  system.  The  same 
process  was  again  followed  and  the  rotated  histograms  were  fit  with  parabolas  and  the 
noise  discarded. 

Determining  the  position  of  the  wall  relative  to  the  measurement  volume  is  critical  to 
obtaining  meaningful  results.  When  the  measurement  volume  is  positioned  on  the  wall  a 
strong  signal  is  sent  from  the  spectrum  analyzer  to  the  oscilloscope.  This  signal  can  then 
be  utilized  to  find  the  position  of  the  wall.  Although  this  method  is  extremely  useful,  it  is 
difficult  to  consistently  determine  the  wall  position  to  less  than  ±30  pm  uncertainty. 
Factors  that  lead  to  this  discrepancy  are  associated  with  the  traverse  movements  as  well 
as  the  size  of  the  measurement  volume,  50  pm.  Consequently,  it  is  necessary  to  have  a 
more  refined  approach  to  finding  the  actual  wall  position.  This  is  done  by  applying  a 
least  square  fit  to  the  theoretical  mean  velocity  sublayer  equation,  as  used  by  Olcmen  and 
Simpson  (1995)  and  Kuhl  (2001), 

Q  =  C,y+C2y4  (2.1) 

where  Q  -  ^^U~2  +  W~2  ,  and  Cj  and  C2  are  coefficients  to  be  solved  for.  A  curve  is  fit 
through  Q  =  0  at  y  =  0  using  at  least  4  points  in  the  viscous  sublayer  y+<10.  This  fit  was 
done  in  an  iterative  fashion  by  changing  the  y-shift  value  in  order  to  maximize  the 
correlation  coefficient.  All  y-shifts  that  were  obtained  were  on  the  order  of  the 
measurement  volume  diameter,  ±50  pm,  which  is  the  uncertainty  in  finding  the  wall. 
Nearly  all  y-shift  values  were  in  the  range  of  20-30  pm.  Using  the  curve  fit  equation,  it  is 
possible  to  calculate  the  friction  velocity  via  the  following  relationship  with  the  wall 
shearing  stress p  =  v\dQldy\wall  =U2r ,  where  Cxv  =  U ].  When  the  skin  friction 

velocity  becomes  zero,  or  even  negative,  the  previous  method  of  finding  a  y-shift  value 
does  not  work.  Consequently,  in  the  near  element  measurement  regions  where  there  is 
separated  flow  and  backflow  an  average  y-shift  value  was  used.  As  was  discussed 
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previously,  all  y-shift  values  obtained  for  the  unseparated  profiles  were  between  20  and 
30  microns.  Thus,  an  average  value  of  25  microns  was  used  as  the  y-shift  value  for  the 
near  element  profiles. 

With  all  LDV  systems  it  is  impossible  to  achieve  perfect  alignment  of  the  LDV 
coordinate  system  with  the  wind  tunnel  coordinate  system.  Consequently,  it  is  important 
to  perform  a  series  of  final  rotations  on  the  data  sets.  A  slight  rotation  about  the  tunnel 
coordinate  system  x,  y,  and  z  axes  (roll,  yaw,  and  pitch  angles)  was  performed  on  the 
two-dimensional  mean  flow  data  in  order  for  it  to  reach  minimum  values  within  a  set  of 
tolerances.  Final  rotations  of  the  data  were  always  on  the  order  of  1°.  For  more 
information  on  these  rotations,  see  Kuhl  (2001). 


2.4  Quality  of  Flow  Conditions 
2.4.1  Spanwise  Two-Dimensionality 

After  the  initial  setup  of  the  small  boundary  layer  wind  tunnel  it  was  necessary  to 
determine  the  two-dimensionality  of  the  flow  field.  An  initial  series  of  pitot  tube 
measurements  were  taken  without  any  trip  arrangement  inserted  into  the  tunnel.  The  sole 
purpose  of  these  measurements  was  simply  to  indicate  that  there  was  a  constant  velocity 
core  within  the  test  section.  Thus,  the  tree-stream  speed  is  not  set  exactly  to  the  desired 
27.5  m/s  during  these  tests.  Pressure  readings  were  taken  with  a  calibrated  Series  475 
Mark  III  Digital  Manometer  by  Dwyer  with  an  uncertainty  of  ±0.01  inches  of  water. 
Figure  2.10  shows  spanwise  velocity  profiles  at  a  constant  height  of  2.23  inches  (5.66 
cm)  from  the  tunnel  floor  and  a  distance  of  50.38  inches  (127.95  cm)  from  the 
contraction  exit.  This  height  is  the  vertical  center  of  the  wind  tunnel  test  section. 
Measurements,  for  this  data  set  only ,  were  taken  starting  at  the  starboard  side  wall  (z  = 
0).  All  other  data  sets  presented  use  the  designated  axis  system  with  z  =  0  being  the 
centerline  of  the  wind  tunnel.  The  profiles  are  very  repeatable  and  show  a  constant 
velocity  core  of  about  12  cm.  The  profiles  tend  to  follow  a  smooth  curve  except  that 
some  edge  effects  can  be  seen  near  both  tunnel  walls.  LDV  measurements  are  taken  well 
inside  these  areas  and  therefore  these  regions  will  not  affect  the  data  in  any  way. 
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2.4.2  Trip  Arrangement  and  Zero  Pressure  Gradient 

Previous  research  directly  related  to  the  present  research  has  been  conducted  in  Virginia 
Tech’s  Boundary  Layer  Wind  Tunnel.  Consequently,  to  obtain  meaningful  results  that 
can  be  related  to  previous  work  it  is  imperative  that  the  flow  conditions  be  matched  in 
Virginia  Tech’s  Small  Boundary  Layer  Wind  Tunnel.  In  order  to  assist  the  production  of 
the  necessary  turbulent  boundary  layer  characteristics,  the  first  45.72  cm  of  the  test 
section  floor  were  modified.  A  0.32  cm  square  rod  was  placed  directly  at  the  contraction 
exit  followed  by  another  0.32  cm  square  rod  5.40  cm  downstream  of  the  contraction  exit. 
Silicon  carbide,  20  grit  sandpaper  with  a  width  of  20.32  cm,  manufactured  by  Norton, 
was  inserted  everywhere  on  the  first  45.72  cm  of  the  test  section  floor  not  occupied  by  a 
square  rod.  The  placement  of  the  second  0.32  cm  square  rod  was  determined  by  giving 
the  flow  sufficient  distance  to  reattach  before  applying  another  step.  Thus,  with  the 
discussed  arrangement,  a  thicker  boundary  layer  with  a  greater  momentum  deficit  was 
created  in  order  to  obtain  the  necessary  turbulent  boundary  layer  characteristics  for  the 
particular  tests.  For  more  information  on  various  trip  arrangements  and  relevant  data, 
please  refer  to  Appendix  A:  Trip  Arrangements. 

Following  the  selection  and  insertion  of  the  trip  arrangement,  a  zero  pressure  gradient 
flow  was  obtained.  The  wind  tunnel  ceiling  height  was  adjusted  the  entire  length  of  the 
test  section  to  ensure  a  constant  velocity  flow  [dUe/dx  =  0 .'.  dp/dx  =  0].  Velocity 

measurements  were  made  at  each  of  the  port  #’s,  see  section  2.1.6,  along  the  length  of  the 
test  section.  The  ceiling  height  was  adjusted  until  the  velocity  differences  between  each 
port  was  ±0.01  inches  of  water.  The  height  of  the  ceiling  was  adjusted  to  obtain  a  zero 
pressure  gradient  flow  for  all  single  element  tests  as  well  as  the  fetch  of  Gaussian 
roughness.  Ceiling  height  measurements  as  a  function  of  streamwise,  x,  distance  from 
the  contraction  exit  can  be  found  in  Table  2.2  below. 
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Table  2.2:  Ceiling  height  as  a  function  of  distance  from  contraction  exit  (for 

dp/dx=0  flow) 


Sin; 

gle  Elements 

Gaussian  Fetch 

Ceiling 
Height  (cm) 

nmm 

Ceiling 
Height  (cm) 

10.87 

49.53 

10.72 

49.53 

10.91 

58.42 

10.87 

58.42 

10.95 

73.66 

11.11 

73.66 

11.03 

86.36 

11.31 

86.36 

11.15 

96.52 

11.35 

96.52 

11.27 

109.22 

11.39 

109.22 

11.35 

124.46 

11.47 

124.46 

11.43 

137.16 

11.55 

137.16 

11.55 

165.1 

11.67 

165.1 

11.67 

182.88 

11.91 

182.88 

2.4.3  Profile  Comparison  and  Boundary  Layer  Characteristics 

LDV  boundary  layer  profiles  were  taken  at  various  locations  within  the  flow  field  to 
determine  two-dimensionality  and  agreement  with  direct  numerical  simulation  (DNS) 
data.  Figure  2.1 1  shows  profiles  located  at  the  center  of  the  glass  disc  (x  =  1 16.72  cm,  z 
=  0  cm),  2.26  cm  on  either  side  of  center,  and  5.08  cm  upstream  of  center.  The  different 
streamwise  locations  were  taken  because  all  axisymmetric  elements  were  placed  on  the 
glass  disc  at  the  upstream  location  whereas  the  cube  element  was  placed  at  the  center  of 
the  glass  disc.  This  different  element  location  was  to  ensure  that  all  of  the  desired  LDV 
profiles  could  be  taken  for  the  elements  within  the  glass  disc.  The  spanwise,  z-direction, 
locations  for  the  boundary  layer  profiles  were  chosen  to  be  outside  of  the  bounds  that 
LDV  measurements  would  be  taken  on  the  roughness  elements.  Meaning  that  all  LDV 
profiles  associated  with  the  roughness  elements  are  taken  within  this  spanwise  range  of 
±2.26  cm.  If  two-dimensionality  was  satisfied  at  these  locations  as  is  shown  in  Figure 
2.1 1,  then  it  is  assumed  to  be  satisfied  within  the  bounds  as  well.  The  two  sets  of  DNS 
data  are  for  a  turbulent  boundary  layer  over  a  flat  plate  and  are  from  Spalart  (1988). 
Differences  between  all  LDV  profiles  are  small  (less  than  5%)  and  are  within  the 
experimental  uncertainties,  see  Appendix  C.  Therefore,  the  flow  conditions  are  very  near 
two-dimensional. 
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The  boundary  layer  characteristics  for  the  single  elements  tested  are  as  follows:  Ue=27.5 
m/s,  Ree~7500,  5=39.3  mm,  H=5*/0=1.33,  and  UT=0.95  m/s.  Here  UT  is  determined 
from  Equation  2.1,  0  from  Equation  A.l,  and  5*  from  Equation  A.2.  The  boundary  layer 
thickness,  5,  is  defined  as  the  distance  from  the  wall  where  U=0.99Ue.  The  boundary 
layer  characteristics  for  the  fetch  of  Gaussian  roughness  are  Ue=27.5  m/s,  Ree~10500, 
5=52.2  mm,  and  UT=1.43  m/s.  All  parameters  are  determined  in  a  similar  fashion  as 
above  except  for  the  value  of  UT.  In  the  case  of  the  Gaussian  fetch,  UT  is  defined  via 
Equation  D.l.  Due  to  the  fact  that  there  are  not  any  other  DNS  data  sets  available  that 
have  a  higher  and  more  representative  Ree  value  with  which  to  compare  to  the  current 
two-dimensional  boundary  layer  data,  the  two  DNS  data  sets  from  Spalart  (1988)  are 
used.  Thus,  the  difference  between  the  DNS  profiles  and  the  experimental  profiles  is 
attributed  to  the  large  difference  in  the  Ree  values.  All  profiles,  of  the  mean  streamwise 

velocity  (U),  Reynolds  normal  stresses  (u^,v2,w2  j,  and  Reynolds  shear  stress  (uv),  do 
show  the  same  trends  as  the  DNS  data.  It  is  known  from  AGARD  (1996)  that  differences 
in  the  Ree  have  an  effect  on  the  peak  values  of  the  Reynolds  stresses.  Thus  it  is  not 
expected  that  the  experimental  profiles  match  exactly  with  the  DNS  profiles  due  to  the 
large  difference  in  the  value  of  Ree.  The  very  small  differences  between  the  LDV 
profiles  themselves,  along  with  the  same  trends  as  DNS  data,  give  confidence  in  the 
quality  of  the  flow  within  the  wind  tunnel. 


2.5  Description  of  Roughness  Elements  Tested 
2.5.1  Single  Elements 

All  single  elements  tested  were  machined  from  aluminum  and  made  in  the  Virginia  Tech 
Aerospace  and  Ocean  Engineering  machine  shop.  All  elements  were  machined  so  as  to 
have  a  post  slightly  less  than  1/16  inches  (0.159  cm)  in  diameter  protruding  from  the 
bottom  of  the  element.  This  post  was  1/8  inches  (0.318  cm)  long  in  order  to  secure  it  to 
the  glass  disc.  The  post  at  this  diameter  is  necessary  to  leave  adequate  room  for  the  glue 
used  to  secure  it  to  the  glass  disc  via  a  1/16  inch  (0.159  cm)  diameter  hole.  Seven  single 
elements  were  tested  in  total  and  are  listed  here:  cone,  Gaussian  bump,  cone  with  spatial 
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variations  equal  to  the  smallest  sublayer  structure  length  scale  {\2v/U  T )  (hereafter  fine 

grooved),  cone  with  spatial  variations  equal  to  2.5  times  the  smallest  sublayer  structure 
length  scale  (hereafter  large  grooved),  hemisphere,  cube  aligned  90°  relative  to  the  flow, 
and  a  cube  aligned  45°  relative  to  the  flow.  The  spatial  variations  of  (\2v/Ur)  were 

chosen  via  an  average  value  for  the  smallest  turbulent  length  scale.  This  value  has  been 
shown  by  Jimenez  and  Moin  (1991)  to  give  sufficient  resolution  in  order  to  resolve  the 
turbulent  scales.  Jimenez  and  Moin  showed  a  grid  resolution  in  the  spanwise  and 
streamwise  directions  equal  to  Az+=5-10  and  Ax+=8-16.  Figure  2.12  shows  three- 
dimensional  drawings  of  each  element  and  Figure  2.13  shows  close-up  views  of  the 
actual  elements.  Dimensions  of  these  elements  and  other  details  are  discussed  below. 

The  cone  has  a  height  of  0.1  inches  and  a  base  diameter  equal  to  0.1  inches.  The 
Gaussian  bump  has  a  height  of  0.1  inches,  a  base  diameter  equal  to  0.1  inches  and  has  its 
profile  defined  by  the  following  equation, 

y  =  0. 1  exp(- 1 843z 2 )  (2.2) 

where  y  is  in  the  vertical  direction  and  z  is  in  the  spanwise  direction.  The  fine  grooved 
element  has  the  same  frontal  projected  area  as  the  cone,  a  height  equal  to  0.1  inches,  and 
a  base  diameter  equal  to  0.106  inches.  The  large  grooved  element  also  has  the  same 
frontal  projected  area  as  the  cone,  a  height  of  0.1  inches,  and  a  base  diameter  of  0.1 18 
inches.  The  hemisphere  is  defined  by  a  base  diameter  equal  to  0.1 10  inches.  Finally,  the 
cube  has  an  edge  length  of  0.068  inches.  Element  profile  plots,  showing  the  variation  of 
element  width  with  respect  to  element  height,  of  the  more  complex  elements  (Gaussian, 
fine  grooved,  and  large  grooved)  can  be  seen  in  Figure  2.14.  Similarly,  Table  F.l  in 
Appendix  F  gives  the  y-z  coordinates  that  define  the  shape  of  both  the  large  and  fine 
grooved  elements. 

2.5.2  Fetch  of  Gaussian  Elements 

The  fetch  of  Gaussian  elements,  shown  in  Figures  2.15  through  2.17,  was  constructed  by 
means  of  a  vacuum  bagging  technique.  A  specially  made  Gaussian  shaped  end-mill 
cutter  with  the  same  profile  as  described  by  Equation  2.2  was  used  to  drill  holes  at  a 
constant  spacing  of  0.216  inches  into  a  12  inch  x  12  inch  x  1/8  inch  Teflon  mold.  The 
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spacing  was  determined  to  give  a  desired  value  of  A.=0.088  where  [(A.=projected  frontal 
area  to  flow)/( total  surface  area)].  Vent  holes  with  a  diameter  equal  to  0.005  inches  were 
also  drilled  at  the  center  of  each  Gaussian  shaped  cavity.  These  holes  continue  the  rest  of 
the  way  through  the  mold.  The  vent  hole  allows  the  escape  of  excess  resin  as  well  as  air 
bubbles  that  can  create  voids. 

A  layer  of  vacuum  bagging  film  was  secured  to  a  table  top  with  duct  tape.  The  mold  was 
placed  on  this  film  and  raised  off  of  the  vacuum  bagging  surface  by  metal  shims  that 
facilitate  the  escape  of  resin  and  air.  A  fiber  compound  was  then  placed  over  the  mold 
(see  note  below).  An  epoxy  mixture  was  used  as  the  base  resin  for  the  elements.  This 
mixture  was  poured  over  the  mold  and  fiber  compound.  A  layer  of  plastic  film  was  then 
placed  over  the  epoxy  in  order  to  create  a  protective  interface  between  the  vacuum 
bagging  film  and  resin.  Finally,  the  outermost  layer  of  vacuum  bagging  film  was  secured 
to  the  table.  All  leaks  and  connections  within  the  vacuum  bagging  setup  were  sealed  with 
vacuum  bagging  sealant  tape.  The  vacuum  pump  was  left  running  until  a  sample  of  the 
epoxy  resin  was  hard,  a  minimum  of  5  hours.  Once  the  resin  hardened,  the  Gaussian 
fetch  was  removed  from  the  mold.  The  ‘spikes’  created  by  the  vent  holes  could  easily  be 
removed  by  brushing  a  ruler  over  the  top  of  the  fetch.  The  fetch  was  then  cut  into  8  inch 
x  12  inch  sections  and  placed  in  the  wind  tunnel  via  rubber  cement.  The  fetch  for  the 
measurement  region  was  cut  into  an  8  inch  diameter  circle.  This  circle  was  then  aligned 
with  the  upstream  and  downstream  fetches  and  secured  to  the  glass  disc  using  a  thin  layer 
of  rubber  cement.  The  fetch  of  roughness  covered  the  tunnel  floor  from  the  end  of  the 
trip  arrangement  to  20.32  cm  past  the  downstream  edge  of  the  glass  disc. 

The  carbon- fiber  substrate  was  used  for  the  fetch  everywhere  on  the  tunnel  floor  except 
over  the  measurement  volume.  The  measurement  volume  region  utilized  a  thinner 
cheesecloth  substrate.  Various  sections  of  this  cheesecloth  substrate  needed  to  be 
removed  in  order  for  data  acquisition  to  occur.  A  thin  substrate  was  required  in  order  for 
the  ‘lip’  between  the  glass  and  the  fetch  to  be  minimized. 
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2.6  Alignment  of  LDV  Head  with  Wind  Tunnel  and 
Roughness  Elements 

Alignment  of  the  LDV  head  with  the  roughness  elements  tested  is  absolutely  essential  in 
obtaining  meaningful  and  repeatable  results.  Consequently,  very  vigilant  care  was  taken 
with  this  aspect  of  the  measurement  process.  First  the  LDV  head  and  traverse  system 
base  had  to  be  aligned  with  the  tunnel  axis.  Following  this  procedure  it  was  then  possible 
to  align  the  LDV  head  with  the  roughness  elements.  The  following  two  procedures  are 
discussed  below. 

2.6.1  Alignment  of  LDV  Head  with  Wind  Tunnel 

The  traverse  system  was  designed  and  constructed  to  have  three-orthogonal  axes.  Thus, 
it  was  necessary  to  align  only  one  of  these  axes  with  the  desired  axis  in  the  wind  tunnel. 
The  resolution  of  the  traverse  system  in  the  horizontal  (x-z  plane)  is  5  pm  whereas  the 
resolution  in  the  y-direction  is  ~ 2  pm.  These  fine  resolutions  give  significant  confidence 
in  the  final  alignment  and  position  of  the  traverse  system.  Refer  to  Figure  2.1  for  a 
description  of  the  wind  tunnel  axis  system.  The  streamwise  x-axis  was  the  axis  chosen 
for  the  alignment  of  the  traverse  system  and  the  wind  tunnel.  A  12  inch  long  1/8  inch 
thick  metal  engineer’s  ruler  was  placed  on  the  centerline  of  the  wind  tunnel  over  the  glass 
disc  measurement  area.  The  measurement  volume  of  the  LDV  head  was  then  placed  on 
the  surface  of  the  glass  disc  at  the  edge  of  the  ruler.  The  LDV  head  was  traversed  in  the 
x-direction  along  the  edge  of  the  ruler  over  the  entire  measurement  range.  If  the 
measurement  volume  strayed  in  any  one  direction  over  this  range,  the  traverse  base  was 
moved  and  the  procedure  repeated.  After  a  number  of  iterations  it  was  possible  to  get  the 
measurement  volume  to  follow  the  centerline  of  the  tunnel  over  the  entire  measurement 
range.  A  secondary  check  was  made  along  the  spanwise  z-axis  of  the  wind  tunnel  to 
ensure  the  orthogonality  of  the  traverse  system. 

2.6.2  Alignment  of  LDV  Head  with  Roughness  Elements 
2.6.2.1  Single  Elements 

All  single  elements  tested  were  placed  on  a  20.32  cm  diameter  0.64  cm  thick  piece  of 
float  glass  through  which  LDV  measurements  could  be  taken.  Each  element  was 
machined  so  as  to  have  a  post  slightly  less  than  1/16  inches  (0.159  cm)  in  diameter 
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protruding  from  the  bottom  of  the  element.  This  post  was  1/8  inches  (0.318  cm)  long  in 
order  to  secure  it  to  the  glass.  The  float  glass  disc  had  a  1/16  inch  (0.159  cm)  diameter 
hole  drilled  1/8  inches  (0.318  cm)  deep  into  it  with  a  diamond  coated  glass  drill  bit  from 
Precision  Diamond  Drills.  The  roughness  element  was  then  secured  to  the  glass  disc  with 
super  glue.  All  axisymmetric  elements  were  placed  2  inches  (5.08  cm)  upstream  of 
center  on  the  glass  disc.  This  extra  downstream  space  was  necessary  in  providing  the 
LDV  head  access  to  all  measurement  locations.  The  cube  element  was  placed  at  the 
center  of  the  disc  in  order  to  facilitate  the  taking  of  measurements  at  the  90°  and  45° 
orientations  relative  to  the  oncoming  flow.  As  is  discussed  previously,  boundary  layer 
profiles  were  taken  at  both  locations  (center  and  5.08  cm  upstream)  on  the  glass  disc  and 
there  was  no  discemable  difference  between  the  boundary  layer  characteristics. 

The  alignment  procedure  for  the  cone,  Gaussian,  fine  grooved,  and  large  grooved  element 
was  identical.  The  measurement  volume  was  placed  at  the  surface  of  the  glass  and  the 
LDV  head  was  aligned  along  the  centerline  of  the  tunnel.  With  the  LDV  head  aligned 
with  the  centerline  of  the  tunnel,  the  probe  was  traversed  in  the  vertical  direction  while 
moving  the  glass  disc  so  the  beams  would  be  at  the  edges  and  center  of  the  element.  The 
LDV  head  was  then  traversed  in  the  vertical  direction  until  the  measurement  volume  was 
at  the  height  of  the  roughness  element.  Using  laser  goggles  it  was  possible  to  follow  the 
beams  up  the  side  of  the  element  until  it  reached  the  peak.  The  glass  disc  was  rotated  in 
the  tunnel  in  order  to  do  the  final  alignment  of  the  peak  of  the  element  with  the 
measurement  volume.  By  watching  the  beam  profiles  change  along  the  sides  of  the 
element  and  the  wall  of  the  tunnel  (when  the  measurement  volume  was  higher  than  the 
element  height)  it  was  possible  to  have  an  uncertainty  in  positioning  of  ±60  microns  in 
the  x  and  z  directions.  After  the  element  was  aligned  with  the  probe  volume,  the  glass 
disc  was  taped  to  the  wind  tunnel  floor  ensuring  that  no  positioning  change  occurred. 

The  hemisphere  was  aligned  in  a  similar  fashion  as  that  described  above.  The  only 
difference  being  that  the  probe  volume  could  not  be  aligned  with  the  element  peak.  The 
head  was  traversed  back  and  forth  along  the  base  of  the  hemisphere  and  the  glass  disc 
rotated  until  the  measurement  volume  was  moved  the  same  distance  from  the  center  of 
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the  tunnel  to  each  spanwise  edge  of  the  element.  A  secondary  check  that  ensured 
alignment  was  traversing  the  LDV  head  in  the  streamwise  direction  and  checking  that  the 
distance  from  the  spanwise  located  element  center  was  also  equidistant  from  the 
streamwise  edges  of  the  hemisphere. 

The  cube  alignment  method  is  as  follows.  With  the  glass  disc  outside  of  the  tunnel,  a  fine 
scaled  ruler  (1/100  inch  divisions)  was  used  to  find  the  center  of  one  edge  of  the  cube.  A 
jeweler’s  eye  loop  was  utilized  for  this  procedure.  Once  the  center  of  the  cube  was 
found,  the  ruler  was  taped  to  the  glass  disc.  A  line  on  the  edge  of  the  glass  disc  was  then 
made  with  a  permanent  marker  that  signified  the  center  of  the  cube’s  edge  (hereafter  the 
90°  line).  Using  this  90°  line  another  line  was  drawn  at  an  angle  of  45°.  The  centerline 
of  the  tunnel  was  then  marked  and  the  glass  disc  inserted  into  the  tunnel  floor.  After 
aligning  the  90°  line  on  the  disc  with  the  tunnel  centerline  a  series  of  secondary  checks 
was  performed.  The  measurement  volume  was  traversed  along  the  ruler’s  edge  (note  that 
the  cube’s  edge  is  much  too  short  to  ensure  proper  alignment)  ensuring  that  no  deviation 
of  the  measurement  volume  could  be  detected  from  the  edge  of  the  ruler.  Finally,  the 
distance  from  the  now  known  element  center  to  each  edge  of  the  cube  was  checked  to 
guarantee  that  the  center  of  the  element  was  known.  If  anything  did  not  check  out  with 
the  previous  measurements  and  alignments,  the  entire  process  was  repeated. 

Alignment  for  the  cube  rotated  45°  relative  to  the  flow  consisted  of  rotating  the  glass  disc 
until  the  45°  line  on  the  disc  was  aligned  with  the  centerline  of  the  wind  tunnel. 
Secondary  checks  were  again  made  to  ensure  that  all  comers  were  aligned  with  the  tunnel 
coordinate  system  and  equidistant  from  the  center  of  the  element. 

Once  the  center  of  each  roughness  element  was  located,  the  traverse  system  zero  location 
was  determined.  All  measurement  locations  were  relative  to  this  center  location. 

2.6.2.2  Fetch  of  Gaussian  Elements 

The  majority  of  the  Gaussian  roughness  was  made  on  a  carbon-fiber  substrate  in  sections 
that  were  8  inches  (20.32  cm)  wide  by  12  inches  long.  These  sections  were  centered  in 
the  wind  tunnel  and  rubber  cemented  to  the  floor.  Measurements  from  each  side  of  the 
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roughness  fetch  to  its  respective  tunnel  wall  were  made  to  ensure  that  the  fetch  was 
centered  in  the  test  section.  The  roughness  extended  from  the  end  of  the  trip  arrangement 
to  8  inches  (20.32  cm)  past  the  end  of  the  glass  disc.  Figure  2.16  shows  a  view  with  the 
wind  tunnel  ceiling  removed  looking  upstream  at  the  roughness  fetch. 

As  is  discussed  previously,  the  fetch  of  Gaussian  roughness  that  was  made  to  cover  the 
measurement  region  has  a  cheesecloth  substrate.  In  order  to  obtain  measurements  from 
this  setup,  sections  of  cheesecloth  had  to  be  removed  for  the  passage  of  the  laser  beams, 
see  Figure  2.17.  Please  note  that  all  sections  that  were  removed  upstream  were  replaced 
before  taking  measurements  downstream.  Very  few  Gaussian  elements  were  missing 
from  the  fetch  of  roughness.  The  ones  that  were  missing  were  far  enough  away  from  the 
measurement  region  in  the  upstream,  downstream,  and  spanwise  directions  so  as  to  not 
affect  the  obtained  measurements.  This  conclusion  was  determined  from  previous  results 
and  experience  done  by  George  and  Simpson  (2002  and  2004).  A  He-Ne  laser  located  on 
the  optical  table  was  used  to  align  the  probe  with  the  desired  Gaussian  element.  The  red 
beam  coming  from  the  laser  was  launched  into  the  receiving  optical  fiber  and  thus 
directed  out  of  the  receiving  lens  and  into  the  wind  tunnel.  It  was  then  a  matter  of 
iteration  to  center  the  LDV  head  on  the  Gaussian  element.  The  probe  was  traversed 
vertically  until  the  red  beam  could  be  seen  through  the  element.  The  probe  was  then 
traversed  in  the  streamwise  and  spanwise  directions  until  the  center  of  the  beam 
extending  from  the  receiving  lens  was  at  the  center  of  the  peak  of  the  Gaussian  element. 
The  traverse  system  zero  location  was  defined  and  measurement  locations  were  relative 
to  the  center  of  the  element. 


Chapter  2:  Apparatus  and  Instrumentation 


-34- 


2.7  Chapter  2  Figures 


Figure  2.1 :  Schematic  of  small  boundary  layer  wind  tunnel  setup  (Adapted  from 

Smith,  et  al.  1990) 


1 .  Air  Conditioner 

2.  Filter  Box 

3.  Blower 

4.  Drip  Pan 

5.  Baffle  Plate 

6.  Seeding  Pipe 

7.  Plenum  Chamber 

8.  Contraction 

9.  Honeycomb 


10.  Screens 

11.  Contraction 

12.  Trip  Arrangement 

13.  Glass  Side  Walls 

14.  Bottom  Wall 

15.  Top  Wall 

16.  Return  Ducting 

17.  Rubber  Joint 
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Figure  2.2:  Small  boundary  layer  wind  tunnel  setup  (view  from  front) 


Figure  2.3:  Small  boundary  layer  wind  tunnel  setup  (view  from  back) 
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Figure  2.4:  View  of  baffle  plate  in  plenum  chamber  (looking  from  blower) 


Figure  2.5:  Seeding  tube  and  drip  pan  on  downstream  side  of  baffle  plate 
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Figure  2.6:  Optical  table  set-up  shown  near  laser  exit 


Figure  2.7:  Full  view  of  optical  table 
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Figure  2.8:  Close-up  view  of  LDV  head 


Figure  2.9:  LDV  head  during  data  acquisition  (view  is  above  the  tunnel  ceiling 

looking  upstream) 
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Figure  2.10:  Spanwise  velocity  profile  at  center  of  wind  tunnel 
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Figure  2.12:  Three-Dimensional  drawings  of  single  element  cases 
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Figure  2.14:  Element  profile  plots  showing  element  width  varying  with  respect  to  element  height 
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Figure  2.15:  Close-up  view  of  Gaussian  roughness  fetch  (center-to-center 

spacing=0.216  inches) 


Figure  2.16:  Top  view  of  test  section,  looking  upstream,  with  wind  tunnel  ceiling 

removed 
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Chapter  3  Results 


This  chapter  discusses  the  results  of  the  data  obtained  within  the  analyzed  flow  fields  of 
the  single  elements  and  Gaussian  fetch  of  roughness.  Comparisons  between  elements  are 
made  for  all  mean  velocities,  turbulence  quantities,  and  skin  friction.  Roughness  element 
comparisons  can  also  be  found  for  other  fluid  dynamic  quantities  derived  from  the 
previous  parameters.  Oil  flow  visualizations  are  also  presented  and  described  within  this 
chapter.  Please  note  that  the  spanwise  normalization  scheme  uses  the  parameter  (d). 
This  value  for  all  axisymmetric  elements  is  the  base  diameter  of  the  element.  For  brevity, 
this  same  parameter  is  the  edge  length  of  the  cube  (used  to  normalize  the  cube  at  90°)  and 
the  diagonal  of  the  cube  (used  to  normalize  the  cube  at  45°). 


3.1  Measurement  Locations 

The  twenty-six  boundary  layer  profile  measurement  locations  for  the  single  elements  can 
be  seen  in  Figures  3.1  and  3.2.  In  an  effort  to  keep  the  data  and  figures  of  the  isolated 
roughness  elements  separate  from  the  Gaussian  fetch  of  roughness,  the  figures  associated 
with  the  Gaussian  fetch  can  be  found  following  the  figures  related  to  the  isolated 
roughness  elements.  As  a  result,  the  thirteen  measurement  locations  for  the  fetch  of 
Gaussian  shaped  roughness  can  be  seen  in  Figure  3.113.  Single  point  LDV 
measurements,  18  to  24  in  number,  were  taken  in  logarithmically  spaced  increments  at 
each  shown  profile  location.  All  measurements  were  taken  in  the  downstream  flow  field 
of  the  elements.  The  coordinate  systems  used  are  also  defined  in  the  figures  listed  above. 
The  origin  of  the  coordinate  systems  is  located  at  the  center  of  each  roughness  element 
under  investigation.  The  x-direction  is  taken  to  be  along  the  centerline  of  the  wind  tunnel 
with  (+)x  being  downstream  of  the  element.  The  y-direction  is  normal  to  the  wind  tunnel 
floor  (wall)  where  (+)y  is  up  from  the  floor.  The  coordinate  system  is  completed 
following  the  right  hand  rule.  Please  note  that  all  axisymmetric  single  elements  were 
placed  at  a  distance  of  1 1 1.64  cm  from  the  contraction  exit  and  the  cube  elements  were 
placed  1 16.72  cm  from  the  contraction  exit.  There  was  no  discemable  difference  in  the 
boundary  layer  characteristics  at  these  two  locations  (see  section  2.4.3). 
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3.1.1  Single  Element 

Centerline  measurements  were  taken  at  six  downstream  x/d  locations  in  order  to 
determine  the  distance  the  boundary  layer  needed  to  recover  from  the  presence  of  the 
roughness  element.  These  centerline  locations  are  x/d=1.36,  2.75,  5.0,  10.0,  20.0,  and 
40.0.  Please  note  that  the  last  downstream  measurement  location  for  the  cube  at  45°  is 
located  at  x/d=31.0.  This  is  due  to  insufficient  room  to  traverse  the  LDV  probe 
downstream  and  still  take  measurements  through  the  glass  disc.  Two  different  planes  of 
data  were  taken  at  x/d  locations  of  2.75  and  10.0  for  all  elements  except  the  fine  grooved 
element.  Data  for  the  fine  grooved  element  were  taken  at  an  x/d=1.36  plane.  The 
x/d=1.36  plane  was  initially  thought  to  provide  further  insight  into  the  vortical  structures 
emanating  from  the  peak  of  the  roughness  element.  All  planes  of  data  make  it  possible  to 
determine  contours  of  both  measure  and  derived  quantities.  Hundreds  of  data  points  were 
used  to  formulate  each  grid  used  to  make  the  contour  plots.  Measurements  were  mainly 
taken  in  the  -z  direction  for  symmetry  reasons.  One  profile  was  taken  in  the  +z  direction 
at  the  centerline  locations  in  order  to  assure  that  the  flow  was  symmetric.  Results  for  the 
single  elements  are  shown  with  the  reference  smooth  wall  results.  The  smooth  wall 
profile  was  taken  at  the  location  of  the  element  without  the  element  being  present. 

3.1.2  Gaussian  Fetch 

Boundary  layer  profiles  for  the  Gaussian  fetch  of  roughness  are  similar  to  those  for  the 
single  element  cases  described  previously.  Three  centerline  locations,  x/d=0.79,  1.08, 
and  1.37,  were  taken  downstream  of  the  chosen  roughness  element.  A  plane  of  data  at 
the  x/d=l  .08  location  (halfway  between  elements)  was  again  taken  in  order  to  develop 
contours  of  the  measured  and  derived  quantities.  Two  profiles  were  taken  in  the  +z 
direction  at  this  plane  location  in  order  to  assure  flow  symmetry.  A  profile  directly 
between  roughness  elements  was  also  taken  at  x/d=2.16  and  z/d=-1.08.  Finally  a 
centerline  profile  at  x/d=1.08  was  taken  at  a  location  8  cells  downstream  of  the  shown 
configuration.  For  data  and  a  brief  analysis  related  to  the  Gaussian  fetch  of  roughness, 
see  Appendix  D.  A  more  thorough  analysis  of  this  data  set  will  be  done  by  Stewart 
(2004). 
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3.2  Surface  Oil  Flow  Visualizations 

The  oil  flow  visualization  technique  helps  to  give  a  qualitative  idea  of  what  the  global 
flow  field  looks  like  on  the  wall.  The  surface  oil  flow  visualizations  were  done  with  a 
mixture  of  1 5  ml  titanium  dioxide,  40  ml  kerosene,  and  1  ml  oleic  acid.  This  mixture  was 
brushed  onto  a  piece  of  black  contact  paper  that  was  attached  to  each  element’s 
respective  glass  disc.  The  tunnel  was  then  turned  on  and  left  running  until  the  mixture 
was  dry.  Digital  pictures  of  the  oil  flows  done  for  each  element  can  be  seen  in  Figures 

3.3  and  3.4.  The  vertical  lines  shown  in  Figure  3.3  and  Figure  3.4  are  the  streamwise 
locations  where  planes  of  data  were  obtained. 

In  front  of  all  elements  it  is  possible  to  see  the  flow  separation  and  related  formation  of  a 
horseshoe  vortex.  For  each  cube  element  there  is  evidence  of  high  scouring,  denoted  by 
the  black  regions,  directly  behind  the  element.  Consequently,  this  suggests  the  formation 
of  an  arch  vortex  which  promotes  backflow  in  the  downstream  regions  near  the  element. 
This  vortex  looks  to  be  similar  in  both  cases  of  when  the  cube  is  rotated  45°  relative  to 
the  flow  and  when  the  cube  is  oriented  at  90°  relative  to  the  flow.  It  can  also  be  seen  that 
there  is  a  significant  region  directly  behind  the  cube  elements  where  the  wake  region 
experiences  a  ‘necking’  down  effect.  In  the  close-up  views  of  the  elements,  Figures  3.5 
and  3.6,  flow  separation  is  easily  noticed  at  the  sharp  edges  of  the  cubes.  There  is  a 
larger  region  of  flow  separation  in  the  cube  at  45°  as  would  be  expected  at  this  particular 
orientation.  Further  downstream  of  the  cube  elements,  lines  of  particle  deposits  can  be 
seen.  These  deposits  give  rise  to  the  conclusion  that  there  is  a  significantly  lower  amount 
of  wall  shear  in  these  areas  as  compared  to  the  areas  along  the  centerline  of  the  element. 
Flow  reattaches  behind  the  arch  vortex  region  as  the  accelerating  fluid  over  the  top  and 
around  the  sides  of  each  element  converges  toward  one  another  near  the  centerline.  A 
significant  amount  of  scouring  is  also  seen  directly  behind  the  hemisphere  element.  Fluid 
rushing  over  the  top  of  the  element  impinges  on  the  wall  and  causes  the  noticeable  black 
region.  Flow  separation  can  also  be  seen  on  the  back  half  of  the  hemisphere,  indicated  by 
the  line  of  tracer  particles  extending  downstream  from  the  edges  of  the  element.  Again 
shear  layer  roll-up  is  found  in  front  of  the  element  which  leads  to  the  formation  of  a 
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horseshoe  vortex  structure.  The  wake  behind  the  hemisphere  looks  to  be  spreading  in  the 
spanwise  direction  fairly  rapidly. 

All  peaked  elements  in  Figure  3.4  seem  to  have  a  similar  wake  structure.  The  formation 
of  a  horseshoe  vortex  structure  can  be  seen  in  front  of  the  elements  which  in  turn  extends 
downstream  of  the  elements.  All  elements  have  a  region  of  flow  separation  and  backflow 
directly  behind  the  element  followed  by  a  region  of  particle  deposits.  Directly 
downstream  of  these  deposits  flow  reattaches  due  to  the  fluid  coming  around  the  sides 
and  over  the  peaks  of  the  element  to  converge  towards  the  floor.  The  wake  behind  the 
‘peaked’  elements  tends  to  be  more  of  a  wedge  shape  as  compared  to  the  cubes  and 
hemisphere.  Regions  of  tracer  particles  are  also  noticeable  away  from  the  centerline  of 
the  elements  which  again  show  a  lower  region  of  wall  shear. 

As  discussed  previously,  away  from  the  centerline  in  the  spanwise  z-direction  there  are 
noticeable  deposits  of  tracer  particles  which  indicate  a  region  of  lower  shear  as  compared 
to  the  centerline.  The  regions  of  lower  shear  are  attributed  to  less  mixing  due  to  the 
counter-rotating  vortices  creating  a  region  of  upwash.  Whereas,  in  the  high  regions  of 
shear  there  will  be  greater  mixing  due  to  the  downwash  of  the  counter-rotating  vortex 
pair.  There  will  also  be  some  fluid  acceleration  along  the  centerline  that  will  have  a 
slight  affect  on  the  increase  of  wall  shear.  The  presence  of  the  wake  can  be  seen  in  all 
elements  to  around  ten  diameters  downstream.  The  axisymmetric  elements’  wakes  tend 
to  dissipate  first  followed  by  the  cube  at  90°  and  then  the  cube  at  45°.  This  phenomenon 
suggests  that  the  strength  of  the  vortex  structure  produced  by  the  elements’  presence  in 
the  turbulent  boundary  layer  is  significantly  stronger  in  the  cube  at  45°. 


3.3  Mean  Velocities 

Centerline  profiles  of  the  mean  streamwise  velocity  are  presented  in  non-dimensional 
form  as  U/Ur  vs.  yUjv  in  Figure  3.7.  The  value  used  for  the  skin  friction  velocity,  UT, 


is  from  the  approach  boundary  layer  just  upstream  of  a  given  element.  Contours  of  the 
same  mean  velocity  for  x/d  locations  of  2.75  (1.36  for  the  fine  grooved  cone)  and  10.0 
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can  be  seen  in  Figures  3.8  through  3.11.  At  the  first  measurement  location  of  x/d=1.36, 
all  elements  have  a  backflow  region  present  except  for  the  Gaussian  element.  The  largest 
magnitude  of  backflow  can  be  seen  in  the  cube  elements.  Here  the  cube  at  45°  has  a  U+ 
value  very  near  -7.  The  fine  grooved  element  and  hemisphere  are  very  comparable 
initially,  but  the  hemisphere  recovers  to  the  reference  smooth  wall  sooner  due  to  its  lower 
height.  The  cone  and  large  grooved  elements  follow  with  still  smaller  values  of 
backflow.  It  is  reasonable  to  surmise  that  the  large  grooved  element  causes  less  blockage 
to  the  oncoming  flow,  due  to  its  sizeable  grooves,  and  in  return  has  a  smaller  backflow 
region.  All  elements  have  a  high  velocity  gradient  near  their  respective  heights  as  all 
profiles  show  a  collapse  to  the  reference  2D  boundary  layer  above  the  element  heights. 
Each  element  collapses  in  the  order  of  its  height  except  the  cube  at  45°.  This  element 
tends  to  a  later  collapse  than  the  cube  at  90°.  A  drastic  change  is  seen  in  the  profiles  by 
x/d=2.75.  At  this  location  the  separated  flow  downstream  of  all  elements  has  reattached 
and  the  flow  for  most  of  the  axisymmetric  elements  has  collapsed  to  one  another.  The 
cube  elements  still  have  the  most  significant  velocity  deficit  due  to  the  amount  of 
blockage  and  flow  separation  that  they  create.  At  x/d=5.0  the  mean  U /U T  profiles  for  all 

elements  have  collapsed  to  one  another  except  for  the  cube  at  90°.  Profiles  for  all 
elements  collapse  to  the  2D  reference  boundary  layer  slightly  above  their  height.  The 
cube  at  90°  is  the  only  element  that  shows  signs  of  a  velocity  defect  at  x/d=10.0,  and  by 
x/d=20.0  the  profiles  for  all  elements  have  collapsed  to  the  reference  boundary  layer 
values.  At  x/d=40.0  there  are  no  effects  of  the  elements’  presence  in  the  boundary  layer. 
The  U+  velocity  contours  in  Figures  3.8  through  3.11  show  the  spanwise  effect  that  the 
elements  have  on  the  oncoming  flow.  As  is  expected,  the  cubes  show  the  most 
significant  effect  at  the  upstream  location.  By  x/d=10.0  most  discemable  effects  have 
dissipated  for  all  elements.  Please  note  that  in  all  contour  plots,  crude  outlines  of  the 
tested  elements  are  shown.  Due  to  the  fact  that  the  plots  are  in  log  scale,  the  outlines  are 
not  drawn  to  scale.  These  outlines  are  simply  presented  in  order  to  give  a  general 
understanding  of  the  height  of  each  element  with  respect  to  the  presented  contours. 

Centerline  profiles  of  the  mean  normal-to-wall  velocity  are  presented  in  non-dimensional 
form  as  V/Ur  vs.  yUjv  in  Figure  3.12.  Contours  of  the  same  mean  velocity  for  x/d 
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locations  of  2.75  (1.36  for  fine  grooved)  and  10.0  can  be  seen  in  Figures  3.13  through 
3.16.  At  the  first  measurement  location,  x/d=1.36,  there  is  an  extremely  significant 
amount  of  fluid  being  brought  toward  the  wall  in  all  cases.  In  all  cases  the  peak  value  is 
reached  at  the  roughness  element  height  or  slightly  below.  Most  elements  have  a  peak 
value  at  the  forward  most  location  of  around  V+~-3.3  with  the  cube  at  90°  being 
significantly  less  at  V+~-1.9.  This  large  difference  can  be  attributed  to  the  increased 
blockage  effects  of  the  cube.  There  will  be  a  region  of  separated  and  more  stagnant  flow 
directly  behind  the  element.  The  profiles  downstream  of  the  elements  do  not  collapse  to 
the  reference  boundary  layer  until  well  above  their  heights.  No  effects  can  be  seen  above 
three  element  heights.  At  x/d=2.75  there  is  a  very  drastic  difference  in  the  behavior  of 
the  normal-to-wall  mean  velocity  between  elements.  The  cube  at  45°  has  a  peak  value  of 
~-3.4  which  is  over  twice  the  amount  seen  by  any  other  element.  This  drastic  difference 
is  caused  by  the  large  downwash  of  fluid  coming  over  the  top  of  the  element  along  with 
the  strong  horseshoe  vortex  structure  bringing  a  large  amount  of  fluid  toward  the  wall. 
The  cube  at  90°  has  the  next  largest  magnitude  followed  by  a  close  conglomeration  of  the 
peaked  elements  and  finally  the  hemisphere  follows  with  the  least  amount  of  normal-to- 
wall  velocity.  All  elements  present  a  noticeable  peak  that  occurs  below  the  height  of  the 
element. 

The  magnitudes  of  velocity  continue  to  decrease  at  x/d=5.0  with  the  cube  at  45°  still 
having  the  largest  effect  on  the  boundary  layer.  At  x/d=10.0  the  hemisphere  has 
collapsed  to  the  reference  boundary  layer  and  all  elements,  save  the  cube  at  45°,  have 
collapsed  to  each  other.  At  x/d=20.0  and  40.0  the  elements  have  collapsed  back  to  the 
reference  2D  boundary  layer.  As  is  evidenced  by  the  contour  plots  in  Figures  3.13  and 
3.14,  the  wall-ward  rush  of  fluid  is  greatest  near  the  centerline  with  a  peak  value  slightly 
below  the  element  height.  The  velocity  magnitude  decreases  with  spanwise  distance 
from  the  centerline.  The  most  extensive  effects  in  the  spanwise  direction  can  be  seen  in 
the  cubes  followed  by  the  cone,  Gaussian,  fine  grooved  (estimate  due  to  measurement 
location),  large  grooved,  and  finally  the  hemisphere.  A  noticeable  peak  of  normal-to- 
wall  velocity  can  still  be  seen  near  the  element  height,  albeit  a  small  amount,  in  all 
elements  except  for  the  hemisphere.  The  streamwise  decay  of  the  V  velocity  is  fairly 
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drastic  near  the  element.  This  decay  smoothes  out  substantially  as  the  streamwise 
distance  is  increased. 

Contours  of  the  mean  normalized  spanwise  velocity  component,  W+,  can  be  seen  in 
Figures  3.17  through  3.20.  These  plots  show  a  significant  amount  of  velocity  in  the 
spanwise  direction  for  the  cube  elements  as  compared  to  the  peaked  elements  and  the 
hemisphere.  Flow  symmetry  can  be  confirmed  by  these  plots,  thus  giving  confidence  to 
the  quality  of  flow  obtained  for  these  extremely  small  roughness  elements.  More 
discussion  concerning  the  spanwise  velocity  and  related  vorticity  can  be  found  in  section 
3.4.1. 


3.4  Derived  Fluid  Dynamic  Quantities 


3.4.1  Streamwise  Vorticity 

In  order  to  analyze  the  horseshoe  vortex  structure  more  thoroughly,  contour  plots  of 
streamwise  vorticity  (Qx)  for  all  single  elements  are  presented  in  Figures  3.21  through 
3.24.  The  streamwise  vorticity  for  all  elements  was  calculated  using  Equation  3.1 . 


Q  =  cur  IV  = 


d/dx  d/dy  djdz 
U  V  W 


n  _  dW  dv 

dy  dz 


(3.1) 


The  boundary  layer  profiles  used  to  calculate  the  derivatives  involved  with  the 
streamwise  vorticity  were  first  interpolated  to  a  fine  spacing.  This  interpolation  step  is 
important  in  order  to  have  sufficient  resolution  to  perform  numerical  differentiation  on 
the  data.  The  contours  of  streamwise  vorticity  are  in  the  y-z  plane  at  x/d  locations  of  2.75 

(1.36  for  fine  grooved)  and  10.0.  The  vorticity  is  normalized  by  JaJ jlJT  where  Af  is 
the  projected  frontal  area  of  each  element.  George  and  Simpson  (2004)  found  that  the 
length  scale  was  the  most  appropriate  parameter  to  use  to  normalize  the  vorticity. 

The  normalized  secondary  flow  vectors  V/Ut  and  WjUT  are  also  shown  with  the 
contours  of  streamwise  vorticity. 
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At  the  x/d=2.75  element  locations,  the  secondary  velocity  vectors  show  a  large  amount  of 
fluid  being  brought  toward  the  wall  by  the  counter-rotating  horseshoe  vortex  structures. 
All  elements  at  this  location  show  a  significant  normal-to-wall  (V)  velocity  at  the 
centerline.  These  nearly  perpendicular  velocity  vectors  also  show  flow  symmetry  about 
the  centerline  plane.  The  velocity  vectors  curve  sharply  toward  and  away  from  the 
centerline  near  the  height  of  the  elements  and  at  the  wall,  respectively.  The  flow  shown 
near  the  wall  is  nearly  parallel  to  the  wall  in  the  negative  z-direction.  The  secondary 
velocity  vectors  for  the  fine  grooved  element  at  the  x/d=1.36  plane  show  fluid  being 
swept  in  around  the  sides  of  the  element.  At  this  location  the  horseshoe  vortex  structure 
has  not  yet  taken  over  the  secondary  flow  field.  Therefore,  the  dominant  secondary  flow 
structure  appears  to  be  the  flow  coming  around  the  sides  of  the  element  rather  than  the 
counter-rotating  structure  shown  at  the  x/d=2.75  plane.  The  velocity  vectors  along  the 
centerline  directly  behind  the  element  are  very  small  which  also  denotes  that  the 
horseshoe  vortex  structure  has  not  yet  started  to  bring  large  amounts  of  fluid  toward  the 
wall.  The  highest  magnitude  of  secondary  flow  can  be  seen  in  the  two  cube  elements, 
followed  by  the  hemisphere,  cone,  Gaussian,  and  the  large  grooved  element.  As  is 
evidenced  by  the  V+  plots  discussed  previously  (Figure  3.12),  the  fine-grooved  element’s 
secondary  flow  at  the  x/d=2.75  plane  is  very  near  the  magnitude  of  the  Gaussian  element. 
Further  downstream  at  x/d=10.0,  the  secondary  flow  has  lost  a  significant  amount  of 
strength  as  compared  to  the  x/d=2.75  plane.  The  cube  at  45°  shows  the  largest  magnitude 
of  secondary  flow  followed  by  the  cube  at  90°.  All  other  elements’  secondary  flow 
vectors  are  very  similar  in  magnitude.  However,  the  fine-grooved  element  is  slightly 
larger.  Flow  is  still  being  directed  toward  and  away  from  the  centerline  for  all  elements 
in  a  similar  fashion  as  discussed  previously.  A  noticeable  difference  in  the  secondary 
flow  structure  is  the  increase  in  distance  away  from  the  wall  and  centerline  as  compared 
to  the  upstream  measurement  plane.  This  incident  is  due  to  the  induced  velocity  from  the 
‘wall  image  vortex’  on  the  actual  horseshoe  vortex  structure. 

Readily  identifiable  regions  of  positive  streamwise  vorticity  can  be  seen  at  the  upstream 
measurement  plane  for  all  elements.  The  cube  elements  have  the  highest  pockets  of 
vorticity  as  well  as  the  most  spanwise  skewed  contour  lines.  The  presence  of  negative 
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vorticity  very  close  to  the  wall  is  also  present  at  this  measurement  location.  This  near¬ 
wall  opposite  sign  vorticity  is  due  to  the  ‘no  slip’  condition  of  the  wall  creating  large 
velocity  gradients  in  the  spanwise,  W,  velocity.  Furthermore,  the  presence  of  these 
negative  regions  originates  from  the  convection  of  wall  induced  vorticity  caused  by  the 
interaction  of  the  horseshoe  vortex  at  the  wall.  The  negative  regions  of  vorticity  at  the 
upstream  measurement  locations  are  smaller  than  and  not  as  developed  as  the  primary 
positive  vortex  structures.  Although  not  as  developed  as  the  positive  regions,  these 
negative  regions  of  vorticity  impact  the  primary  vortex  by  generating  an  exchange  of 
vorticity  that  is  created  by  near  wall  induced  velocities.  As  the  vortices  propagate 
downstream  to  the  x/d=10.0  plane,  the  negative  regions  of  vorticity  develop  and  become 
larger.  As  is  shown  in  Figures  3.23  and  3.24,  the  size  of  the  negative  regions  of  vorticity 
are  very  close  to  the  same  size  as  their  respective  positive  regions.  At  the  downstream 
location,  the  negative  vorticity  is  of  greater  magnitude  than  that  of  the  positive  vorticity. 
These  negative  vorticity  values  are  comparable  to  the  upstream  positive  vorticity  values. 

Following  previous  studies  such  as  Wendt  and  Greber  (1992),  the  location  of  the  vortex 
center  is  taken  to  be  the  point  of  maximum  streamwise  vorticity,  Qx_max-  However,  it  is 
important  to  note  that  in  general  the  V=W=0  location  is  not  the  same  as  the  Qx_max 
location.  Locations  of  the  vortex  centers  for  each  element  can  be  found  in  Table  3.1 
below.  Also  presented  are  the  maximum  vorticity,  ju r)  ,  and  streamwise 

velocity  at  the  vortex  center  location. 
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Table  3.1:  Locations  of  maximum  streamwise  vorticity 


x/d=2.75  ( fine  grooved  x/d=1.36) 

x/d=10.0  for  all  elements 

Element 

+ 

Zc 

IHi 

a (s'1) 

+ 

Zc 

+ 

yc 

uc+ 

l  u-  L 

a (s'1) 

Cube  at  45° 

-115.4 

13.6 

9.3 

11.35 

5263 

44.0 

14.2 

2.46 

1143 

Cube  at  90° 

-81.6 

8.8 

7.0 

11.59 

6394 

32.9 

12.9 

1.63 

899 

Cone 

-120.0 

12.0 

9.3 

6.80 

3607 

-200.7 

45.7 

14.8 

0.68 

361 

Gaussian 

-79.9 

13.0 

9.1 

5.03 

2940 

-160.7 

44.5 

14.4 

0.68 

398 

-142.6 

13.9 

9.5 

4.44 

2353 

-159.1 

43.9 

14.1 

0.68 

362 

-87.8 

14.5 

10.3 

7.64 

4159 

-221.5 

57.9 

16.0 

0.73 

395 

Fine  Grooved 

-84.3 

5.6 

5.3 

23.63 

12533 

-170.4 

43.7 

14.2 

0.99 

527 

At  the  x/d=2.75  location,  the  horseshoe  vortex  structures  are  located  -9-15  wall  units  off 
of  the  wall  and  -80-143  wall  units  from  the  centerline.  In  order  to  relate  the  location  of 
maximum  streamwise  vorticity  to  the  location  of  the  sides  of  each  element,  the 
normalized  half  diameter  values  (d+/2)  for  each  element  are  listed  below.  Each  value  is 
in  the  same  order  as  its  respective  element  shown  in  Table  3.1.  The  values  of  d+/2  for 
each  element  are  69.6,  49.2,  72.4,  72.4,  85.5,  79.7,  and  79.8,  respectively.  As  can  be  seen 
in  Table  3.1,  the  hemisphere  and  Gaussian  elements  show  a  location  of  maximum 
streamwise  vorticity  that  is  located  very  close  to  the  half  diameter  location  of  each 
element.  On  the  other  hand,  all  of  the  other  elements  tend  to  push  the  point  of  Qx_max 
further  away  from  the  centerline.  For  the  fine  grooved  element  at  x/d=1.36  the  vortex 
structure  is  located  —6  units  off  of  the  wall  and  -84  units  away  from  the  centerline.  The 
vortex  structure  produced  by  the  cube  at  90°  is  centered  the  closest  to  the  wall  whereas 
the  structure  produced  by  the  hemisphere  is  the  farthest  from  the  wall.  The  maximum 
value  of  streamwise  vorticity  at  x/d=2.75  is  seen  in  the  cube  elements.  The  cube  at  90° 
shows  the  largest  magnitude  with  the  cube  at  45°  slightly  less.  The  rest  of  the  elements 
tend  to  be  fairly  close  to  one  another  with  the  large  grooved  element  having  the  smallest 
vorticity  maximum. 


At  the  x/d=10.0  location,  the  vortex  structures  produced  by  the  cube  elements  still 
maintain  a  significant  amount  of  strength,  whereas  all  other  elements’  vortex  structures 
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have  dissipated  greatly.  These  axisymmetric  elements  experience  streamwise  vorticity 
nearly  an  order  of  magnitude  less  than  the  upstream  location.  This  tells  us  that  due  to  the 
sharp  edges  on  the  cube  elements,  there  is  a  significantly  stronger  vortex  structure 
produced  that  has  the  ability  to  maintain  its  strength  greater  streamwise  distances  within 
the  turbulent  boundary  layer.  Due  to  the  fact  that  the  vortices  produced  by  the  elements 
are  located  very  close  to  the  wall  (Table  3.1),  viscous  forces  will  have  a  major 
contribution  in  the  decay  of  the  vortices  as  they  propagate  downstream.  Another  aspect 
that  is  related  to  the  diffusion  and  transport  of  the  vortices  is  the  entrainment  of  fluid 
from  the  boundary  layer.  As  is  given  by  Cutler  and  Bradshaw  (1993b)  and  stated  in 
George  (2004),  the  equation  governing  the  transport  of  streamwise  vorticity  is, 


, ,  dC2x  dQ  hn. 
U^JL  +  V^-*-  +  W  - 


dx 


dy 


dz 


„  dU  _  dU  _  OU 
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(3.2) 

Terms  4  and  5  on  the  right  hand  side  of  this  equation,  the  gradients  of  the  normal  and 
shear  stresses,  are  the  main  contributors  to  the  diffusion  of  the  vortices. 


Table  3.2  presents  the  change  in  location  of  the  primary  vortex  structure  as  it  proceeds 
downstream.  The  Azc+  value  is  the  normalized  change  of  the  peak  Qx  vorticity  location  in 
the  spanwise  direction  whereas  the  Ayc+  value  is  the  normalized  change  in  the  normal-to- 
wall  direction.  All  values  are  differences  from  the  upstream  measurement  plane  x/d=2.75 
(1.36  fine  grooved)  and  the  downstream  x/d=10.0  measurement  plane. 


Table  3.2:  Movement  of  primary  vortex 


Element 

Cube  at  45° 

-39.1 

30.4 

Cube  at  90° 

-55.1 

24.1 

Cone 

-80.7 

33.7 

Gaussian 

-80.8 

31.6 

-16.5 

30.1 

-133.7 

43.4 

Fine  Grooved 

-86.2 

38.1 
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As  the  vortex  structures  proceed  downstream,  they  diffuse  and  are  transported  away  from 
the  wall  and  centerline.  The  movement  and  strength  of  the  vortex  structure  will  be 
affected  a  great  deal  by  the  generation  of  secondary  vorticity.  In  this  region  the  near  wall 
viscous  forces  will  cause  the  loss  of  strength  in  the  primary  vortex.  All  vortex  structures 
move  a  very  comparable  distance,  -24-43  wall  units,  in  the  y-direction  away  from  the 
wall.  The  vortex  transport  in  the  spanwise  z-direction  for  all  elements  is  not  nearly  as 
comparable  as  it  was  in  the  y-direction.  The  peaked  elements’  vortices  all  move  a  similar 
distance  except  for  the  large  grooved  element  whose  initial  vortex  structure  was  already  a 
significant  distance  from  the  centerline  at  the  upstream  measurement  location  (x/d=2.75). 
The  cube  elements  also  have  a  Azc+  value  that  is  reasonably  comparable  with  respect  to 
one  another.  The  peak  Qx  location  produced  by  the  cube  at  45°  does  not  move  as  much 
in  the  spanwise  direction  as  does  the  vortex  produced  by  the  cube  at  90°.  The  peak 
location  of  for  the  hemisphere  moves  the  greatest  distance  in  the  spanwise  direction, 
as  it  is  at  least  55%  greater  than  all  other  elements. 


3.4.2  Circulation 

The  circulation,  T,  was  calculated  by  integrating  the  velocity  vector  along  a  closed 
rectangular  circuit  in  the  y-z  plane  using  the  following  equation. 


where  V  is  the  velocity  vector  and  ds  is  an  element  of  the  circuit  C.  The  centerline 
(z=0)  and  the  wall  (y=0)  are  two  sides  of  the  circuit.  The  other  two  sides  are  formed  by 
y=constant  and  z=-constant  values  where  the  boundary  layer  shows  no  effect  of  the 
presence  of  the  roughness  element.  Due  to  the  fact  that  W=0  at  all  planes  along  the  path 
of  integration  (flow  symmetry  at  z=0,  no-slip  condition  at  the  wall,  and  two-dimensional 
flow  at  the  two  remaining  planes),  the  circulation  is  only  dependent  on  the  V  velocity 
profile  at  the  centerline.  Table  3.3  gives  the  circulation  of  the  vortex  produced  by  each 
element  at  successive  streamwise  locations.  This  table  is  presented  for  raw  comparisons 
between  the  respective  roughness  elements  only. 
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Cutler  and  Bradshaw  (1986)  have  defined  a  non-dimensional  parameter  involving  the 
circulation  of  a  vortex  as  Y/{UeS),  where  8  is  the  boundary  layer  thickness  and  Ue  refers 

to  the  free-stream  speed  of  the  flow.  According  to  this  notation,  vortices  in  the  range  of 
0.08  to  2  are  considered  weak.  For  the  vortices  generated  by  the  elements  under  current 
investigation,  the  largest  value  obtained  using  the  Cutler  and  Bradshaw  notation  is  0.009. 
Consequently,  the  vortices  under  scrutiny  are  much  weaker  than  those  defined  by  Cutler 
and  Bradshaw.  Even  though  the  present  vortices  are  an  order  of  magnitude  weaker  than 
those  defined  as  ‘weak’  by  Cutler  and  Bradshaw,  they  still  have  a  very  significant  effect 
on  altering  the  structure  of  the  turbulent  boundary  layer  (especially  in  the  near  wall 
region). 


Table  3.3:  Circulation,  T,  as  a  function  of  streamwise  distance,  xl d 


T  (m2/s) 

Element 

x/d=1.36 

x/d=2.75 

x/d=5.0 

x/d=10.0 

x/d=20.0 

x/d=3 1 .0 

x/d=40.0 

Cube  at  45° 

5.43E-03 

6.19E-03 

3.27E-03 

1.22E-03 

3.28E-04 

— 

N/A 

Cube  at  90° 

2.39E-03 

2.64E-03 

1.41E-03 

2.18E-04 

— 

N/A 

Cone 

9.67E-03 

6.53E-03 

1.82E-03 

8.53E-04 

4.23E-04 

N/A 

2.57E-04 

Gaussian 

8.23E-03 

4.37E-03 

1.66E-03 

9.55E-04 

4.65E-04 

N/A 

2.97E-04 

6.27E-03 

2.72E-03 

1.50E-03 

9.57E-04 

2.58E-04 

N/A 

5.56E-05 

4.07E-03 

3.12E-03 

4.02E-04 

9.86E-05 

4.42E-05 

N/A 

Fine  Grooved 

5.18E-03 

2.92E-03 

1.67E-03 

6.77E-04 

3.38E-04 

N/A 

— 

It  is  readily  identifiable  that  the  two  elements  with  the  ‘sharpest’  peaks  (cone  and 
Gaussian)  have  the  highest  amount  of  circulation  at  the  x/d=1.36  location.  By  a 
downstream  distance  of  x/d=40.0,  the  circulation  has  more  or  less  dissipated  for  all 
elements.  Please  note  that  a  dashed  line  in  Table  3.3  refers  to  a  circulation  value  of  zero. 
In  each  of  these  instances,  the  local  V  velocity  profile  has  collapsed  to  the  reference  2D 
profile  thus  giving  rise  to  a  zero  value  of  circulation.  Figure  3.25  shows  the  circulation, 
T,  non-dimensionalized  by  UtJaJ  as  a  function  of  streamwise  distance  x/d.  The  cone 

and  Gaussian  elements  do  have  the  highest  circulation  values,  but  they  also  drop  off 
fairly  rapidly  and  meet  up  with  the  other  elements  around  x/d=5.0  for  the  Gaussian 
element  and  at  x/d=10.0  for  the  cone. 
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The  cube  elements  do  not  have  the  same  near  element  trend  as  the  axisymmetric 
elements.  The  x/d=1.36  location  for  the  cubes  shows  a  lower  value  of  circulation  than  is 
present  at  the  x/d=2.75  location.  This  phenomenon  can  be  explained  through  brief 
analysis  of  the  V+  plots  shown  in  Figure  3.12.  At  the  x/d=1.36  location  there  is  a 
significant  amount  of  downwash  near  the  element  height  only.  This  is  due  to  the 
blockage  that  the  cube  elements  present  to  the  flow.  The  measurement  location  is 
relatively  close  to  the  cube  and  there  is  still  a  significant  region  of  lower  velocity  fluid 
present.  The  x/d=2.75  profile  shows  significant  downwash  through  the  boundary  layer 
profile  as  well  as  above  the  height  of  the  element.  This  downstream  location  not  only 
sees  more  fluid  coming  over  the  top  of  the  cube  elements  that  is  being  directed  toward  the 
wall,  but  also  more  of  an  effect  from  the  horseshoe  vortex  structures  that  are  present. 
Subsequently,  this  larger  flow  magnitude  can  be  seen  in  the  high  circulation  values  over 
the  rest  of  the  measurement  locations.  This  effect  can  be  seen  in  particular  with  the  cube 
at  45°  whose  value  of  circulation  is  the  greatest  for  all  locations  between  x/d=2.75  and 
x/d=10.0.  All  elements  show  a  sharp  drop  in  circulation  until  x/d~10.0.  From  this 
location  there  is  a  more  gradual  decrease  in  magnitude  as  the  circulation  dissipates  and 
eventually  reaches  zero.  The  decay  of  circulation  can  primarily  be  attributed  to  the 
turbulent  diffusion  of  the  vortex  structure  as  it  propagates  downstream  as  well  as  the  near 
wall  viscous  forces  in  the  spanwise  direction  that  act  in  retarding  the  flow.  Viscous 
diffusion  is  also  a  contributor  as  the  vortex  structure  is  located  in  the  near-wall  region  of 
the  turbulent  boundary  layer. 

3.4.3  Wall  Shear  Stress 

For  information  on  the  scheme  used  to  determine  the  friction  velocity,  UT,  please  refer  to 
Section  2.3:  Post  Processing.  In  all  figures  shown,  the  wall  shear  stress,  zw  =U^p,  is 
presented  as  a  ratio  of  the  local  wall  shear  value  to  the  reference  2D  smooth  wall 
boundary  layer,  x0.  Figure  3.26  shows  the  variation  of  the  normalized  wall  shear  (xw/x0) 
along  the  centerline  whereas  Figures  3.27  and  3.28  show  the  spanwise  variation  in  wall 
shear  at  the  x/d=2.75  and  10.0  planes,  respectively.  As  is  discussed  in  Appendix  C,  the 
uncertainty  in  UT  is  ±5%.  Consequently,  it  is  important  to  bear  in  mind  that  a  small 
percentage  of  the  variations  in  the  presented  wall  shear  stress  data  is  due  to  uncertainties. 


Chapter  3:  Results 


-61  - 


Focusing  on  Figure  3.26,  the  centerline  variation  in  the  wall  shear  stress,  the  first 
measurement  location  after  reattachment  for  all  elements  is  x/d=2.75.  Therefore,  this 
location  is  the  first  available  position  at  which  local  xw  values  are  available.  All 
axisymmetric  elements  are  very  near  the  same  value  of  xw  at  x/d=2.75,  whereas  the  cube 
elements  have  a  drastically  lower  value.  The  lower  wall  shear  stress  values  exhibited  by 
the  cubes  are  due  to  the  fact  that  the  sharp  edges  on  the  cubes  create  a  defined  separation 
point.  This  separation  leads  to  an  area  of  lower  momentum  fluid  directly  behind  the 
elements.  Furthermore,  this  location  is  just  aft  of  flow  reattachment  and  is  in  a  newly 
formed  shear  layer  that  extends  downstream  of  the  reattachment  point.  All  elements 
show  a  distinct  peak  of  high  local  wall  shear  values  at  x/d=5.0.  All  ‘peaked’  elements  are 
located  at  nearly  identical  values  with  the  hemisphere  slightly  below  and  the  cube 
elements  at  the  lowest  value.  This  point  of  peak  wall  shear  is  caused  by  the  horseshoe 
vortex  structures  entraining  outer  layer  fluid  and  creating  a  wall-ward  rush  of  fluid.  A 
significant  drop  can  be  seen  in  the  values  between  x/d=5.0  and  10.0,  but  the  wall  shear 
value  is  still  larger  than  the  reference  smooth  wall  value.  This  increased  wall  shear  with 
respect  to  the  reference  smooth  wall  is  seen  through  the  remainder  of  the  measurement 
locations  until  x/d=40.  At  this  location,  all  elements  seem  to  have  collapsed  back  to  the 
reference  wall  shear  value.  This  indicates  that  the  horseshoe  vortices  exhibit  their 
presence  in  the  boundary  layer,  albeit  very  small,  for  a  significant  distance. 

Focusing  attention  on  the  spanwise  wall  shear  stress  variation  at  the  x/d=2.75  plane, 
Figure  3.27,  we  see  that  all  elements  exhibit  a  peak  value  of  wall  shear  at  the  centerline 
location.  Again  this  is  due  to  the  outer  layer  fluid  being  entrained  and  brought  toward  the 
wall,  creating  an  increase  in  wall  shear  relative  to  the  reference  wall  shear  value.  In  all 
cases  away  from  the  centerline  at  zAHt0.28,  there  is  a  decrease  in  wall  shear  which  leads 
to  a  minimum  value.  From  this  minimum  value,  the  wall  shear  gradually  increases  until 
around  z/d— 1.75-2.0.  In  this  region  of  increasing  wall  shear  the  secondary  flow  vortices 
are  nearly  parallel  to  the  wall.  Once  the  spanwise  edge  of  the  wake  is  reached,  the  wall 
shear  once  again  reaches  the  smooth  wall  value  within  the  experimental  uncertainty  of 
±5%.  The  region  of  lower  wall  shear  is  related  to  the  flow  being  directed  away  from  the 
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wall,  as  is  also  visible  in  the  oil  flow  images  previously  discussed.  The  cube  elements 
show  a  significantly  lower  wall  shear  value  as  compared  to  the  rest  of  the  elements  due  to 
the  previously  discussed  reasons.  There  is  a  sharp  change  in  the  wall  shear  stress  value 
from  z/d=-0.55  to  z/d=-0.83.  This  dramatic  change  is  due  to  the  effects  caused  by  the 
sharp  edges  of  the  cubes.  Once  past  the  sides  of  the  cubes,  the  wall  shear  reaches  similar 
values  as  the  other  elements. 

As  the  flow  proceeds  downstream  to  the  x/d=10.0  plane,  there  is  still  fluid  being  brought 
toward  the  wall  via  the  counter-rotating  horseshoe  vortex  structure.  For  all  of  the 
axisymmetric  elements  we  see  similar  tendencies  in  the  wall  shear  as  we  saw  at  the 
x/d=2.75  plane.  Although  the  tendencies  are  similar,  the  magnitudes  have  been  reduced 
due  to  the  diffusion  of  the  horseshoe  vortex  structure  as  it  progresses  downstream. 
Again,  there  is  a  peak  of  higher  wall  shear  due  to  the  downwash  of  fluid  toward  the  wall. 
Directly  on  either  side  of  this  peak  we  see  a  minimum  value  and  then  an  increase  to  the 
reference  2D  flow  wall  shear  value.  As  the  outer  regions  of  the  spanwise  direction  are 
reached  there  is  a  noticeable  difference  between  the  wall  shear  values  of  the  elements. 
This  difference  can  be  attributed  to  the  differences  in  the  horseshoe  vortex  structures 
produced  by  each  element.  At  the  z/d=±0.28  locations,  the  large  grooved  element  is  seen 
to  have  a  different  wall  shear  value  as  compared  to  the  other  axisymmetric  elements. 
This  difference  is  attributed  to  the  weaker  vortex  that  is  present,  at  the  x/d=10.0  location, 
in  the  large  grooved  element  case  as  is  shown  in  Table  3.1.  There  is  a  noticeable 
difference  in  the  wall  shear  behavior  of  the  cube  elements  around  the  centerline  when  the 
two  measurement  planes  are  compared.  The  x/d=10.0  plane  exhibits  a  rise  in  wall  shear 
at  the  closest  off-centerline  locations  and  then  a  decrease  to  near  smooth  wall  values. 
The  slightly  lower  values  at  the  centerline  for  the  cube  cases  can  be  related  to  the  fact  that 
the  flow  diverges  from  the  plane  of  symmetry  in  the  near  wall  regions.  This  divergence 
in  turn  causes  the  reduction  in  wall  shear  at  the  centerline  location.  As  is  evidenced  by 
Tables  3.1  and  3.2,  the  vortices  created  by  the  cube  elements  have  not  moved  as  far  in  the 
spanwise  direction  as  the  other  elements.  This  lack  of  movement  creates  a  different 
effect  on  the  centerline  wall  shear  value  as  more  fluid  diverges  away  from  the  centerline 
at  the  downstream  measurement  plane.  This  phenomenon  is  in  contrast  to  the  greater 
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amount  of  fluid  converging  toward  the  centerline  at  the  upstream  measurement  plane  thus 
creating  a  larger  wall  shear  stress  at  the  centerline  measurement  location. 


3.4.4  Drag 

The  drag  for  each  element  was  calculated  using  a  control  volume  and  momentum  balance 
approach.  The  undisturbed  upstream  velocity  profile  as  well  as  the  wake  profiles  behind 
each  element  at  the  x/d=2.75  plane  are  taken  into  account  in  this  analysis.  The  derivation 
of  the  drag  equation  and  related  discussion  can  be  found  in  Appendix  B.  Equation  3.4 
below  was  used  to  calculate  the  drag  increment  due  to  the  presence  of  a  roughness 
element  and  is  taken  from  Equation  (B.ll)  in  Appendix  B.  The  subscripts  ‘E’  and 
‘w/oE’  refer  to  boundary  layer  profiles  taken  with  and  without  the  element  being  present 
in  the  boundary  layer,  respectively.  The  numbered  subscripts  refer  to  the  control  volume 
plane  number  at  which  the  profile  was  taken.  Similarly,  the  subscript  ‘top’  refers  to  the 
quantity  at  the  top  of  the  control  volume.  Finally,  -z2d  is  related  to  the  spanwise  distance 
in  the  (-)z-direction  at  which  the  local  boundary  profile  has  collapsed  to  the  reference  2D 
smooth  wall  profile. 


2A De 


incremental 


(3.4) 


The  incremental  drag  consists  of  the  drag  due  to  the  element’s  presence  in  the  boundary 
layer  as  well  as  the  skin  friction  drag  created  by  the  surface  area  of  the  floor  contained 
within  the  control  volume.  The  control  volume  used  in  the  derivation  of  Equation  3.4 
deems  the  contribution  of  the  pressure  term  to  form  drag  to  be  negligible  since  the 
upstream  and  downstream  control  volume  planes  should  be  at  the  same  uniform  pressure. 
This  result  stems  from  the  work  of  Matinuzzi  and  Tropea  (1993)  who  show  a  pressure 
recovery  at  a  downstream  distance  of  x/d=2.25  for  a  surface  mounted  cube.  Thus,  due  to 
the  roughness  element  geometries  under  investigation,  all  elements  should  experience  full 
pressure  recovery  by  x/d=2.75.  Please  note  that  the  analysis  used  to  calculate  the  drag 

also  includes  the  streamwise  Reynolds  normal  stress,  u2 ,  term.  Consequently,  this 
approach  is  more  rigorous  than  a  method  that  only  takes  into  account  the  mean  velocity. 
The  addition  of  the  fluctuating  term  will  cause  a  slight  decrease  in  the  drag  as  compared 
to  the  drag  calculated  with  only  the  mean  velocity.  This  is  due  to  the  fact  that  some  mean 
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flow  momentum  is  lost  when  turbulence  is  created.  As  was  mentioned  previously,  it  is 
necessary  to  account  for  the  drag  on  the  control  volume  surface  in  order  to  obtain  the 
actual  drag  due  to  each  roughness  element.  More  discussion  related  to  the  procedure 
used  to  obtain  the  actual  drag  can  be  found  in  Chapter  4.  Table  3.4  presents  the  actual 
drag  calculated  on  each  element  as  well  as  the  related  elemental  areas,  where  Af  and  Aw 
correspond  to  the  frontal  (Af)  and  wetted  (Aw)  areas  of  each  roughness  element.  There 
are  no  drag  results  presented  for  the  fine  grooved  element  due  to  the  fact  that 
measurements  were  taken  at  a  plane  equal  to  x/d=l  .36  rather  than  2.75. 


Table  3.4:  Actual  experimental  drag  and  elemental  areas 


Element 

Af  (m2) 

Aw  (m2) 

Aw/Af 

Cone 

3.23E-06 

1.13E-05 

3.51 

1.69E-04 

52.49 

Gaussian 

2.66E-06 

1.02E-05 

3.85 

1.22E-04 

45.81 

3.23E-06 

1.94E-05 

6.01 

1.98E-04 

61.62 

Cube  at  90° 

2.98E-06 

1.49E-05 

5.00 

2.52E-04 

84.63 

Cube  at  45° 

4.22E-06 

1.49E-05 

3.54 

2.13E-04 

50.64 

3.07E-06 

1.23E-05 

4.00 

1.15E-04 

37.70 

Fine  Grooved 

3.23E-06 

2.12E-05 

6.58 

N/A 

N/A 

The  largest  value  of  drag  is  seen  on  the  cube  at  90°  followed  by  the  cube  at  45°,  large 
grooved  element,  cone,  Gaussian  element,  and  the  hemisphere.  The  higher  value  of  drag 
on  the  large  grooved  element  can  be  related  not  only  to  its  frontal  projected  area  but  also 
its  wetted  surface  area.  With  the  highest  wetted  surface  area,  next  to  the  fine  grooved 
element,  the  large  grooved  element  will  experience  more  viscous  drag  than  the  other 
elements  under  scrutiny.  This  result  is  somewhat  obvious  and  is  shown  by  Table  3.4. 
The  drag  on  each  roughness  element  appears  to  be  in  direct  correlation  with  a 
combination  of  the  form  drag,  related  to  the  frontal  area,  and  the  viscous  drag,  which  can 
be  related  to  the  wetted  surface  area,  of  each  element.  More  discussion  related  to  the 
drag  of  each  element  and  possible  correlations  between  the  roughness  elements  is 
discussed  in  Chapter  4. 

3.4.5  Eddy  Viscosities  and  Anisotropy  Factor 

The  eddy  viscosities  in  the  x  and  z-directions  were  computed  using  Equations  3.5  and  3.6 
respectively. 
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_  -uv 
S*  ~  dUJdy 


(3.5) 


-  vw 
1  ~  dW/dy 


(3.6) 


The  centerline  variation  of  normalized  streamwise  eddy  viscosity,  sx/v  vs.  'v  can 

be  found  in  Figure  3.108.  The  eddy  viscosity  is  an  indicator  of  how  much  of  the  mean 
flow  gradient  goes  into  the  production  of  turbulence.  In  the  near  wall  region  all  profiles 
behave  in  the  same  fashion  until  y+~15.  In  this  region  there  is  an  increase  in  the 
streamwise  eddy  viscosities  behind  the  elements  due  to  the  high  amount  of  turbulence 
present  in  this  region.  At  the  closest  measurement  location  to  the  element,  x/d=1.36, 
there  are  a  few  points  that  have  negative  eddy  viscosities.  These  events  can  be  linked  to 
the  fact  that  there  is  still  a  significant  amount  of  backflow  in  this  area  behind  the  element. 
The  increase  in  the  streamwise  eddy  viscosities  can  be  seen  until  just  above  the  roughness 
element  heights,  y+~150,  where  the  profiles  tend  to  collapse  back  to  the  reference  smooth 
wall  values. 


The  parameter  N,  is  the  anisotropy  factor  of  a  flow  and  is  defined  as, 

n=£j_  =  -vw/(dW/dy)  _  tan{SSA) 

£x  -uv/(dU/dy)  tan  (FGA) 

where  SSA  and  FGA  are  the  shear  stress  angle  and  flow  gradient  angle.  Further 
discussion  related  to  these  parameters  can  be  found  in  Appendix  E.  Figures  3.109 
through  3.112  show  contours  of  N  in  the  different  measurement  planes  associated  with 
each  element.  It  is  largely  assumed  that  in  three  dimensional  boundary  layer  flows  that 
the  eddy  viscosity  is  isotropic.  This  assumption  is  generally  not  valid  and  is  also  shown 
to  be  the  case  in  the  present  study.  In  order  for  the  flow  to  be  isotropic  the  eddy  viscosity 
ration  N  would  be  equal  to  one  and  any  deviation  from  unity  would  mean  that  the  flow  is 
anisotropic.  As  is  shown  in  Figures  3.109  through  3.1 12,  the  flow  is  largely  anisotropic 
and  any  values  larger  than  unity  indicate  a  more  significant  spanwise  interaction  as 
compared  to  the  streamwise  direction.  Exceedingly  large  values  in  the  magnitude  of  N  in 

the  plots  should  not  be  taken  as  actual  results  due  to  the  scatter  in  the  vw  Reynolds 
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shearing  stress  at  theses  particular  locations.  It  is  also  important  to  note  that  any  negative 
regions  of  N  can  be  attributed  to  the  decorrelation  of  the  streamwise  and  spanwise 
Reynolds  shearing  stresses. 


3.5  Reynolds  Stresses  and  Turbulent  Kinetic  Energy 

An  effective  means  of  analyzing  the  turbulent  flow  field  behind  the  elements  is  through 
the  inspection  of  the  Reynolds  stress  behavior.  To  this  end,  the  centerline  variation  of  the 

normalized  streamwise  Reynolds  normal  stress,  u2  jl)] ,  plotted  against  yUr/v  is  shown 

in  Figure  3.29.  At  the  x/d=1.36  location,  all  elements  show  a  peak  in  the  streamwise 
Reynolds  normal  stress  at  or  slightly  below  the  respective  element  height.  All 

axisymmetric  elements  are  shown  to  have  a  greater  u 2  Reynolds  normal  stress 
throughout  the  early  stages  of  the  profile  as  compared  to  the  reference  smooth  wall 
boundary  layer.  These  elements  then  show  a  sharp  drop  off  in  the  streamwise  Reynolds 
normal  stress  and  collapse  to  the  reference  smooth  wall  profile  above  the  height  of  the 
elements  near  y+~200.  In  contrast  to  the  axisymmetric  elements,  the  cube  elements  show 

a  reduction  in  u 2  near  the  wall.  This  reduction  can  be  attributed  to  the  blockage  that  the 
cubes  present,  therefore  limiting  the  amount  of  mixing  that  occurs  directly  behind  the 
element.  The  higher  values  of  normal  stresses  are  due  to  the  intense  mixing  created  by 
the  shear  layers  that  are  rolling  up  and  converging  toward  the  centerline  of  the  element. 

The  cube  at  45°  shows  the  highest  magnitude  of  u2  as  compared  to  the  rest  of  the 
elements.  The  cube  at  45°  is  followed  by  a  close  conglomeration  of  the  hemisphere,  fine 

grooved,  and  cone  element  data.  The  lowest  value  of  u2  at  the  earliest  measurement 
location  is  seen  in  the  Gaussian,  large  grooved,  and  cube  at  90°  elements.  At  the 

x/d=2.75  location  the  maximum  value  in  u2  is  still  from  the  cube  at  45°;  the  other 
elements  have  collapsed  to  approximately  the  same  peak  value.  The  peak  in  the 
streamwise  Reynolds  normal  stress  is  still  seen  at  or  below  the  element  heights  with  a 
collapse  to  the  reference  boundary  layer  at  approximately  the  same  location  as  the 

x/d=1.36  measurement  location.  The  reduction  in  peak  magnitude  in  u 2  of  about  20- 
60%  is  very  large  as  the  flow  propagates  downstream  to  the  x/d=2.75  location.  As  the 
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flow  continues  to  proceed  farther  downstream,  there  is  a  small  reduction  in  the 
streamwise  Reynolds  normal  stress  due  to  slight  flow  acceleration  along  the  centerline 
behind  the  elements.  By  x/d=T0.0  the  profiles  have  nearly  collapsed  to  the  reference 
boundary  layer  profile  and  at  x/d=20.0  and  40.0  there  is  no  effect  of  the  roughness 
elements’  presence  in  the  boundary  layer. 

Contour  plots  of  the  normalized  streamwise  Reynolds  normal  stress  at  the  x/d=2.75  and 
10.0  measurement  planes  can  be  found  in  Figures  3.30  through  3.33.  These  plots  show 

that  for  all  elements  there  is  a  peak  value  in  u2  at  or  slightly  below  the  height  of  the 
element.  All  elements  show  that  this  peak  value  is  also  connected  to  a  higher  magnitude 

pocket  located  to  the  side  of  each  element,  z/d~l.  This  high  region  of  the  u 2  Reynolds 
stress  can  be  attributed  to  the  shear  layers  that  are  rolling  up  and  surrounding  the  trailing 
leg  of  the  horseshoe  vortex  structure.  Consequently,  the  vertical  locations  of  the  pockets 

of  high  u2  Reynolds  stress  are  similar  to  the  peak  Qx  locations  shown  in  Table  3.1.  The 
distance  from  the  wall  for  the  peak  values  of  w2  is  near  y+;=20  whereas  the  peak  Qx 

locations  were  found  to  be  near  y+~13.  The  magnitude  of  this  u 2  region  is  related  to  the 
magnitude  of  the  vortex  structure  produced  by  each  element,  in  this  case  the  cube  at  45° 
is  the  highest  followed  by,  cube  at  90°,  cone,  hemisphere,  Gaussian,  and  large  grooved 
elements.  At  the  x/d=10.0  plane  the  values  for  all  elements  are  nearly  equal  to  the 
reference  two-dimensional  boundary  layer  values.  There  is  some  reduction  near  the 
center  plane  due  to  slight  acceleration,  but  for  the  most  part  the  influence  of  the  element 
has  drastically  dissipated. 

Figure  3.34  shows  the  centerline  variation  of  the  normalized  normal-to-wall  Reynolds 
normal  stress,  v2/u2  ,  plotted  with  respect  to  yUjv .  The  same  trends  can  be  seen  in 

these  plots  as  was  seen  in  the  plots  of  u2 jl)2  .  There  is  a  peak  value  at  or  below  the 
element  heights,  with  the  cube  at  45°  generally  having  the  largest  magnitude.  The  cube 

at  90°  at  the  x/d=1.36  location  shows  a  significantly  lower  value  of  v2  again  due  to  the 
lower  speed  fluid  created  by  the  ‘blockage  effect’.  The  elements  drop  off  sharply  and 
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collapse  to  the  reference  boundary  layer  values  just  above  the  height  of  the  element  near 
y+~200.  At  x/d=2.75,  lower  magnitudes  in  the  v2  stress  are  seen  as  compared  to  the 
upstream  location,  but  there  is  not  as  much  of  a  decrease  in  the  v2  stress  as  was  seen  in 

the  m2  stress.  By  x/d=5.0  it  appears  that  the  cube  element  profiles  have  collapsed  nearly 
to  one  another,  while  the  axisymmetric  elements  have  done  the  same.  At  x/d=10.0, 
profiles  for  all  elements  have  collapsed  nearly  to  the  same  values  and  by  x/d=20.0  and 
40.0  there  is  no  influence  by  the  elements  seen  in  the  profiles. 

Contour  plots  of  the  normalized  v2  Reynolds  normal  stress  can  be  seen  in  Figures  3.35 
through  3.38  for  the  x/d=2.75  and  10.0  planes,  respectively.  At  the  x/d=2.75  plane, 
nearly  circular  contours  surround  a  peak  stress  value  that  is  located  slightly  below  the 
height  of  the  element.  These  peaks  in  the  Reynolds  normal  stress  can  again  be  attributed 
to  the  mixing  created  by  shear  layer  roll-up  due  to  the  horseshoe  vortex  structures.  The 
rush  of  fluid  toward  the  wall  brought  in  by  the  common  downwash  of  the  vortices  is  the 
primary  factor  that  will  create  and  maintain  these  higher  stress  values.  Consequently,  the 

highest  regions  of  the  v2  Reynolds  normal  stress  are  found  near  the  height  of  each 

element.  The  ‘peaked’  roughness  elements  tend  to  see  the  largest  values  in  the  v2  stress 
at  a  y-location  near  y+~70,  whereas  the  hemisphere  and  cube  elements  are  closer  to  the 

wall  at  y+=40-50.  These  locations  differ  from  the  locations  of  the  peak  u2  stress  and 
peak  streamwise  vorticity,  Qx.  in  that  they  are  along  the  centerline  and  significantly 
farther  from  the  wall.  The  contour  plots  at  the  x/d=10.0  plane  show  a  small  peak  of  stress 
that  is  still  located  near  the  element  heights.  This  area  is  of  considerably  less  magnitude 
than  the  upstream  location.  The  higher  values  tend  to  be  stretched  in  the  spanwise 
direction  as  the  vortices  move  away  from  the  centerline  and  dissipate. 

Figure  3.39  shows  the  centerline  variation  of  the  normalized  spanwise  Reynolds  normal 
stress,  w2/u2 ,  plotted  with  respect  to  yUjv .  These  centerline  plots  are  very 

comparable  and  show  nearly  the  same  trends  as  do  the  u2 /u2T  and  v2/f/2  plots.  A  few 
slight  differences  are  that  the  w2  stresses  are  larger  in  the  near-wall  region,  they  keep 
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their  magnitude  for  a  wider  range  of  y  values,  and  the  peak  values  are  seen  farther  below 
the  element  heights  as  compared  to  the  other  stresses.  Also,  there  are  no  local  peaks  in 

w2  off  of  the  centerline  at  x/d=2.75.  Contour  plots  of  the  normalized  w2  Reynolds 
normal  stress  can  be  seen  in  Figures  3.40  through  3.43  for  the  x/d=2.75  and  10.0  planes, 
respectively.  These  plots  show  the  same  tendency  for  the  peak  stress  value  below  the 

element  height  at  the  x/d=2.75  plane.  There  is  also  a  higher  region  of  w2  stress  off  of 
the  centerline  in  the  spanwise  direction  at  the  x/d=10.0  plane  that  can  be  linked  to  the 
propagation  of  the  vortices  in  the  spanwise  direction  as  they  dissipate  and  proceed 
downstream. 

The  normalized  Reynolds  shear  stress  profiles  of  -  uv/u r2  along  the  centerline  can  be 

found  in  Figure  3.44.  These  plots  show  a  significant  peak  at  the  x/d=1.36  location  that  is 
an  order  of  magnitude  larger  for  the  cube  at  45°,  hemisphere,  and  fine  grooved  element. 
As  the  flow  proceeds  downstream  to  the  x/d=2.75  location,  the  cube  elements  still  show 

the  largest  Reynolds  shear  stress  values,  whereas  the  -  uv  stress  for  the  other  elements 
has  deteriorated  to  nearly  the  same  values.  Overall  there  is  a  reduction  in  the  Reynolds 
shear  stress  value  in  the  range  of  40-50%.  The  profiles  for  all  elements  have  almost 
collapsed  to  one  another  by  x/d=5.0  with  the  cubes  being  the  outliers.  At  x/d=10.0  the 
centerline  plots  show  no  evidence  of  the  roughness  element  as  all  the  profiles  have 
collapsed  to  the  reference  profile. 

Contour  plots  of  the  normalized  Reynolds  shear  stress,  -  uv ,  can  be  seen  in  Figures  3.45 
through  3.48  for  the  x/d=2.75  and  10.0  planes,  respectively.  The  x/d=2.75  plane  shows 
large  peaks  near  the  element  heights.  This  peak  region  is  due  to  the  ability  of  the  vortex 
structures  to  create  significant  sweeping  motions  that  promote  mixing.  A  lower  region  of 
shear  stress,  as  compared  to  the  peak  region,  is  also  seen  extending  away  from  the 
elements  in  the  spanwise  direction.  By  the  time  the  flow  has  reached  the  x/d=10.0  plane, 
the  Reynolds  shear  stress  has  undergone  very  substantial  decay  and  is  very  close  to  two- 
dimensional  values  with  small  regions  of  higher  shear  stress  extending  in  the  spanwise 
direction. 
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The  plots  contained  in  Figures  3.49  through  3.102  are  the  centerline  and  contour  plots 

associated  with  the  uw  and  vw  Reynolds  shear  stresses  as  well  as  the  remaining  triple 
products.  These  figures  show  good  symmetry  for  all  plots  involving  a  single  spanwise 
fluctuating  term.  This  result  is  expected  due  to  the  fact  that  the  spanwise  (W)  velocity  is 
zero  along  the  centerline  of  the  element.  All  parameters  involving  a  single  normal-to- 
wall  fluctuating  term  tend  to  show  peak  values  along  the  centerline  and  near  the  height  of 
the  element.  This  is  the  area  that  the  most  significant  amount  of  downwash  is  seen. 
Other  areas  of  peak  triple  products  can  also  be  seen  in  the  high  spanwise  velocity  regions 
associated  with  the  secondary  flow  produced  by  the  horseshoe  vortex  structures. 
Consequently,  areas  near  the  element  heights  and  in  the  high  secondary  flow  regions  will 
have  significant  levels  of  triple  products  and  stresses  due  to  the  intense  mixing  and  the 
rolling  up  of  shear  layers  that  will  be  present. 

Plots  of  the  centerline  variation  of  the  TKE,  {cj2  / 2^1/ 2  ,  with  respect  to  y+  can  be  seen  in 

Figure  3.103.  The  highest  values  of  TKE  are  always  found  below  the  height  of  each 
element.  The  x/d=1.36  location  shows  the  cube  at  45°  with  the  highest  TKE  followed 
closely  by  all  of  the  axisymmetric  elements  and  then  the  cube  at  90°.  Significant  decay 
between  the  first  and  second  measurement  locations  can  also  be  seen  with  the  cube 
elements  now  both  having  the  highest  amount  of  TKE.  There  is  still  a  noteworthy 
amount  of  TKE  present  at  x/d=5.0  for  all  elements.  At  x/d=10.0  there  is  very  little  TKE 
present  as  was  discussed  previously.  By  the  time  the  flow  reaches  a  downstream  distance 
of  x/d=20.0  the  TKE  has  dissipated  and  reached  the  reference  two-dimensional  turbulent 
boundary  layer  values. 

Figures  3.104  through  3.107  show  contour  plots  of  the  turbulent  kinetic  energy  (TKE) 
normalized  on  UT,  ^  jl^LJ2  ,  where  q2  =^i2  +  v2  +w2j.  These  plots  contain  the 
contours  of  the  x/d=2.75  and  10.0  measurement  planes,  respectively.  These  plots  also 
show  the  normalized  TKE  transport  velocity  vectors,  Vq /UT  ={yq2 / q2  ]/t/r  and 
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Wq fU x  =  ^ I q2]lUT .  At  the  x/d=2.75  plane,  a  core  of  higher  TKE  can  be  found 

below  the  element  height  along  the  centerline.  The  large  regions  of  TKE  near  the 
element  heights  are  due  to  the  horseshoe  vortices  that  not  only  create  sweeping  motions 
toward  the  wall  but  also  a  convergence  of  associated  shear  layers  toward  the  center  plane. 
The  largest  magnitude  of  TKE  is  produced  by  the  cube  at  45°  followed  by  the  cube  at 
90°,  fine  grooved  (from  Figure  3.103  at  the  x/d=2.75  location),  cone,  large  grooved, 
Gaussian,  and  hemisphere.  The  TKE  tends  to  diffuse  outward  in  a  radial  pattern  above 
the  element  heights.  This  radial  pattern  can  only  be  observed  in  the  plots  having  log  y 
coordinates.  Behind  the  element  and  below  its  height  there  is  still  a  large  amount  of  TKE 
present.  The  larger  elements  that  protrude  into  the  boundary  layer  and  have  sharp  edges 
tend  to  produce  more  TKE.  The  hemisphere  which  is  shorter  and  has  smooth  sides  has 
the  lowest  amount  of  TKE  production. 

The  TKE  diffuses  outward  from  the  core  of  highest  TKE  as  is  evidenced  by  the  diffusion 
velocity  vectors.  The  highest  diffusion  is  seen  in  the  cone  and  cubes  whose  vectors  have 
very  similar  magnitudes.  The  Gaussian,  fine  grooved,  large  grooved,  and  hemisphere 
elements  follow  with  decreasing  diffusion  velocities.  The  basis  of  putting  the  fine 
grooved  element  ahead  of  the  large  grooved  element  is  due  to  the  element’s  ‘peakiness’. 
As  can  be  seen  in  Figures  3.104  and  3.105,  the  diffusion  velocities  are  highest  near  sharp 
edges  (cubes)  and  near  sharp  peaks  (cone,  Gaussian,  etc.).  Subsequently,  this  leads  to  the 
conclusion  that  both  the  element  shape  and  the  orientation  play  very  significant  roles  in 
the  production  and  diffusion  of  TKE.  All  axisymmetric  elements  show  a  higher  diffusion 
of  TKE  near  their  respective  heights  than  anywhere  else.  The  cube  elements  tend  to  have 
large  diffusion  along  the  sides  as  well  as  near  the  element  height.  It  is  also  worth  noting 
that  even  though  the  hemisphere  is  smooth,  there  are  still  significant  diffusion  velocities 
that  emanate  from  the  element  height. 

The  x/d=10.0  measurement  plane  (Figures  3.106  and  3.107)  shows  a  very  extreme  decay 
of  the  TKE  produced  by  each  element.  The  most  noticeable  regions  of  higher  TKE  are 
seen  away  from  the  elements  in  the  spanwise  direction.  These  areas  are  mainly  due  to  the 
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streamwise  vortices  that  have  propagated  downstream.  The  diffusion  velocity  vectors 
have  also  gone  through  a  drastic  decay  and  are  for  the  most  part  negligible. 
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3.6  Chapter  3  Figures:  Single  Element 


Figure  3.1:  Measurement  locations  for  single  element  cases  (except  fine  grooved) 


Figure  3.2:  Measurement  locations  for  fine  grooved  element 
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Figure  33:  Oil  flows  for  the  cube  and  hemisphere  elements 


Chapter  3:  Results _ -  75  - 


"O 


c 

lOO  > 
T—  O  O 
-  <1> 
“O  Q. 
<D  (/) 
>  0) 
O  u- 

2  £ 

o  | 

<D  ^ 
a*  to 
csi 


■O 

<1) 

> 

o 

o 

6 


o 

o 

CO 

G) 

CD 

c\i 


<1) 

ifs 

gti 

55  "o 
</)  ~ 
co  = 

3  CO 
CO 

O  £ 
d  £ 

£  5 

P  <N 


Figure  3.4:  Oil  flows  for  the  Gaussian,  fine  grooved,  large  grooved,  and  cone  elements 
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Figure  3.5:  Close-up  views  of  oil  flows  for  the  cube  and  hemisphere  elements 
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Figure  3.6:  Close-up  views  of  the  oil  flows  for  the  Gaussian,  large  grooved,  fine  grooved,  and  cone  elements 
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Figure  3.7: 


streamwise  mean  velocity  profiles  along  the  centerline 
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Figure  3.8:  Contours  of  U/UT  in  the  y-z  plane  at  a  location  of  x/d=2.75  (Note:  Fine  Grooved  x/d=1.36) 
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Figure  3.9:  Contours  of  U/UT  in  the  y-z  plane  at  a  location  of  x/d=2.75 
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Figure  3.10:  Contours  of  U/U T  in  the  y-z  plane  at  a  location  of  x/d=10.0 


Chapter  3;  Results _ -  82  - 


ooooooooooooo 

ppqpppppppppo 
cocoTrcviocoaD^tcvioodccTr 
CNJOJCNCNCN*—  *—*—*— 


OOOOOOOOOOOOO 

ooooooooooooo 


Q. 


o 


Hemisphere— k=1 .40  mm 

Figure  3.11:  Contours  of  UjUr  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.12: 


normal  to  wall  mean  velocity  profiles  along  the  centerline 
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Figure  3.13:  Contours  of  V jU T  in  the  y-z  plane  at  a  location  of  x/d=2.75  (Note:  Fine  Grooved  x/d=1.36) 
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Figure  3.14:  Contours  of  V/U T  in  the  y-z  plane  at  a  location  of  x/d=2.75 
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Figure  3.15:  Contours  of  V/U r  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.16:  Contours  of  V/Ut  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.17:  Contours  of  W/U T  in  the  y-z  plane  at  a  location  of  x/d=2.75  (Note:  Fine  Grooved  x/d=1.36) 
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Figure  3.18:  Contours  of  W/U r  in  the  y-z  plane  at  a  location  of  x/d=2.75 
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Figure  3.19:  Contours  of  W/Ur  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.20:  Contours  of  WjUT  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.21:  Contours  of  normalized  streamwise  vorticity,  QxyjA^  lu r ,  with  normalized  secondary  velocity  vectors  V/V T 
and  W jU r  in  the  y-z  plane  at  a  location  of  x/d=2.75  (Note:  Fine  Grooved  x/d=1.36) 


Figure  3.22:  Contours  of  normalized  streamwise  vorticity,  Qxyj~A~f  lu r ,  with  normalized  secondary  velocity  vectors  V/U r 

and  WjU r  in  the  y-z  plane  at  a  location  of  x/d=2.75 
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Figure  3.23:  Contours  of  normalized  streamwise  vorticity,  QxJaJ  ju  t  ,  with  normalized  secondary  velocity  vectors  V jU  r 

and  W/U r  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.24:  Contours  of  normalized  streamwise  vorticity,  Q.x^Af  jUT ,  with  normalized  secondary  velocity  vectors  V/U  T 

and  W/UT  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.25:  Variation  of  circulation,  r/(t/r^4^),  with  streamwise  distance,  x/d 
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Figure  3.28:  Variation  of  normalized  wall  shear  (xw/x0)  at  x/d=10.0  plane 
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Figure  3.29 
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Figure  3,43:  Contours  of  w2/f/r2  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.46:  Contours  of  uv/u r2  in  the  y-z  plane  at  a  location  of  x/d=2.75 
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Figure  3.56:  Contours  of  vw/u2  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.67:  Contours  of  v2w/u3t  in  the  y-z  plane  at  a  location  of  x/d=2.75 
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Figure  3.73:  Contours  of  uv2  jU]  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.74:  Contours  of  uv2  jU)  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.77:  Contours  of  uw1  jU]  in  the  y-z  plane  at  a  location  of  x/d=2.75 
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Figure  3.85:  Contours  of  uvwfu]  in  the  y-z  plane  at  a  location  of  x/d=2.75  (Note:  Fine  Grooved  x/d=1.36) 
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Figure  3.86:  Contours  of  uvw/u ]  in  the  y-z  plane  at  a  location  of  x/d=2.75 
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Hemisphene--k=1 .40  mm 

Figure  3.88:  Contours  of  uvwju ]  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.94:  v 3  vs. 


profiles  along  the  centerline 
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Figure  3.97:  Contours  of  v3/t/3  in  the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.103:  Turbulent  Kinetic  Energy  (TKE),  |^/2yt/r2J  vs. 
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the  y-z  plane  at  a  location  of  x/d=10.0 
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Figure  3.108:  Centerline  variation  of  normalized  streamwise  eddy  viscosity,  ex/v  vs. 
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Figure  3.111:  Contours  of  eddy  viscosity  ratio,  Ne  =  sjsx ,  in  the  y-z  plane  at  x/d=10.0 
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Figure  3.112:  Contours  of  eddy  viscosity  ratio,  Ne  =  sjsx ,  in  the  y-z  plane  at  x/d=10.0 
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3.7  Chapter  3  Figures:  Gaussian  Fetch 


Note; 

A  single  profile  at  x/d=1.08,  z/d=0.0  was  also  taken 
at  a  location  eight  (8)  cells  downstream  of  the  shown 
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Figure  3.113:  Measurement  locations  for  the  Gaussian  fetch  of  roughness 


Chapter  4  Discussion  of  Results 

4.1  Introduction 

The  purpose  of  this  chapter  is  to  unify  the  results  from  the  previous  chapters  into  a  useful 
and  applicable  understanding  of  the  experimental  results.  Schematics  of  the  turbulent 
flow  field  behind  each  roughness  element  are  presented  and  comparisons  are  made  with 
respect  to  each  element.  Work  associated  with  possible  correlations  between  the  drag  of 
all  elements  is  examined.  An  effort  related  to  the  implementation  of  a  circulation 
correlation  to  connect  an  element’s  circulation  to  its  respective  ‘bluntness  factor’  is  also 
presented.  The  streamwise  decay  of  circulation  for  each  element  is  discussed  along  with 
possible  relationships  to  all  tested  roughness  elements.  Through  this  analysis  and 
discussion,  a  more  complete  understanding  of  the  interrelation  between  each  flow  field 
and  its  associated  physics  will  be  brought  to  light.  It  is  of  significance  to  note  that  all 
experimental  data  acquired  throughout  this  investigation  meet  the  realizability  conditions 
set  forth  by  Schumann  (1977).  The  requirements  for  a  realizable  data  set  as  stated  by 

Schumann,  where  Rap  =  uaup  are  the  appropriate  Reynolds  stresses,  are  R«p>  0  for  a  = 

(3,  (Rap)  <  RaaRpp  for  a  f  (3,  and  the  matrix  R«p  be  positive  semi-definite.  Meeting  these 
requirements  signifies  that  the  data  are  realizable  in  the  current  coordinate  system,  as  well 
as  in  any  orthogonal  coordinate  system. 


4.2  Drag  Correlations 

As  was  discussed  briefly  in  Chapter  3,  the  drag  for  each  element  was  calculated  using  a 
control  volume  and  momentum  balance  approach.  The  differences  between  the 
undisturbed  upstream  velocity  profile  and  the  wake  profiles  at  a  downstream  location  of 
x/d=2.75  were  taken  into  account  throughout  the  investigation.  Thus,  the  fine  grooved 
element  is  the  only  element  not  analyzed  due  to  its  plane  of  measurements  being  taken  at 
x/d=1.36.  The  derivation  of  the  drag  equation  and  related  discussion  can  be  found  in 
Appendix  B.  Equation  4.1  below  was  used  to  calculate  the  drag  increment  due  to  the 
presence  of  a  roughness  element  and  is  taken  from  Equation  (B.l  1)  in  Appendix  B. 
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2A De 


incremental 


^ top  (^2  |  w /  0£ 


(4.1) 


The  incremental  drag  consists  of  the  drag  due  to  the  element’s  presence  in  the  boundary 
layer  as  well  as  the  skin  friction  drag  created  by  the  surface  area  of  the  floor  contained 
within  the  control  volume.  The  control  volume  used  in  the  derivation  of  Equation  4. 1 
assumes  that  the  contribution  of  the  pressure  forces  on  the  control  volume  surface  to  be 
negligible.  This  assumption  stems  from  the  work  of  Martinuzzi  and  Tropea  (1993)  who 
show  a  pressure  recovery  at  a  downstream  distance  of  x/d~2.25  for  a  surface  mounted 
cube.  Thus,  due  to  the  roughness  element  geometries  under  investigation,  all  elements 
should  experience  full  pressure  recovery  by  x/d=2.75.  It  is  also  of  significance  to  be 
aware  that  the  analysis  used  to  calculate  the  drag  also  includes  the  stream-wise  Reynolds 

normal  stress,  u 2  ,  term.  Consequently,  this  approach  is  more  rigorous  than  a  method 
which  only  takes  into  account  the  mean  stream-wise  velocity.  The  addition  of  the 
fluctuating  term  will  cause  a  slight  decrease  in  the  calculated  drag  as  compared  to  the 
drag  calculated  using  only  the  mean  velocity.  This  is  due  to  the  fact  that  some  mean  flow 
momentum  is  lost  when  turbulence  is  created.  As  a  result,  the  total  momentum  loss  is 
due  to  the  drag  on  the  element  as  well  as  the  production  of  turbulence  caused  by  the 
element’s  presence  in  the  boundary  layer. 


In  order  to  obtain  the  drag  due  only  to  the  element’s  presence  within  the  boundary  layer 
one  must  take  into  account  the  drag  due  to  the  floor  surface  area  of  the  control  volume. 
To  this  end,  it  is  necessary  to  subtract  off  the  drag  due  to  the  floor  without  the  element 
being  present.  There  are  a  few  critical  factors  to  keep  in  mind  when  performing  this  part 
of  the  analysis.  First,  it  is  imperative  that  knowledge  of  the  wake  decay  in  the  spanwise 
direction  be  known.  The  edge  of  the  control  volume  is  determined  to  be  the  location 
where  the  wake  region  has  decayed  to  the  two-dimensional  smooth  wall  conditions.  If 
this  location  is  not  considered,  it  is  quite  possible  that  too  large  or  too  small  a  quantity 
will  be  subtracted  from  the  incremental  drag  calculation.  The  other  important  aspect  is 
that  the  separated  flow  in  the  wake  region  directly  behind  the  roughness  element  will 
create  a  drastically  different  value  of  wall  shear  as  compared  to  the  regions  farther  away 
from  the  centerline  of  the  element  and  the  reference  smooth  wall.  In  order  to  take  this 
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phenomenon  into  account,  the  area  of  the  wake  region  directly  behind  the  roughness 
element  is  approximated  by  a  triangular  section  extending  from  the  base  of  the  roughness 
element  to  the  x/d=2.75  measurement  plane.  Due  to  the  fact  that  wall  shear  data  within 
this  region  is  extremely  difficult  to  obtain,  a  reasonable  approach  would  be  to  assume  a 
linear  trend  in  the  wall  shear  value  from  the  downstream  edge  of  the  roughness  element 
to  the  x/d=2.75  measurement  plane.  An  average  wall  shear  value  within  this  range  would 
then  be  used  to  find  the  drag  due  to  the  floor  in  this  triangular  section  with  the  roughness 
element  being  present.  This  area  is  in  turn  subtracted  from  the  entire  area  of  the  control 
volume  so  as  not  to  be  accounted  for  twice.  The  drag  on  each  separate  area  of  the  floor  is 
then  calculated  and  subtracted  from  the  incremental  drag  result.  Equation  4.2  below  was 
used  to  obtain  the  drag  due  only  to  the  element,  ADe  actuai, 


actual  incremental 


_  outer 


\pv< 

1  J 


A  +C 

1  lproj  _  outer  _  inner 


proj  __  inner 


(4.2) 


where  Cf_outer  and  Cfjnner  are  the  skin  friction  coefficients  in  the  far  spanwise  wake  and 
near  element  centerline  wake  regions  within  the  control  volume,  respectively.  Here 
Cf  outer  is  taken  to  be  the  reference  smooth  wall  value  and  Cfjnner  is  found  using  the 
averaging  method  previously  discussed.  Similarly,  Apr0j_0uter  and  Aproj  inner  are  the 
projected  floor  surface  areas  of  the  respective  wake  regions.  Table  4.1  below  gives  the 
contributions  of  the  various  terms  used  to  calculate  the  actual  drag  on  each  roughness 
element  as  well  as  the  q  and  Cwwa  terms  in  the  model  equation  used  to  correlate  the 
experimental  drag  data.  A  perturbation  analysis  was  employed  to  calculate  the 
uncertainty  associated  with  ADEjncrementai-  The  uncertainty  in  ADEjncrementai  was  found  to 
be  ±6%.  Discussion  related  to  this  uncertainty  analysis  can  be  found  in  Appendix  C. 
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Table  4.1:  Contributions  of  terms  to  related  drag  equations 


Element 

ADe  incremental 

~(N) 

A 

^proj  inner 

(m2) 

Aproj  outer 

(m2) 

Term  A 
(Eqn.  4.2) 

Term  B 
(Eqn.  4.2) 

q 

CwWA 

Cone 

2.66E-04 

1.69E-05 

8.76E-05 

8.74E-05 

9.12E-06 

124.19 

0.60 

Gaussian 

1.96E-04 

1.51E-05 

6.63E-05 

6.62E-05 

8.16E-06 

0.60 

II5BI 

3.32E-04 

2.47E-05 

1.21E-04 

1 .20E-04 

1.34E-05 

0.55 

Cube  at  90° 

3.02E-04 

1.23E-05 

4.59E-05 

4.58E-05 

3.85E-06 

128.61 

N/A 

Cube  at  45° 

3.13E-04 

2.46E-05 

9.17E-05 

9.16E-05 

7.79E-06 

128.61 

N/A 

2.02E-04 

2.68E-05 

7.16E-05 

7.15E-05 

1 .47E-05 

1 14.35 

0.50 

In  order  to  obtain  a  more  complete  understanding  of  the  effects  that  each  roughness 
element  presents  to  a  turbulent  boundary  layer  flow,  efforts  have  been  made  in  devising  a 
scheme  that  will  correlate  the  calculated  drag  data  from  all  of  the  elements.  This  scheme 
takes  into  account  the  form  (pressure)  drag  related  to  the  frontal  area  of  each  roughness 
element  as  well  as  the  viscous  (friction)  drag  which  corresponds  to  the  windward  wetted 
surface  area  of  each  element.  This  scheme  is  not  designed  to  be  an  a  priori  drag 
calculation  method  or  model  due  the  fact  that,  as  discussed  earlier,  knowledge  of  the 
wake  region  is  necessary  in  order  to  obtain  the  actual  drag  due  to  each  element.  Equation 
4.3  presents  the  model  equation  utilized  to  correlate  the  drag  results  where  C<jP  and  CdV 
are  coefficients  acquired  via  a  least  squares  fit. 

~  calculated  ^  dp  qAf  +C dvCWWA  ju(y)E{y)fdy  (4.3) 

'  I  '  _ _ _ . 

2 

Term  1  takes  into  account  the  pressure  form  drag  on  the  element  where  q  is  the 

stagnation  pressure  averaged  over  the  frontal  area,  Af,  of  the  roughness  element.  The  q 
term  is  taken  from  previous  studies  done  by  Wieghardt  (1946)  and  is  defined  as, 

9  =  2  +  k'  <4-4> 

Af  0  1 

where,  ((/ 2  +  u 2 ]  is  the  approach  streamwise  two-dimensional  mean-square  velocity 
profile,  and  E(y)  is  the  y-variation  of  the  frontal  projected  width  of  the  roughness 
element.  The  E(y)  function  is  known  for  each  element  and  is  integrated  over  the 
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roughness  height,  k,  in  order  to  account  for  the  projected  frontal  area  of  the  roughness 
element. 


Term  2  in  Equation  4.3  considers  the  viscous  contribution  to  the  total  drag.  A  viscous 
stress  is  taken  to  be  of  the  form  /uAU/An  where  An  in  this  instance  is  an  incremental 
distance.  Furthermore,  if  the  smallest  scale  spanwise  structures  for  a  smooth  wall  are 
taken  into  account,  the  following  relationship  can  be  obtained 
AU/An  =  U(y)/(l2v/Ur).  Rearrangement  of  this  relation  leads  to  the  following, 


/jAU/An  -  pU(y)U T  /\2 ,  which  is  what  is  present  in  Term  2  of  Equation  4.3.  The 
E\y)fdy  factor  in  Term  2  is  the  variation  in  the  y-direction  of  the  wetted  area.  In  order 
to  obtain  the  wetted  surface  area  of  each  axisymmetric  element  it  is  necessary  to  integrate 
the  variation  of  the  incremental  circumference  over  the  roughness  height  (k).  Term  2  for 

pUtnk 


the  axisymmetric  elements  would  then  be  CdvC] 


dv^WWA 


12 


\U{y)E{y)fdy ,  where  Cwwa 


0 


Term  2  axisymmetric 


is  a  coefficient  that  is  associated  with  the  fraction  of  total  surface  area  that  sees  attached 
flow,  i.e.  the  windward  wetted  area.  Here,  /  is  a  dimensionless  function  particular  to 
each  roughness  element’s  incremental  variation  in  the  length  of  its  surface  as  a  function 

of  y.  For  example,  the  cone  would  have  a  function  equal  to  fcone  =  ^1  +  ( Sjuo1 ) 

in  order  to  obtain  the  incremental  variation  in  its  surface  area.  To  determine  the  values  of 
Cwwa,  it  is  necessary  to  know  the  separation  locations  on  the  surface  of  the  axisymmetric 
elements.  The  extremely  small  size  of  the  roughness  elements  under  scrutiny  make  it 
very  hard,  if  not  impossible,  to  obtain  these  exact  separation  locations.  Consequently,  it 
is  necessary  to  glean  information  from  previous  studies. 


No  previous  study  has  been  performed  on  roughness  elements  having  comparative  size 
and  shape  characteristics  to  the  ones  under  current  investigation.  An  investigation  was 
done  by  Okamoto  et  al.  (1977)  who  performed  experiments  on  surface  mounted  cones 
with  varying  vertex  angles  (a).  These  experiments  investigate  surface  pressure 
distributions,  wake  properties  and  separation  lines  associated  with  each  cone. 
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Experiments  were  carried  out  on  a  similarly  shaped  cone  as  in  the  current  research.  The 
Okamoto  investigation  examined  a  cone  with  a  vertex  angle  2a=60°  whereas  the  current 
research  considers  a  cone  with  2a=53.1°.  The  60°  vertex  angle  cone  had  a  height  of 
86.43  mm  and  was  placed  in  a  7.6  mm  thick  turbulent  boundary  layer.  Although  the 
parameters  in  the  Okamoto  investigation  are  not  the  same  as  the  current  setup,  which  has 
a  2.54  mm  tall  cone  placed  in  a  39.3  mm  thick  turbulent  boundary  layer,  useful 
comparisons  can  still  be  gathered.  From  the  data  gather  by  Okamoto  et  al.,  it  was  found 
that  the  flow  stayed  attached  to  approximately  45%  of  the  total  surface  area  of  the  60° 
vertex  angle  cone.  It  is  reasonable  to  assume  that  a  roughness  element  of  the  present  size 
(height  <  7%  of  the  boundary  layer  thickness),  fully  immersed  in  a  turbulent  boundary 
would  experience  a  higher  percentage  of  attached  flow  as  compared  to  the  Okamoto 
study.  This  assumption  stems  from  the  fact  that  a  turbulent  boundary  layer  will  have 
more  ‘energy’  and  thus  be  able  to  stay  attached  to  the  surface  longer  as  compared  to  the 
uniform  flow  seen  in  the  ffeestream.  After  thoroughly  scrutinizing  the  close-up  oil  flows 
in  Figure  3.6,  a  value  of  approximately  55-60%  of  the  total  surface  area  was  selected  to 
be  the  attached  flow  surface  area  for  the  peaked  elements.  Consequently,  Cwwa  is  taken 
to  be  0.55-0.6  for  the  peaked  elements. 

A  series  of  experiments  performed  by  Gaudet  (1987)  on  rivet  heads  immersed  in  a 
turbulent  boundary  layer  provide  some  information  on  the  separation  location  on  the 
hemisphere.  As  in  the  Okamoto  case  above  for  the  cones,  the  work  done  by  Gaudet  was 
on  elements  that  were  much  larger  than  the  boundary  layer  thickness  and  the  elements 
tested  were  more  ‘flat’  having  a  height  to  diameter  ratio  of  1:5.  Gaudet,  however,  did 
show  a  trend  in  his  data  that  leads  to  ~50%  of  the  total  surface  area  having  attached  flow 
in  similar  flow  conditions  to  the  current  research.  The  oil  flow  pattern  of  the  near  wall 
flow  around  the  hemisphere,  Figure  3.5,  also  shows  the  flow  on  the  wall  staying  attached 
only  to  the  front  half  of  the  hemisphere.  The  combination  of  these  results  leads  one  to 
use  Cwwa=  0.5  for  the  drag  analysis  done  on  the  hemisphere  in  Equation  4.3. 

When  performing  the  drag  analysis  using  Equation  4.3,  the  cube  elements  are  much 
simpler  due  to  the  fact  that  the  width  of  each  element  is  constant.  Therefore  there  is  no 
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need  for  an  /  function  as  described  above.  The  wetted  surface  area  of  the  cube  at  90°  is 
taken  to  be  the  surface  area  of  three  sides  and  the  top  whereas  the  wetted  surface  area  of 
the  cube  at  45°  is  taken  as  the  surface  area  of  two  sides  and  the  top.  The  surfaces  on  the 
downstream  side,  or  sides,  of  the  elements  will  not  contribute  to  the  viscous  drag  and 
were  deemed  negligible  in  this  analysis. 

Following  the  examination  and  application  of  the  previously  discussed  issues,  a  least 
squares  fit  of  the  ADe  actual  data  to  Equation  4.3  was  employed  to  correlate  the 
incremental  drag  of  the  six  elements  under  scrutiny.  The  two  coefficients  CdP  and  Cdv, 
were  obtained  in  this  fashion  and  found  to  be  0.1  and  14.25,  respectively.  Figure  4.1 
shows  the  results  of  the  incremental  drag  analysis.  The  solid  symbols  are  the  calculated 
values  from  the  least  squares  fit,  whereas  the  open  symbols  are  the  actual  experimental 
values  calculated  from  Equation  4.2.  The  open  symbols  are  plotted  with  respect  to 
themselves  in  order  to  give  an  idea  of  where  the  ideal  calculated  values  should  fall.  The 
calculated  drag  value  for  each  element,  from  Equation  4.3,  is  in  turn  plotted  vs.  its 
respective  experimental  data  point.  The  dashed  line  is  a  linear  fit  to  these  calculated  data 
points  and  has  an  R2  correlation  coefficient  value  equal  to  0.96.  A  better  way  to  analyze 
how  well  the  calculated  values  compare  to  the  experimental  values  is  to  look  at  the 
percent  difference  between  the  two  values.  The  uncertainty  bars  presented  on  the  open 
experimental  symbols  are  ±10%  of  the  respective  experimental  value.  All  calculated 
drag  values  are  within  this  ±10%  range  except  for  the  hemisphere  which  is  ~1 1%  higher 
than  its  respective  experimental  value.  Table  4.2  presents  the  results  for  the  experimental 
(actual)  and  the  calculated  drag  on  each  roughness  element  from  Equation  4.2  and 
Equation  4.3,  respectively.  The  percent  difference  between  the  actual  and  calculated  drag 
value  associated  with  each  roughness  element  is  also  given.  The  discrepancies  between 
the  calculated  and  experimental  values  can  primarily  be  attributed  to  the  uncertainties 
associated  with  the  actual  separation  locations  on  the  surface  of  the  elements  as  well  as 
the  area  of  the  inner  wake  region  and  its  associated  value  of  wall  shear. 
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Table  4.2:  Comparison  of  actual  and  calculated  drag  values 


Element 

ADe  actual  (N) 

ADg  calculated  (N) 

%  Difference 

Cone 

1.69E-04 

1.55E-04 

8% 

Gaussian 

1.22E-04 

1.17E-04 

4% 

1.98E-04 

1.84E-04 

7% 

Cube  at  90° 

2.52E-04 

2.39E-04 

5% 

Cube  at  45° 

2.13E-04 

2.05E-04 

4% 

msmm 

1.15E-04 

1.28E-04 

11% 

The  results  obtained  from  Equation  4.3  help  to  bring  a  certain  sense  of  unity  to  the  drag 
values  obtained  from  the  experiments.  The  largest  value  of  drag  is  seen  on  the  cube  at 
90°  followed  by  the  cube  at  45°,  large  grooved  cone  element,  cone,  Gaussian  element, 
and  the  hemisphere.  The  cube  at  90°  presents  more  of  a  blockage  effect  to  the  flow  and 
is  the  least  ‘streamlined’  element  being  tested.  Thus,  it  is  not  surprising  that  it  has  the 
highest  value  of  drag.  The  cube  at  45°  does  have  a  larger  amount  of  frontal  area  than  the 

cube  at  90°  by  a  factor  of  y/2  ,  but  its  orientation  is  more  aerodynamic  than  that  of  the 
cube  at  90°.  The  larger  value  of  drag  on  the  large  grooved  element,  as  compared  to  the 
other  peaked  elements,  can  be  related  not  only  to  its  frontal  projected  area  but  also  its 
wetted  surface  area.  With  the  highest  wetted  surface  area,  next  to  the  fine  grooved 
element,  the  large  grooved  element  will  experience  more  viscous  drag  than  the  other 
elements  under  scrutiny.  The  rest  of  the  elements  follow  an  apparent  trend  in  that  the 
cone  has  the  next  largest  amount  of  frontal  and  wetted  area  of  the  peaked  elements 
followed  by  the  Gaussian  element.  The  hemisphere  has  the  smoothest  surface  of  all  the 
elements  and  is  also  submerged  the  lowest  in  the  boundary  layer.  Consequently,  the 
lower  speed  fluid  moving  over  the  hemisphere  in  conjunction  with  the  hemisphere’s 
smooth  shape  will  generate  the  least  amount  of  drag  of  all  elements. 


Other  efforts  have  been  made  to  find  possible  drag  correlations  using  just  the  frontal  and 
wetted  areas  of  the  elements.  The  results  of  these  efforts  can  be  seen  in  Figures  4.2  and 
4.3.  The  drag  coefficient  Cd  is  denoted  by  a  (frontal)  or  a  (wetted)  subscript  indicating 
that  either  the  frontal  or  wetted  area  of  the  element  was  used  in  the  normalization  process. 
The  drag  coefficient  is  obtained  using  the  following  equation. 
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C 


d  ( frontal  _  or  _  wetted ) 


_  actual 

\pU'.A,_„ 


(4.5) 


Similarly,  results  are  also  plotted  in  Figure  4.2  using  a  drag  coefficient  defined  by 
Equation  4.6.  This  drag  coefficient  was  first  used  by  Wieghardt  (1946). 


C 


d  _  Wieghardt 


E  _  actual 
frontal 


(4.6) 


Each  drag  coefficient  is  presented  as  a  ratio  with  respect  to  the  reference  smooth  wall 
skin  friction  coefficient,  Cf.  These  drag  ratios  are  plotted  vs.  the  normalized  circulation 
TfaT(Af  or  Sf2)  which  utilizes  the  circulation,  f,  from  the  x/d=2.75  plane  and  the 


reference  smooth  wall  value  of  UT. 


The  relationship  between  the  various  drag  data  normalized  with  respect  to  the  element 
frontal  area  show  a  fairly  flat  trend  as  Cd(fromai/Cf  is  somewhat  constant  with  increasing 
r/(lfT  ( Af  or  wJ/2 ).  The  cube  at  90°  tends  to  be  the  outlier  in  view  of  the  fact  that  it  has 

the  highest  amount  of  drag  with  the  next  to  lowest  (Gaussian  has  least)  frontal  area.  This 
result  is  not  extremely  surprising  due  to  the  fact  that  the  drag  of  the  roughness  element  is 
related  to  both  the  viscous  drag  and  the  form  drag  as  was  discussed  earlier.  With  this  in 
mind,  looking  at  Figure  4.3  reveals  that  the  cube  at  90°  is  still  separated  from  the  other 
elements  even  though  the  wetted  area  of  the  cube  is  the  third  largest  behind  both  of  the 
grooved  elements.  The  blockage  created  by  the  cube  at  90°  decreases  the  amount  of 
circulation  produced,  thus  offsetting  it  from  the  other  values  once  more.  A  somewhat 
linear  trend  can  be  seen  in  the  data  of  Figure  4.3  as  the  normalized  drag  increases  so  does 
the  normalized  circulation.  After  studying  Figures  4.1  through  4.3,  it  is  obvious  that  in 
order  to  grasp  the  physics  of  what  is  happening  related  to  the  drag  of  each  element  it  is 
essential  that  the  analysis  encompass  both  the  form  drag  of  the  element  as  well  as  the 
viscous  drag. 


Chapter  4:  Discussion  of  Results 


-  193- 


4.3  Downstream  Flow  Field  Comparisons 

In  an  attempt  to  bring  the  immense  amount  of  data  gathered  in  this  study  into  a  more 
condensed  and  useful  form,  drawings  for  the  flow  around  each  roughness  element  were 
completed  and  can  be  seen  in  Figures  4.4  through  4.10.  For  conciseness,  a  legend 
common  to  all  drawings  can  be  found  in  Figure  4.4.  Each  drawing  shows  the  measured 
features  at  the  x/d=2.75  plane  for  all  elements  except  the  fine  grooved  element  where  the 
plane  at  x/d=1.36  is  shown.  Regions  of  vorticity  and  turbulent  kinetic  energy  (TKE) 
maxima  are  shown  at  these  measurement  planes.  The  TKE  diffusion  velocity  vectors, 

VjUr=\^/7]/ur  and  WjUt^/7]/Ut  are  also  presented  in  the  region 

nearest  the  maximum  TKE.  These  vectors  are  only  presented  in  the  peak  TKE  region  for 
each  element  in  order  for  one  to  grasp  the  magnitude  and  direction  of  the  greatest 
diffusion  of  TKE. 

As  is  evidenced  by  the  drawings,  all  elements  produce  a  streamwise  horseshoe  vortex 
structure  that  is  on  the  order  of  the  size  of  the  roughness  element  from  which  it  emanates. 
At  the  x/d=2.75  plane,  the  secondary  velocity  produced  by  the  horseshoe  vortex 
structures  is  very  significant.  The  magnitude  of  the  secondary  velocity  at  this  location 
causes  a  large  amount  of  fluid  to  be  brought  toward  the  wall.  All  elements  at  this 
location  show  a  significant  normal-to-wall  (V)  velocity  at  the  centerline  (Figure  3.12). 
The  highest  magnitude  of  negative  V  can  be  seen  in  the  two  cube  elements,  followed  by 
the  hemisphere,  cone,  Gaussian,  and  the  large  grooved  element.  From  Figure  3.12,  the 
fine  grooved  element’s  V  at  the  x/d=2.75  plane  is  very  near  the  magnitude  of  the 
Gaussian  element.  Thus  it  is  reasonable  to  assume  that  the  horseshoe  vortex  structure 
produced  by  the  fine  grooved  element  at  x/d=2.75  is  close  to  the  same  strength  as  the 
structure  produced  by  the  Gaussian  element.  All  counter  rotating  horseshoe  vortices 
move  away  from  the  wall  and  centerline  and  dissipate  as  they  propagate  downstream. 
The  movement  and  strength  of  the  vortex  structure  will  be  affected  a  great  deal  by  the 
generation  of  secondary  vorticity  near  the  wall.  In  this  region  the  near  wall  viscous 
forces  will  cause  the  loss  of  strength  in  the  primary  vortex. 
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Figure  4.7  for  the  fine  grooved  element  shows  that  at  the  x/d=1.36  measurement  plane 
there  is  fluid  being  swept  in  around  the  sides  of  the  element.  This  location  precedes  the 
reattachment  point  and  the  horseshoe  vortex  structure  has  not  yet  taken  over  the 
secondary  flow  field.  Therefore,  the  dominant  secondary  flow  structure  appears  to  be  the 
flow  coming  around  the  sides  of  the  element  rather  than  the  counter-rotating  structure 
shown  at  the  x/d=2.75  plane.  All  peaked  elements,  Figures  4.4  to  4.7,  are  thought  to  have 
weak  vortices  emanating  from  the  top  of  the  element.  Although  the  presence  of  these 
vortices  could  not  be  determined  from  the  experimental  data,  they  have  been  shown  to  be 
present  in  a  study  by  Fontaine  and  Deutsch  (1996)  on  a  smaller  Gaussian  shaped 
element.  The  strength  of  these  smaller  ‘peak’  vortices  is  orders  of  magnitude  less  than 
the  primary  horseshoe  vortex  structure. 

Regions  of  positive  mean  streamwise  vorticity,  Qx,  can  be  located  in  all  drawings.  The 
magnitude  of  maximum  streamwise  vorticity  is  given  in  each  drawing  as  well  as  its 
associated  location.  The  location  of  the  vortex  center  is  taken  to  be  the  point  of 
maximum  vorticity  as  was  also  done  in  previous  studies  by  Wendt  and  Greber  (1992). 
All  spanwise  (z)  position  values  are  normalized  by  the  parameter  d  for  each  element. 
This  parameter  is  simply  the  element’s  base  diameter  for  all  axisymmetric  elements,  the 
edge  length  for  the  cube  at  90°,  and  the  length  of  the  diagonal  for  the  cube  at  45°.  The 
vertical  (y)  position  values  are  normalized  by  the  height,  h,  of  each  element. 

The  location  of  maximum  streamwise  vorticity  for  all  elements  is  quite  close  to  the  wall 
at  x/d=2.75.  The  cube  at  45°  and  the  hemisphere  show  the  largest  normalized  y/k 
locations  from  the  wall  at  y/k=0.138  and  0.183  respectively  where  the  rest  of  the 
elements  are  in  the  range  of  y/k^. 83-0.96.  It  is  worth  noting  that  the  distance  of 
maximum  streamwise  vorticity  for  the  fine  grooved  element  at  x/d=1.36  is  very  near  the 
wall  at  y/k=0.039.  These  results  show  the  significance  of  having  near  wall  data  due  to 
the  large  amount  of  activity  that  is  seen  in  this  locale.  The  major  source  of  vorticity  at 
the  x/d=2.75  location  is  below  y+=15  for  all  elements.  The  spanwise  distance  from  the 
wall  for  most  elements  is  similar  and  around  z/d~-0.83.  This  is  not  the  case  for  the 
Gaussian  element  and  the  hemisphere  which  have  maximum  vorticity  locations  in  the 
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spanwise  direction  equal  to  z/d=  -0.552  and  -0.551  respectively.  This  result  is  expected 
due  the  shape  of  these  elements.  The  Gaussian  element  is  fairly  slender  along  its  height, 
thus  causing  the  vortex  structures  to  be  closer  to  the  centerline.  The  hemisphere  on  the 
other  hand  is  very  smooth  and  has  more  flow  coming  over  its  top  and  rushing  toward  the 
wall;  thus  bringing  the  horseshoe  vortex  structures  closer  to  the  centerline  in  the  near 
element  regions. 

The  magnitude  of  the  vorticity  is  directly  related  to  the  strength  of  the  horseshoe  vortex 
structure.  As  a  result,  the  maximum  value  of  normalized  streamwise  vorticity, 

{flxyj~A^ ju r)  ,  at  x/d=2.75  is  observed  in  the  cube  elements.  The  cube  at  90°  shows 

the  largest  magnitude  equal  to  1 1 .59  and  the  cube  at  45°  is  slightly  less  at  1 1 .35.  The  rest 
of  the  elements  tend  to  be  fairly  close  to  one  another  with  the  large  grooved  element 
having  the  smallest  normalized  vorticity  maximum  at  4.44.  The  presence  of  the  larger 
grooves  (2.5  times  the  smallest  sublayer  structure  length  scale)  seems  to  have  a 
noticeable  effect  on  the  flow  field  around  the  roughness  element.  The  irregularities  that 
are  2.5  times  the  smallest  sublayer  structure  length  scale  not  only  increase  the  viscous 
drag  on  the  element,  but  they  also  tend  to  ‘straighten’  the  flow,  produce  weaker  vortices 
and  less  vorticity  than  the  cone  element  with  equal  frontal  area  and  similar  shape.  The 
finer  grooves  that  are  the  size  of  the  smallest  sublayer  structure  length  scale  also  increase 
the  viscous  drag  on  the  element,  but  do  not  have  as  much  effect  on  the  production  of 
vorticity  as  the  larger  grooves.  Due  to  the  fact  that  the  vortices  produced  by  the  elements 
are  located  very  close  to  the  wall,  as  was  discussed  above,  viscous  forces  will  have  a 
major  contribution  in  the  decay  of  the  vortices  as  they  propagate  downstream.  Another 
aspect  that  is  related  to  the  diffusion  and  transport  of  the  vortices  is  the  entrainment  of 
fluid  from  the  boundary  layer. 

The  zone  of  separated  flow  for  each  element  is  shown  in  Figures  4.4  through  4.10  to  be 
directly  behind  each  roughness  element  and  of  varying  breadth  for  the  different  cases. 
The  defined  separation  points  on  the  sharp  edges  of  the  cube  elements  create  a  larger 
region  of  separated  flow  as  compared  to  the  axisymmetric  elements.  The  cube  elements 
also  present  a  significant  ‘blockage’  effect  to  the  flow.  This  blockage  creates  a  region  of 
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stagnant  fluid  and  tends  to  delay  flow  reattachment.  The  presence  of  an  arch  vortex  can 
be  seen  in  the  flow  field  directly  downstream  of  the  cube  elements.  This  vortex  is  formed 
from  the  separation  and  shear  layer  roll-up  created  near  the  sharp  edges  of  the  element. 
The  horseshoe  vortex  structure  tends  to  bring  fluid  in  toward  the  element  along  the 
centerline.  Smaller  vortices  can  also  be  found  on  the  sides  and  top  of  the  cube  at  90°  due 
to  the  shear  layers  rolling  up  at  these  edges.  The  separated  flow  regions  behind  the 
axisymmetric  elements  are  smaller  and  not  as  well  defined  as  those  associated  with  the 
cube  elements.  The  size  of  the  separated  flow  region  obviously  depends  on  the  size  and 
shape  of  the  element.  The  more  slender  shape  of  the  Gaussian  element  lends  itself  to  the 
creation  of  the  smallest  region  of  separated  flow.  Thus,  the  Gaussian  element  is  the  only 
element  not  to  show  backflow  at  the  x/d=1.36  measurement  location.  The  large  grooved 
element  has  the  second  to  smallest  portion  of  separated  flow  due  to  the  fact  that  the  base 
of  the  element  does  not  extend  all  the  way  to  the  wall.  The  first  groove  machined  into  the 
element  starts  at  the  base  of  the  element,  thus  creating  an  opening  for  the  flow  to  proceed 
through.  The  hemisphere,  cone  and  fine  grooved  elements  have  comparable  separated 
regions  but  the  hemisphere  reattaches  sooner  due  to  the  flow  accelerating  over  the  top  of 
the  element  and  being  directed  toward  the  wall.  The  near  wall  separation  point  for  all 
peaked  elements  tends  to  be  on  the  downstream  side  of  the  element.  Accordingly,  it  is  a 
reasonable  assumption  due  to  the  size  of  the  element  and  associated  boundary  layer 
parameters  that  the  flow  stays  attached  to  approximately  55-60%  of  the  element’s 
surface.  The  hemisphere  element  on  the  other  hand  seems  to  have  a  separation  point 
nearer  to  its  midpoint  owing  to  the  reasoning  that  approximately  50%  of  the  surface  sees 
attached  flow.  As  was  discussed  in  Section  4.2,  these  conclusions  stem  from  the  oil  flow 
visualizations  done  on  the  roughness  elements  as  well  as  the  work  done  by  Okamato  et 
al.  (1977)  and  Gaudet  (1987). 

The  magnitude  and  location  of  the  peak  TKE  at  each  element’s  respective  measurement 
plane  is  presented  in  the  flow  field  drawings  as  well.  Although  all  maximum  values 
occur  along  the  centerline  of  the  element,  the  normalized  height  ratios  vary  significantly 
between  the  axisymmetric  elements  and  the  cube  elements.  The  large  regions  of  TKE  are 
due  to  the  horseshoe  vortices  that  not  only  create  sweeping  motions  toward  the  wall  but 
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also  a  convergence  of  associated  shear  layers  toward  the  center  plane.  The  largest 
magnitude  of  normalized  TKE  at  the  x/d=2.75  measurement  plane  is  produced  by  the 
cube  at  45°  with  a  value  equal  to  21.29.  The  cube  at  90°  follows  with  a  peak  value  of 
18.36.  The  fine  grooved  element  also  has  a  rather  large  peak  value  (17.59),  as  shown  in 
Figure  3.103  at  the  x/d=2.75  location.  The  cone,  large  grooved,  Gaussian,  and 
hemisphere  elements  follow  with  peak  TKE  values  of  15.31,  13.96,  13.82,  and  13.42, 
respectively.  The  larger  elements  that  protrude  higher  into  the  boundary  layer  and  have 
sharp  edges  tend  to  produce  more  TKE.  The  hemisphere  which  is  shorter  and  has  smooth 
sides  has  the  lowest  amount  of  TKE  production.  The  fine  grooved  element  shows  a  peak 
value  at  the  x/d=1.36  plane  equal  to  33.15  which  is  nearly  twice  as  much  as  the  peak 
value  at  the  x/d=2.75  centerline  measurement  location. 

The  diffusion  velocity  vectors  are  shown  in  the  immediate  area  surrounding  the  peak 
value  of  TKE  for  all  elements.  The  highest  diffusion  velocities  are  produced  by  the  cone 
and  cubes  whose  vectors  have  very  similar  magnitudes.  The  Gaussian,  fine  grooved, 
large  grooved,  and  hemisphere  elements  follow  with  decreasing  diffusion  velocities.  The 
fine  grooved  element  is  more  ‘peaky’  than  the  large  grooved  element  and  has  greater 
TKE  diffusion.  As  was  discussed  in  Chapter  3,  the  diffusion  velocities  are  highest  near 
sharp  edges  (cubes)  and  near  sharp  peaks  (cone,  Gaussian,  etc.).  Subsequently,  this  leads 
to  the  conclusion  that  both  the  element  shape  and  the  orientation  play  very  significant 
roles  in  the  production  and  diffusion  of  TKE.  Although  the  hemisphere  has  a  smoother 
surface  than  the  other  elements,  there  are  still  significant  diffusion  velocities  near  the 
element  height.  The  height  at  which  the  greatest  magnitude  of  TKE  exists  at  the 
x/d=2.75  location  is  very  similar  for  the  peaked  elements  as  they  all  are  within  the  range 
0.550<y/k<0.555.  This  result  can  be  attributed  to  the  fact  that  the  horseshoe  vortex 
structures  produced  by  the  peaked  elements  are  nearly  the  same  size  at  this  location. 
Therefore,  the  sweeping  motions  toward  the  wall  as  well  as  the  convergence  of  associated 
shear  layers  toward  the  center  plane  by  the  horseshoe  vortices  will  be  at  similar  positions 
within  the  boundary  layer.  The  normalized  y-location  of  the  maximum  TKE  associated 
with  the  hemisphere  is  y/k=0.555,  thus  making  it  nearly  the  same  as  the  peaked  elements. 
The  normalized  y-location  of  the  maximum  TKE  production  at  the  x/d=1.36  plane  for  the 
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fme  grooved  element  is  =38%  higher  than  its  value  at  x/d=2.75.  The  high  TKE  value 
seen  at  the  x/d=1.36  measurement  location  is  significantly  influenced  by  the  peak  of  the 
element  as  it  causes  intense  mixing  near  the  element.  As  the  flow  progresses 
downstream,  the  influence  of  the  horseshoe  vortices  tends  to  have  more  of  an  effect  on 
the  magnitude  of  the  TKE.  The  cube  elements  have  normalized  height  values  that  differ 
from  the  axisymmetric  elements.  The  cube  at  90°  has  its  maximum  TKE  location  at 
y/k=0.714,  whereas  the  cube  at  45°  sees  its  maximum  value  at  y/k=0.479.  These 
differences  can  be  linked  to  the  size  and  shape  of  the  cube  elements  as  they  obviously 
have  sharp  edges  and  are  more  blunt  than  the  axisymmetric  elements.  These  facets  will 
in  turn  produce  horseshoe  vortex  structures  that  are  significantly  different  in  size,  shape 
and  strength. 


4.4  Circulation 
4.4.1  Circulation  Decay 

The  decay  of  circulation  is  mainly  due  to  near  wall  viscous  forces  in  the  spanwise 
direction  that  act  in  retarding  the  secondary  flow  (Wendt  and  Greber,  1992).  Viscous 
diffusion  also  plays  a  role  in  the  dissipation  of  the  circulation  as  the  vortex  structure  is 
located  in  the  near-wall  region  of  the  turbulent  boundary  layer.  The  streamwise  decay  of 

circulation,  r^A^),  plotted  using  a  logarithmic  scale  can  be  seen  in  Figure  4.11. 

The  peaked  elements  show  a  common  trend  in  the  decay  of  circulation  with  streamwise 
distance  x/d.  By  fitting  a  curve  to  the  data  it  is  possible  to  obtain  the  general  behavior  of 
the  decay  of  circulation  with  streamwise  distance.  Consequently,  the  peaked  elements 

have  a  streamwise  decay  in  circulation  that  is  proportional  to  *  (x/d)~1'2 .  The  exponent 
-1.12  was  obtained  by  taking  an  average  value  from  all  peaked  elements.  All  actual 
values  obtained  from  the  individual  decay  of  each  element’s  circulation  are  all  relatively 
close  to  the  stated  average  value.  The  cone  is  the  farthest  from  the  average  value  with  a 
difference  of  1 1  %.  The  fine  grooved  element,  Gaussian  element,  and  large  grooved 
element  follow  with  respective  %  differences  of  7%,  3%,  and  0.4%.  These  rather  small 
percent  differences  lead  to  a  confidence  in  the  general  trend  related  to  the  streamwise 
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decay  of  circulation.  Knowing  a  general  trend  in  the  decay  circulation  is  of  great  benefit 
in  practical  applications.  Such  knowledge  can  give  a  qualitative  idea  of  the  behavior  of 
the  downstream  secondary  flow  field  and  strength  of  the  horseshoe  vortex  structure.  The 
hemisphere  and  cube  elements  do  not  have  a  discernible  trend  that  can  be  related  to  one 
another.  Starting  at  the  x/d=5.0  measurement  location  and  continuing  downstream,  the 
hemisphere  tends  to  have  lower  values  of  circulation  as  compared  to  the  other  elements. 
This  in  turn  will  lead  to  the  faster  decay  of  circulation  due  to  the  near  wall  viscous  effects 
taking  over  as  the  vortex  structure  propagates  downstream.  On  the  other  hand,  the  cube 
elements  have  a  totally  different  near  element  behavior  as  was  discussed  in  Chapter  3. 
Once  the  x/d=2.75  measurement  location  is  reached  though,  the  cube  elements  show  a 
decay  more  along  the  lines  of  the  other  elements. 

4.4.2  Streamwise  and  Spanwise  Circulation  Relationship 

There  is  a  significant  flux  in  vorticity  between  the  upstream  and  downstream  regions  of 
the  roughness  element  in  the  turbulent  boundary  layer.  This  change  can  be  attributed  to 
the  pressure  gradient  induced  by  the  roughness  element  which  in  turn  generates  the 
additional  vorticity  seen  in  the  wake  of  the  element.  The  vorticity  in  the  x-y  plane  along 
the  centerline  of  the  element  can  be  related  to  the  generation  of  additional  vorticity  in  the 
y-z  plane  at  the  x/d=2.75  measurement  location.  In  order  to  discuss  this  phenomenon,  it 
is  necessary  to  define  two  planar  control  volumes.  The  x-y  planar  control  volume  (CV)  is 
defined  with  an  upstream  location  at  the  front  of  the  roughness  element  and  a 
downstream  location  at  x/d=2.75.  The  height  of  the  x-y  planar  CV  is  at  the  location 
where  the  mean  streamwise,  U,  velocity  is  constant  (i.e.  at  y=8).  The  y-z  planar  CV  is 
located  at  x/d-2.75  and  has  the  same  dimensions  as  the  x-y  CV.  The  y-z  CV  has  a  side 
located  along  the  centerline  (x-axis)  of  the  element.  It  is  possible  to  arrive  at  the 
following  equations  below  if  the  upstream  and  far  spanwise  V  velocity  profiles  are  taken 
to  be  negligible  along  with  the  observable  fact  of  flow  symmetry  along  the  centerline 
(W=0).  The  circulation  in  the  x-y  plane,  Tz,  is  defined  as 

Tz  =  \v{y)dy-  \u{y)dx  (4.7) 

Ay  Ax 

where  Ay  and  Ax  are  taken  to  be  the  edges  of  the  control  volume  circuit  in  the  y  and  x 
directions.  Similarly  the  circulation  in  the  y-z  plane,  Tx,  can  then  be  defined  as 
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T,  —  ]v{y)dy  (4.8) 

where  Ay  is  taken  to  be  the  same  as  above. 

From  Equations  4.7  and  4.8  it  is  clear  that  the  streamwise  and  spanwise  circulation  are 
directly  related  through  the  V  velocity  profile  at  the  x/d=2.75  centerline  measurement 
location.  The  only  difference  between  the  streamwise  and  spanwise  circulation  is  seen  in 
the  mean  streamwise  (U)  velocity  term  that  is  present  in  the  spanwise  circulation.  This 
tells  us  that  a  significant  amount  of  the  Tz  in  this  control  volume  is  due  to  the  strong 
downwash  produced  by  the  horseshoe  vortex  structures  created  by  presence  of  the 
roughness  element.  Knowing  this  relationship  is  beneficial  due  to  the  fact  that  it  provides 
a  deeper  understanding  of  how  Tz  is  related  to  rx.  Similarly,  the  streamwise  and 
spanwise  vorticity  values,  Qx  and  £2Z,  can  be  related.  This  in  turn  allows  one  to  have  an 
understanding  of  the  connection  between  the  additional  vorticity  created  in  the  wake  of 
the  element  as  compared  to  the  incoming  spanwise  vorticity. 


4.5  Future  Work 

Future  plans  that  will  be  done  with  the  knowledge  gained  from  the  current  research  will 
include  the  examination  of  various  fetches  of  roughness.  These  fetches  will  have  varying 
values  of  lambda  (X)  where  A=(projected  frontal  area)/(total  surface  area)  in  order  to 
acquire  a  more  complete  data  base  for  different  roughness  regimes.  Sparsely  distributed 
roughness  (small  X  values)  as  well  as  more  densely  packed  roughness  (large  X  values) 
will  be  examined.  Figure  4.12  presents  two  examples  of  smaller  spacing  between 
roughness  elements  which  in  turn  would  create  larger  values  of  X.  Once  a  fairly  complete 
data  base  is  compiled,  the  ‘effectiveness’  of  a  particular  roughness  element  can  be 
determined.  The  ‘effectiveness’  of  a  roughness  element  or  a  roughness  regime  will 
provide  a  summary  of  how  a  given  roughness  element  geometry  affects  the  boundary 
layer  flow  conditions.  This  ‘effectiveness’  parameter  is  based  on  the  size,  shape,  and 
spacing  of  the  roughness  elements  as  well  as  the  on-coming  flow  characteristics  upstream 
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of  the  elements.  Once  the  ‘effectiveness’  factor  is  known,  it  could  be  superimposed  into 
a  certain  flow  field  in  order  to  more  closely  calculate  the  flow  characteristics. 


Calculated  aDe 
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4.6  Chapter  4  Figures 


Experimental  ADE 


Figure  4.1:  Drag  correlation  results 


Chapter  4:  Discussion  of  Results 


-203- 


With  Normal  Stresses  (x/d=2.75) 


r/(U<(Af)1'2) 

Figure  4.2:  Variation  of  normalized  drag  coefficient  based  on  frontal  area  C<j(frontai)/Cr 
With  Normal  Stresses  (x/d=2.75) 
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Figure  4.3:  Variation  of  normalized  drag  coefficient  based  on  wetted  area  Cd(wetted)/Cf 
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Figure  4.4:  Drawing  of  the  flow  around  the  cone  element 
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Figure  4.5:  Drawing  of  the  flow  around  the  Gaussian  element 
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Figure  4.6:  Drawing  of  the  flow  around  the  large  grooved  element 
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Figure  4.7:  Drawing  of  the  flow  around  the  fine  grooved  element 
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Figure  4.8:  Drawing  of  the  flow  around  the  cube  at  the  90°  orientation 
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Chapter  5  Conclusions 


A  detailed  study  has  been  conducted  in  the  Virginia  Tech  Aerospace  and  Ocean 
Engineering  Department  small  boundary  layer  wind  tunnel  with  the  primary  intent  to 
fully  measure  and  understand  the  effects  of  various  shaped  roughness  elements  in  a  two- 
dimensional  high  Reynolds  number  turbulent  boundary  layer.  The  roughness  elements 
under  scrutiny  were  less  than  7%  of  the  boundary  layer  thickness  in  height  and  are  as 
follows:  cone,  cone  with  spatial  variations  equal  to  the  smallest  sublayer  structure  length 
scale,  cone  with  spatial  variations  equal  to  2.5  times  the  smallest  sublayer  structure  length 
scale,  Gaussian-shaped  element,  hemisphere,  cube  aligned  perpendicular  to  the  flow 
(cube  at  90°),  and  a  cube  rotated  45°  relative  to  the  flow.  The  roughness  element  heights, 
k+,  non-dimensionalized  by  the  friction  velocity  (UT)  of  the  approaching  turbulent 
boundary  layer,  are  145,  145,  145,  145,  80,  98,  and  98,  respectively.  A  brief  analysis  of  a 
three-dimensional  fetch  of  the  same  Gaussian-shaped  elements  described  previously  was 
also  undertaken.  In  order  to  analyze  the  complex  flow  fields,  detailed  measurements 
were  obtained  using  a  fine-measurement- volume  (50  micron  diameter)  three-velocity 
component  laser-Doppler  velocimetry  (LDV)  system.  The  LDV  measurement  system 
was  used  to  obtain  boundary  layer  profiles  along  the  centerline  of  the  element  at  varying 
downstream  distances.  Two  planes  of  data  were  also  acquired  at  different  downstream 
locations  for  each  element.  Surface  oil  flow  patterns  were  taken  for  each  isolated 
roughness  element  in  order  to  give  a  qualitative  idea  of  the  near  wall  flow  structure.  This 
is  the  first  time  that  a  detailed  fundamental  study  has  been  done  on  roughness  elements  of 
these  shapes  and  sizes.  Moreover,  it  is  essential  to  note  that  no  detailed  analysis  exists 
concerning  the  effects  of  spatial  smoothing  related  to  these  types  of  roughness  elements. 

Through  this  research  it  was  found  that  a  horseshoe  vortex  structure  forms  in  front  of  all 
elements  and  is  on  the  order  of  the  size  of  the  roughness  element  from  which  it  emanates. 
Evidence  of  the  upstream  formation  of  these  vortices  can  be  seen  in  the  oil  flow 
visualizations  performed  on  each  isolated  roughness  element.  As  the  horseshoe  vortex 
structure  proceeds  downstream,  its  presence  induces  the  downwash  of  higher  momentum 
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fluid  toward  the  wall.  This  ‘sweep’  motion  not  only  creates  high  Reynolds  stresses 

(v2  ,  w2 ,  -«v)  downstream  of  the  element,  but  also  leads  to  higher  skin-friction  drag. 
Triple  products  of  velocity  fluctuation  components  were  also  found  to  be  very  significant 
near  the  height  of  the  element.  These  parameters  are  important  in  regards  to  the  turbulent 
diffusion  of  the  turbulent  kinetic  energy  in  the  flow.  The  larger  elements  that  protrude 
into  the  boundary  layer  and  have  sharp  edges  tend  to  produce  more  TKE.  The  largest 
magnitude  of  TKE  was  produced  by  the  cube  at  45°  followed  by  the  cube  at  90°,  the  fine 
grooved  element,  cone,  large  grooved,  Gaussian,  and  finally  the  hemisphere  element. 
The  cone  and  cubes  have  the  highest  diffusion  velocities  with  their  respective  vectors 
having  very  similar  magnitudes.  This  indicates  that  sharp  edges  and  peaks  produce 
greater  diffusion  of  TKE.  Likewise,  both  the  shape  of  a  roughness  element  as  well  as  its 
orientation  play  vital  roles  in  the  diffusion  of  TKE. 

The  oil  flow  patterns  indicate  the  size  of  the  separated  flow  regions  directly  behind  each 
element,  as  well  as  a  qualitative  idea  of  the  horseshoe  vortex  behavior  in  the  near  wall 
region.  The  cube  elements  were  determined  to  have  an  arch  vortex  that  formed  directly 
downstream  of  the  element.  The  spanwise  variation  in  wall  shear  was  very  similar  for  all 
of  the  axisymmetric  elements,  whereas  the  cube  elements  showed  similar  trends  as 
compared  to  one  another.  A  distinct  peak  value  in  the  centerline  variation  in  wall  shear 
occurs  for  the  axisymmetric  elements  at  x/d=5.0  as  the  flow  from  around  the  elements  as 
well  as  the  flow  from  over  the  top  converged  on  the  centerline.  The  cube  elements 
produced  drastic  reductions  in  wall  shear  in  the  near  element  measurement  locations  due 
to  the  significant  blockage  created  by  the  bluntness  of  the  cubes.  At  the  x/d=2.75  plane, 
all  elements  had  the  spanwise  variation  in  wall  shear  decrease  directly  on  either  side  of 
the  centerline  and  gradually  climb  to  the  smooth  wall  value.  This  is  due  to  the  horseshoe 
vortex  structure  causing  significant  downwash  in  the  centerline  region  and  then  an 
associated  upwash  region  at  the  off-centerline  locations  as  fluid  is  directed  away  from  the 
wall.  This  phenomenon  was  confirmed  by  the  oil  flow  patterns  as  dark  high  shear 
scoured  zones  near  the  centerline  followed  by  deposits  of  tracer  particles  in  the  off- 
centerline  lower  shear  regions. 
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It  was  discovered  that  the  cube  oriented  45°  relative  to  the  flow  produced  a  stronger 
downstream  counter-rotating  horseshoe  vortex  pair  than  the  cube  oriented  90°  relative  to 
the  flow.  As  was  previously  discussed,  the  same  result  was  also  found  to  be  the  case  in 
the  study  performed  by  Castro  and  Robins  (1977).  In  the  present  research,  the  highest 
magnitude  of  secondary  flow  can  be  seen  in  the  two  cube  elements,  followed  by  the 
hemisphere,  cone,  Gaussian,  and  the  large  grooved  element. 

The  strengths  of  the  time  mean  horseshoe  vortex  structures  are  evaluated  through  the 
investigation  of  the  vorticity  and  circulation  of  each  vortex.  The  major  source  of 
vorticity  at  the  x/d=2.75  location  is  below  y+=15  for  all  elements.  It  was  found  that  the 
maximum  value  of  streamwise  vorticity,  at  the  x/d=2.75  measurement  plane,  is  produced 
by  the  cube  elements.  The  cube  at  90°  shows  the  largest  magnitude  with  the  cube  at  45° 
slightly  less.  The  rest  of  the  elements  tend  to  be  fairly  close  to  one  another  with  the  large 
grooved  element  having  the  smallest  vorticity  maximum.  The  presence  of  the  larger 
grooves  (2.5  times  the  smallest  sublayer  structure  length  scale)  seems  to  have  a 
noticeable  effect  on  the  flow  field  around  the  roughness  element.  The  irregularities  that 
are  2.5  times  the  smallest  sublayer  structure  length  scale  not  only  facilitate  the  increase  in 
viscous  drag  on  the  element,  but  they  also  tend  to  ‘straighten’  the  flow  and  produce 
weaker  vortices  than  the  cone  element  with  equal  frontal  area  and  similar  shape.  The  fine 
grooved  element  (spatial  variations  equal  to  the  smallest  sublayer  structure  length  scale) 
also  increase  the  viscous  drag  on  the  element  but  do  not  seem  to  have  near  the  effect  on 
the  production  of  vorticity  as  the  larger  grooves  have  shown. 

The  upstream  measurement  planes  also  reveal  the  presence  of  negative  vorticity  in  the 
near  wall  region.  The  presence  of  these  negative  regions  originates  from  the  convection 
of  wall  induced  vorticity  caused  by  the  interaction  of  the  horseshoe  vortex  with  the  no¬ 
slip  condition  at  the  wall.  These  negative  regions  of  vorticity  are  smaller  and  of  less 
magnitude  as  compared  to  their  corresponding  positive  vorticity  regions.  These  regions 
tend  to  grow  in  size  and  magnitude  as  the  vortex  structure  propagates  downstream.  The 
movement  and  strength  of  the  vortex  structure  is  affected  a  great  deal  by  the  generation 
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of  secondary  vorticity  near  the  wall.  In  this  region  the  near  wall  viscous  forces  cause  the 
loss  of  strength  in  the  primary  vortex. 

The  ‘peakiness’  of  the  roughness  element  was  found  to  have  a  direct  correlation  to  the 
production  of  circulation,  whereas  the  spatial  smoothing  does  not  have  an  immense  effect 
on  this  parameter.  A  trend  in  the  streamwise  decay  of  circulation,  r/(t/r^/i^),  was 
observed  for  the  peaked  elements.  These  elements  showed  a  general  trend  in  their  decay 
that  is  proportional  to  »  \xjd)12 .  Knowledge  of  a  trend  in  the  decay  of  circulation  can 
give  a  qualitative  idea  of  how  the  downstream  flow  field  is  behaving  regarding  the 
strength  of  the  horseshoe  vortex  structure  and  its  related  secondary  flow  field. 

A  model  equation  has  been  proposed  to  correlate  the  drag  of  the  roughness  elements 
under  scrutiny.  This  equation  is  not  an  a  priori  method  for  calculating  the  drag  due  to 
roughness.  However,  the  development  of  the  equation  has  provided  a  deeper 
understanding  of  the  various  parameters  that  affect  the  drag  on  roughness  elements  in  a 
turbulent  boundary  layer. 
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Appendix  A  Trip  Arrangements 

A.l  Introduction 

In  order  to  obtain  the  desired  momentum  thickness  Reynolds  number,  Ree~7500,  and 
boundary  layer  thickness,  6^39  mm,  at  the  tunnel  measurement  location,  various  trip 
arrangements  were  studied  and  carefully  considered.  The  purpose  of  investigating  the 
numerous  and  varied  trip  arrangements  was  to  match  the  flow  conditions  in  the  Virginia 
Tech  Aerospace  and  Ocean  Engineering  (AOE)  Department  small  boundary  layer  tunnel 
with  the  conditions  already  present  in  the  Virginia  Tech  AOE  Department  larger 
boundary  layer  tunnel.  A  few  important  aspects  to  note  about  the  trip  arrangement 
selection  process  are  as  follows.  First,  momentum  thickness,  0,  needs  to  be  large  enough, 
yet  the  boundary  layer  thickness,  5,  and  shape  factor,  H=8*/0,  also  need  to  fit  into 
reasonable  bounds.  Another  factor  to  consider  is  that  the  boundary  layer  needs  to  be 
relaxed  enough  downstream  of  the  trip  to  be  a  two  dimensional  equilibrium  flow  at  the 
measurement  location. 

All  boundary  layer  profiles  were  taken  by  means  of  a  flattened-mouth  boundary  layer 
impact  probe  with  a  y-direction  probe  opening  equal  to  0.127  mm.  Profiles  were  taken  in 
even  logarithmically  spaced  y-direction  increments;  the  nearest  wall  reading  was 
assumed  to  be  at  half  of  the  outside  diameter  of  the  probe,  0.165  mm  from  the  wall.  The 
entire  setup  was  fastened  on  top  of  the  wind  tunnel  to  alleviate  any  wandering  during  the 
measurements.  The  vertical  traverse  system  is  powered  by  an  electric  motor  which  turns 
a  screw  with  a  thread  ratio  equal  to  40  threads/inch.  Pressure  readings  were  taken  with  a 
Series  475  Mark  III  Digital  Manometer  by  the  Dwyer  Co.  with  an  accuracy  of  ±0.01 
inches  of  water.  This  manometer  was  calibrated  and  checked  for  accuracy  with  an  actual 
inclined  manometer  before  measurements  were  taken.  The  flow  conditions  in  the  tunnel 
are  the  same  as  those  described  in  Chapter  2;  a  zero  pressure  gradient  flow  with  a  free- 
stream  speed  equal  to  27.5  m/s  at  a  temperature  of  25°C  ±0.56°C.  Since  this  is  a 
nominally  zero-pressure-gradient  flow,  the  location  of  the  static  pressure  port  will  not 
have  a  significant  effect  on  the  data.  Thus,  all  profiles  used  a  single  static  pressure  port 


Appendix  A:  Trip  Arrangements 


-227- 


that  is  described  in  Chapter  2.  All  profiles  are  measured  in  the  center  of  the  tunnel  at  port 
#20  unless  otherwise  specified.  Please  see  Chapter  2  for  actual  measurement  locations 
related  to  Port  #’s. 

The  results  of  the  measurement  system  used,  although  good  and  very  repeatable,  should 
only  be  taken  as  a  solid  starting  point  in  order  to  obtain  more  refined  data  with  the  LDV 
system.  At  least  one  point  in  the  viscous  sublayer  was  obtained  in  all  boundary  layer 
profiles  taken.  A  y-shift  value  was  applied  to  each  profile  and  was  found  in  an  iterative 
fashion  by  matching  the  viscous  sublayer  relation  of  U+=y+  for  y+<10  with  the  measured 
data.  The  experimental  setup  with  the  boundary  layer  probe  and  traverse  system  can  be 
seen  in  Figure  A.  1 . 

A  total  of  twelve  trip  configurations  were  studied  including  various  combinations  of  a 
grooved  plate,  sandpaper,  and  various  sized  square  bars.  Figures  of  the  trip  arrangements 
and  their  related  abbreviations  can  be  found  in  Figure  A.2  and  Table  A.l,  respectively. 
Please  note  that  in  Table  A.l,  the  word  sandpaper  is  abbreviated  and  denoted  by  s.p.  for 
brevity.  The  same  sandpaper  was  used  in  all  trip  arrangements.  The  sandpaper  was  20 
grit  silicon-carbide,  manufactured  by  Norton,  with  a  width  equal  to  8  inches.  The  final 
selected  trip  arrangement  was  Trip  Arrangement  1 1  (TA  1 1).  This  configuration  and  its 
relation  to  the  flow  field  is  discussed  in  more  detail  in  Chapter  2. 
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Table  A.l:  Trip  arrangements  and  related  abbreviations 


Description  of  Trip  Arrangement  (TA) 

Abbreviation 

13"  s.p.  on  floor 

TA  1 

13"  s.p.  on  floor  +  5"  s.p.  in  contraction 

TA  2 

18"  s.p.  on  floor  +  5"  s.p.  in  contraction 

TA  3 

Slotted  Plate  (19  slots) 

TA  4 

Slotted  Plate  (19  slots)  +  0.25"  bar 

TA  5 

Slotted  Plate  (19  slots)  +  0.125"  bar 

TA  6 

Slotted  Plate  (19  slots)  +  0.1875"  bar 

TA  7 

12"  s.p.  on  floor  +  0.1875"  bar 

TA  8 

13"  s.p.  on  floor  +  0.382"  bar 

TA  9 

18"  s.p.  on  floor  +  0.125"  bar 

TA  10 

18"  s.p.  on  floor  +  two  0.125"  bars 

TA  11 

18"  s.p.  on  floor  +  0.125"  and  0.1875"  bars 

TA  12 

A.2  Configurations  and  Related  Data 

For  each  trip  arrangement,  the  boundary  layer  thickness,  8,  momentum  thickness,  9, 
shape  factor,  H=5*/0,  and  the  friction  velocity,  UT,  are  calculated  and  presented  in  table 
format.  Here  8  is  defined  as  the  height  at  which  the  local  streamwise  velocity  is  equal  to 
0.99Ue.  For  each  measured  profile  the  momentum  thickness,  0,  and  the  displacement 
thickness,  8*,  were  calculated  using  the  following  equations,  respectively  (Kays,  1966). 


e=i\ -  —  \—dy 

0Jl  UJU. 


(A  A) 


u 

S--1 


u  v 


(A.2) 


o  v 

In  order  to  obtain  the  friction  velocity,  the  use  of  the  Ludwieg-Tillman  equation 
(Equation  A.3)  from  Pope  (2000)  for  skin  friction,  Cf,  was  employed. 


^  =  0.1 23(1 0)_o 


-0.678  H 


VJL 

K  v  ) 


-0.268 


(A.3) 
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This  equation  is  well  known  to  give  good  results  in  favorable  and  adverse  pressure 
gradient  flows  with  no  separation  for  103<Ree<30(103).  As  is  discussed  in  Chapter  2  this 
flow  is  a  zero  pressure  gradient  flow  with  an  Ree  within  these  bounds.  The  friction 
velocity  in  a  zero-pressure-gradient  2D  boundary  layer  can  then  be  calculated  using  the 
following  relation. 


V=U. 


(A.4) 


In  order  to  alleviate  unnecessary  measurements  for  each  trip  arrangement,  it  was  helpful 
to  calculate  the  momentum  thickness  along  the  streamwise  direction  using  known  data 
for  Ue(x)  and  0  at  one  x  location.  Kays  (1966)  method  for  finding  the  momentum 
thickness  at  various  streamwise,  x,  locations  within  a  turbulent  boundary  layer  was 
therefore  employed  and  can  be  seen  below. 


0*U 


4.11 


=  p(0-016)7e38V4dh: 


(A.5) 


\X0 

If  Ue(x)  is  measured  between  the  downstream  end  of  the  roughness  (x0)  and  the  desired 
location  of  boundary  layer  measurements  (xi),  then  one  can  calculate  0  at  any  point  x 
between  x0  and  xi,  due  to  the  fact  that  0  is  known  at  xi . 


Please  note  in  the  following  trip  arrangements  with  bars:  The  first  bar  listed  is  located  at 
the  contraction  exit  and  the  second  bar  is  always  located  2  inches  downstream  of  the 
contraction  exit.  In  all  trip  arrangements  with  square  bars,  the  bars  extend  the  entire  way 
across  the  tunnel  floor  and  are  fastened  to  the  floor  with  super  glue. 

A.2.1  Trip  Arrangements  1-3 

A  piece  of  silicon-carbide  Norton  20  grit  sandpaper,  8  inches  wide,  was  used  as  the  only 
method  for  tripping  the  boundary  layer  in  the  first  three  cases  analyzed.  This  8  inch  wide 
section  of  sandpaper  was  placed  in  the  center  of  the  9.5  inch  wide  test  section.  The  TA  1 
was  simply  a  section  of  sandpaper,  13  inches  in  streamwise  length,  attached  to  the  test 
section  floor  at  the  contraction  exit.  Rubber  cement  was  employed  to  attach  the 
sandpaper  to  the  center  of  the  tunnel  floor  for  its  durability  and  easy  cleanup 
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characteristics.  The  TA  2  and  TA  3  used  the  same  sandpaper  as  TA  1  but  the  sandpaper 
was  extended  5  inches  into  the  contraction  and,  in  the  case  of  TA  3,  5  inches  farther  into 
the  test  section.  The  results  for  these  trip  arrangements  can  be  seen  in  Table  A.2  and 
Figure  A. 3.  For  actual  x-locations  associated  with  the  Port  #’s  please  refer  to  Chapter  2. 


Table  A.2:  Results  from  trip  arrangements  1-3 


Trip  Arrangement 

C9ISB 

Ree 

H=8*/9 

UT  (m/s) 

TA  1 (Run  1) 

20 

0.0391 

6632 

1.36 

0.90 

TA  1  (Run  2) 

20 

0.0403 

6782 

1.37 

0.89 

TA  2 

20 

0.0339 

6339 

1.34 

1.01 

TA  3 

20 

0.0376 

6817 

1.33 

1.02 

It  is  important  to  note  that  the  boundary  layer  thickness  for  TA  2  and  TA  3  actually 
dropped  from  the  values  shown  by  TA  1.  This  is  due  to  the  fact  that  there  is  a  thickness 
to  the  sandpaper  which  in  turn  acts  as  a  trip  itself  when  placed  at  the  beginning  of  the  test 
section.  Whereas  for  TA  2  and  TA  3,  the  sandpaper  was  extended  into  the  contraction 
and  actually  provided  a  smoother  transition  into  the  test  section;  thus  the  initial  ‘front 
edge  trip’  was  not  felt,  which  decreased  the  boundary  layer  thickness.  The  shape  of  all 
boundary  layer  profiles  do  match  nicely  with  what  is  expected  thus  giving  added 
confidence  in  the  measurement  scheme. 

A.2.2  Trip  Arrangements  4-7 

Trip  arrangements  TA4-TA7  made  use  of  an  aluminum  plate  located  in  the  spanwise 
center  of  the  test  section  floor  with  the  following  dimensions:  12  inches  long,  8  inches 
wide,  and  0.125  inches  thick.  Spanwise  slots,  0.25  inches  wide  by  .0625  inches  deep, 
were  machined  across  the  entire  width  of  the  plate.  These  19  slots  started  3  inches  from 
the  front  of  the  plate  and  continued  downstream  with  a  spacing  equal  to  0.3125  inches. 
At  the  front  and  rear  sections  of  the  plate,  a  slope  of  7.13°  by  1  inch  long  was  machined 
across  the  width  of  the  plate.  This  ensured  smooth  transition  between  the  tunnel  floor 
and  the  slotted  plate.  This  trip  arrangement  was  placed  at  the  contraction  exit  and  was 
fastened  to  the  tunnel  floor  via  countersunk  screws.  TA  4  was  simply  the  slotted  plate 
itself.  TA  5-7  were  the  slotted  plate  along  with  a  0.25  inch,  0.125  inch,  and  0.1875  inch 
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square  bar,  respectively,  placed  directly  in  front  of  the  plate.  Results  from  the  boundary 
measurements  are  given  in  Tables  A. 3  and  A.4  and  Figures  A.4  and  A.5. 


Table  A.3:  Results  from  trip  arrangements  4-6 


Trip  Arrangement 

Ree 

H=8*/0 

Ux  (m/s) 

TA  4 

20 

0.0307 

5683 

1.34 

1.04 

TA  5 

20 

0.0516 

8466 

1.26 

1.04 

TA  6 

20 

0.0429 

7462 

1.31 

1.02 

Table  A.4:  Results  from  trip  arrangement  7 


Trip  Arrangement 

Reg 

H=8*/0 

RBI 

TA  7 

20 

0.0493 

8271 

1.31 

TA  7 

13 

0.0412 

7049 

1.28 

1.07 

TA7 

15 

0.0486 

7620 

1.28 

1.06 

As  would  be  expected,  TA  5  has  the  greatest  boundary  layer  thickness  due  the  larger  bar 
placed  in  the  flow.  Figure  A.4  also  shows  that  there  is  not  a  wake  region  present  for  this 
particular  flow  regime.  Measurements  were  taken  at  ports  13  and  15  because  the  actual 
measurement  location  for  the  LDV  measurements  is  between  these  two  ports;  see  Chapter 
2.  The  TA  7  is  close  to  the  desired  boundary  layer  properties,  but  the  shape  factor  is  not 
as  close  to  the  desired  value  of  H-1.34  and  the  boundary  layer  thickness  is  slightly  too 
high. 


A.2.3  Trip  Arrangements  8-9 

Trip  arrangements  8-9  both  utilized  the  same  sandpaper  as  discussed  previously.  The  TA 
8  has  12  inches  of  sandpaper  attached  to  the  tunnel  floor  along  with  a  0.1875  inch  bar 
placed  directly  in  front  of  the  sandpaper  at  the  contraction  exit.  The  TA  9  is  slightly 
different  in  that  it  has  13  inches  of  sandpaper  attached  to  the  tunnel  floor  along  with  a 
0.382  inch  bar  placed  directly  in  front  of  the  sandpaper  at  the  contraction  exit.  Results 
from  TA  8-9  can  be  seen  in  Table  A.5  and  Figure  A.6. 


Table  A.5:  Results  from  trip  arrangements  8-9 


BSMIMM 

Ree 

H=S*/0 

UT  (m/s) 

TA  8 

13 

0.0405 

6934 

1.28 

1.07 

TA  9 

20 

0.0546 

7379 

1.23 

1.06 
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Again,  these  configurations  yielded  a  significantly  lower  shape  factor  than  desired  along 
with  a  larger  boundary  layer.  As  can  be  seen  in  Figure  A.6,  TA  9  does  not  exhibit  a 
normal  wake  region,  thus  indicating  a  possible  merging  of  the  ceiling  and  floor  boundary 
layers. 

A.2.4Trip  Arrangements  10-11 

In  order  to  ensure  a  relaxed  and  two-dimensional  turbulent  boundary  layer  in  equilibrium, 
more  sandpaper  was  utilized  as  well  as  smaller  initial  bar  trips.  Sandpaper  is  an  excellent 
boundary  layer  trip  due  to  the  fact  that  it  has  a  significant  amount  of  small  randomly 
distributed  roughness  elements  that  produce  drag  and  ensure  good  mixing.  To  this  end, 
TA  10  consists  of  18  inches  of  the  same  sandpaper  discussed  previously  along  with  a 
0.125  inch  square  bar  placed  directly  in  front  of  the  sandpaper  at  the  contraction  exit.  TA 
1 1  is  exactly  the  same  as  TA  10  but  with  another  0.125  inch  square  bar  placed  2  inches 
downstream  of  the  contraction  exit  on  the  tunnel  floor.  The  placement  of  the  second 
0.125  inch  square  bar  was  determined  by  giving  the  flow  sufficient  distance  to  reattach 
before  applying  another  step.  Results  from  the  boundary  measurements  can  be  seen  in 
Table  A.6  and  Figure  A.7. 


Table  A.6:  Results  from  trip  arrangements  10-11 


Trip  Arrangement 

Ree 

H=s*/e 

UT  (m/s) 

TA  10 

14 

0.0313 

6181 

1.33 

1.04 

TA  1 1 

14 

0.0350 

6500 

1.33 

1.04 

Measurements  of  these  trip  arrangements  were  made  at  Port  #14  due  to  the  fact  that  this 
location  is  the  closest  to  the  actual  measurement  location  for  the  LDV  data.  Figure  A.  7 
shows  the  boundary  layer  profiles  for  TA  10-11  to  be  smooth  as  well  as  have  well 
defined  log  and  wake  regions.  The  shape  factor  for  these  configurations  is  very  near  the 
desired  value  of  H;=1.34.  The  boundary  layer  thickness  and  Ree  are  slightly  low,  but  as 
can  be  seen  in  Chapter  2,  more  refined  LDV  measurements  give  the  desired  results  in  all 
areas  of  describing  the  boundary  layer.  Thus,  TA  1 1  was  the  final  trip  arrangement  that 
was  selected  for  use  in  all  measurements  of  the  roughness  elements. 
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A.2.5  Trip  Arrangement  12 

The  TA  12  is  identical  to  TA  1 1  except  the  second  square  bar  in  TA  1 1  is  replaced  with  a 
0.1875  inch  square  bar  in  TA  12.  The  results  from  the  TA  12  boundary  layer  profiles  can 
be  seen  in  Table  A. 7  and  Figure  A.8. 


Table  A.7:  Results  from  trip  arrangement  12 


Trip  Arrangement 

Location  (Port  #) 

5(m) 

Ree 

H=8*/0 

UT  (m/s) 

TA  12 

14 

0.0383 

6728 

1.31 

1.05 

TA  12 

15 

0.0379 

6610 

1.31 

1.05 

TA  12 

20 

0.0428 

7892 

1.32 

1.01 

TA  12  gave  results  that  were  very  similar  to  TA  1 1,  the  most  noticeable  differences  being 
in  the  shape  factor,  H,  and  boundary  layer  thickness,  8.  Considering  all  effects,  TA  1 1 
and  TA  12  are  remarkably  close  as  should  be  expected  due  to  the  fact  that  only  a  slight 
difference  exists  in  the  second  bar  for  each  case. 
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A.3  Appendix  A  Figures 


Figure  A.l:  Boundary  layer  impact  probe  and  traverse  system 
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13  inches  of  sandpaper  (TA  1) 


Slotted  plate  +  1/4  inch  square  bar  (TA  5)  cont. 
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Slotted  plate  +  1/8  inch  square  bar  (TA  6) 


Slotted  plate  +  3/16  inch  square  bar  (TA  7) 


12inches  of  sandpaper  +  3/16  inch  square  bar  (TA  8)  coni. 
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13  inches  of  sandpaper  +  0.382  inch  square  bar  (TA  9) 


18  inches  of  sandpaper  +  two  1/8  inch  square  bars  (TA  11)  coni. 
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1 8  inches  of  sandpaper  +  1/8  inch  square  bar  +  3/16  inch  square  bar  (TA  12) 


Figure  A.2:  Tested  trip  arrangements:  TA  1-TA  12  (Note  that  TA  2  and  TA  3  are 

not  shown) 


-B -  13"  of  sandpaper~(mn  1) 

A -  1 3”  of  sandpaper— {run  2) 


Figure  A3:  Streamwise  velocity  profiles  of  trip  arrangements  TA  1  -  TA  3 
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-E3 -  Slotted  plate 

■A -  Slotted  plate  +  0.25"  bar 


Figure  A.4:  Streamwise  velocity  profiles  of  trip  arrangements  TA  4  -  TA  6 


■B -  Slotted  plate  +  0.1 875”  bar 

A -  Slotted  plate  +  0.1875”  bar— (port  13) 


Figure  A.5:  Streamwise  velocity  profiles  of  trip  arrangement  TA  7 
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Appendix  B  Single  Element  Drag  Analysis 


B.l  Derivation  of  Drag  Equation 

The  following  derivation  of  the  equation  used  to  compute  the  drag  on  the  single  isolated 
roughness  elements  is  taken  directly  from  George  and  Simpson  (2004).  The  control 
volume  (CV)  used  in  this  analysis  can  be  found  in  Figure  B.l.  The  CV  encloses  one  half 
of  the  element  domain;  thus,  it  is  important  to  keep  in  mind  that  the  drag  obtained  will  be 
Vi  of  the  total  drag  due  to  the  element.  In  order  to  obtain  the  drag  due  to  the  element,  one 
needs  to  start  from  the  mass  conservation  and  momentum  principles.  Each  principle  is 
applied  to  the  CV  with  the  element  (E)  and  again  without  the  element  (w/oE).  The  drag 
due  to  the  element  is  then  determined  to  be  the  difference  between  the  individual  drag 
terms  with  and  without  the  element.  Please  note  that  the  subscripted  numbers  in  the 
equations  refer  to  the  respective  face  of  the  control  volume. 

The  CV  is  selected  so  that  one  side.  Face  4,  is  along  the  line  of  symmetry  of  the  element 
in  order  to  obtain  a  spanwise  velocity  component  equal  to  zero  on  this  face. 

Consequently,  V.n  =  0  on  this  particular  face  which  simplifies  the  analysis, 
where  V  =  (U  +  u)i  +  (V  +  v)j  +  (W  +  w)£.  The  length  of  the  CV  is  taken  to  be  7.75d, 
where  Face  2  is  located  a  distance  of  2.75d  downstream  of  the  element  and  Face  1  is 
taken  to  be  at  the  upstream  side  of  the  element  and  is  referred  to  as  -2D,  i.e.  two- 
dimensional  flow  at  a  negative  x-direction  location  (see  Figure  B.l).  Full  pressure 
recovery  is  expected  for  all  elements  at  this  downstream  distance  based  on  previous  work 
done  by  Martinuzzi  and  Tropea  (1993).  From  this  analysis  it  was  shown  that  in  the  flow 
field  downstream  of  a  cube,  there  was  full  pressure  recovery  at  a  distance  equal  to  2.25d. 
Therefore,  full  pressure  recovery  can  be  expected  in  all  of  the  cases  under  investigation 
due  to  the  fact  that  the  other  elements  should  recover  even  sooner  than  the  cube.  The 
location  of  Face  3,  -Z2D,  is  taken  as  the  distance  at  which  the  boundary  layer  sees  no 
effect  from  the  element’s  presence  and  is  the  same  as  the  smooth  wall  2D  boundary  layer. 
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The  location  of  the  top  face  of  the  CV,  ytoP,  is  located  at  a  distance  far  enough  from  the 
wall  where  the  U  and  V  velocity  values  show  no  effect  of  the  element  being  present.  For 
consistency,  this  distance  is  taken  to  be  the  boundary  layer  thickness,  8. 


Start  by  applying  the  steady  flow  conservation  of  mass  principle  to  the  CV  with  and 
without  the  element. 

$p(v.n}iA  =  0,  (B.l) 

where  dA  is  the  elemental  area. 

After  applying  Equation  (B.l)  to  the  CV  with  and  without  the  element,  the  time  averaged 
results  are  Equations  (B.2)  and  (B.3)  below  respectively. 

With  the  element  present,  Eqn.  (B.2): 


y«>p 

}’iop  ~Z2  p 

y,op  2.1  Sd 

-Z2D  21Sd 

p  J  ju^dzdy 

+  p  J  | U2dzdy 

+  p  j  ^W^dxdy 

+  P  |  fKopdxdz 

0  0 

£.  1  0  0 

£.2  0  -2D 

E%  3  0  -2D 

Without  the  element  present,  Eqn.  (B.3): 


y?op  ~z2p 

ytop  -z2o 

>v  2.15 d 

-Z2D  2.15d 

p  |  J U,dzdy 

+  p  J  | U2dzdy 

+  p  j  ^W^dxdy 

+  P  J  J Kopdxdz 

0  0 

wloE,  1  0  0 

w/oE.2  0  ~2D 

w,oE3  0  -2D 

Next,  apply  the  momentum  principle  to  the  CV  knowing  that  the  change  in  the  x- 
momentum  is  balanced  with  the  sum  of  the  forces  in  the  x-direction. 

-  D,  (B.4) 


Note  that  dA  is  the  elemental  area,  da  is  the  elemental  volume,  and  Dx  is  the  x- 
component  of  the  drag  acting  on  the  control  volume.  Following  the  time-averaging  of 
Eqn.  (B.4),  we  get: 

With  the  element  present,  Eqn.  (B.5): 


ytop  v  y<op  z2  a,  v  y>op  2.i5d  _ , 

-p  \  J p'+ufjjzdy  +p\  J\C22  +  u22 jizdy  +p  J  J(f/3fF3  +w3w3 jixdy 


0  0 


I  £.1 


0  0 


Ef  2 


0  -2  D 


E,  3 


+  p]  \tu,oPKoP  +u,oPv,oP}ixdz  = -^[Volumecv\  -Dx 

0  ~2D  E,top  'E 


Without  the  element  present,  Eqn.  (B.6): 
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yiop-Z2Ix  v 

-p\  \P\2 +u\2pzdy 

o  o 

-Z2D2JSd  _ 

+  P  J  \V'opVtop+U,op 

0  -2  D 

Now,  subtract  Eqn.  (B.6)  from  Eqn.  (B.5)  to  obtain  the  following  equation,  Eqn.  (B.7). 

Keep  in  mind  that  the  U  and  u 2  profiles  are  identical  on  the  front  face  of  the  CV,  Face  1. 
Thus,  the  two  instances  with  and  without  the  element  will  cancel  each  other.  For  reasons 
discussed  previously,  the  pressure  gradient  in  the  streamwise  direction  is  also  negligible. 
This  in  turn,  reduces  the  respective  term  to  zero. 


y>op  -z2  q. 

*p\  jk 

+  u22^lzdy 

y,op  2.75 d 

+  p\  J(^3^3  +uiwipxdy 

/Cl  00 

wl  oE,  1 

w/oE,2  0  ~2D 

y»op 

-p\  JV* 


-U2  pzdy 


\z 


0  0 


y lop  ~z2d 


-p\ \  “ 


- U o 


v!  oE 


fizdyl 


E  J 


y«>p  2.1  Sd 


-P  J 


0  -2D 


-pTUI.  - W|« W  -e 


0  -2D 
-Z2D  2.75, 


=p  )  it"), L«-^vU«)tefe 


0  -2D 


— D  c~\~D  ,  F 

x  E  x  wloE 


top 


(B.7) 


The  next  step  in  the  analysis  is  to  multiply  Eqns.  (B.2)  and  (B.3)  by  their  respective  U 
velocity  at  the  top  face  of  the  CV,  Utop.  Following  this  step,  rearrange  to  obtain  Eqns. 
(B.8)  and  (B.9)  below. 


U„ 


-z2D  2.75 d 

y<op  “22£> 

yiop 

ytop  2.75  d 

\e 

P  f  \V,opdxdz 

=  p  |  \l/ldzdy 

-  p  \  ]u2dzdy 

-  p  J  jfV}dxdy 

0  -2D 

E.iop  ®  ® 

E.\  0  0 

E.2  0  -2D 

£,3_ 

(B.8) 


-zlD  2JSd 

y top  ~z2d 

ytop  ~%2d 

y,op  2.15d 

l  w/  oE 

P  \  \V.oPdxdz 

=  p\  \Uxdzdy 

—  p\  ^U2dzdy 

-  p  j  ^W2dxdy 

0  -2D 

w/oEjop  ®  ® 

w/oE.i  0  0 

v!oE.l  0  -20 

w/ oE,  3 

(B.9) 

Due  to  the  fact  that  the  respective  streamwise  velocities  at  the  top  of  the  CV,  Utop,  are 
constant,  we  can  substitute  Eqns.  (B.8)  and  (B.9)  into  the  fifth  term  on  the  LHS  of  Eqn. 

(B.7)  above.  It  is  also  of  importance  to  note  that  U  top\  /  =  Ulop  |  ^  for  all  faces  of  the  CV 
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and  for  the  remainder  of  the  derivation  will  simply  be  called  Ulop .  In  doing  so,  the 
following  result  is  obtained: 


Stop  —Zirf  \  *‘°p  ~Z2D  / _  _  \  y>op  2JSd 

-p  \  WA„,0E-UA>dy  ~p  \  J  (M^2  -u22  Uzdy  -p  \  J((t/3^)|w/o£-(t/3^3)|£)irrf: 

nn  n  n  V  w/°E  EJ  n  n 


y>op  2.15d  -z2D2J5d _  _  X 

■u~p  1  -p  J 


(B.10) 

In  order  to  simplify  Eqn.  (B.10),  it  is  significant  to  note  that  for  large  distances  away 
from  the  centerline,  the  flow  will  be  two-dimensional  and  thus  there  will  be  no  spanwise 
velocity  component.  Therefore,  terms  3,  4,  and  7  will  be  zero.  Term  5  will  also  reduce 
to  zero  due  to  the  fact  that  the  U  velocity  profiles  on  Face  1  are  identical  for  both  cases  as 
was  discussed  previously.  Term  8  will  also  reduce  to  zero  due  to  the  fact  that  at  the  top 

of  the  CV,  (wv),op  |  =  ( uv\op  .  Thus,  it  can  be  seen  then  that  the  only  velocity  profiles 

needed  to  calculate  the  drag  are  the  profiles  located  at  Face  2.  Following  the 
simplification  we  arrive  at  Eqn.  (B.l  1)  below: 


y top  \  y<op  -Zip  / _  \ 

&de  =  p  f  }V22  .  ~ui2\yzdy  +  p  f  f  u22  ~u22  fzdy 

0  0  0  E  0  0  V  wtoE  E) 

y<op  ~Z7D.  . 

-UmP\ 


(B.  11) 


As  discussed  previously,  the  total  drag  due  to  the  element’s  presence  in  the  boundary 
layer  can  now  be  defined  as  twice  that  of  Eqn.  (B.l  1).  The  drag  due  to  the  element  can 
also  be  presented  as  a  drag  coefficient  in  the  following  way. 


„  _  2A De 

L' d  ~  . 

-pU]A 

2 


WlVL  -M.H-.'L  -u^zdy 


—pu,A 


(B.12) 
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This  drag  coefficient  is  normalized  by  the  tree-stream  dynamic  pressure  multiplied  by  an 
elemental  area  term.  This  area  is  taken  to  be  either  the  frontal  or  wetted  area  of  the 
element,  Af  or  Aw.  Note  that  a  designation  is  made  for  the  frontal  or  wetted  area  when 
presented. 
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B.2  Appendix  B  Figures 


Top  View  of  Control  Volume 
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Figure  B.l:  Views  of  control  volume  used  in  drag  analysis 


Appendix  C  Uncertainty  Analysis 

C.l  Wall  Location 

An  iterative  procedure  was  utilized  in  order  to  find  the  location  of  the  wall  for  each 
profile.  A  spectrum  analyzer  and  an  oscilloscope  were  used  to  find  the  wall  position  with 
respect  to  the  measurement  volume.  When  the  measurement  volume  is  centered  on  the 
surface  of  the  glass  disc  (wall)  a  strong  signal  is  sent  to  the  spectrum  analyzer  and 
subsequently  to  the  oscilloscope  for  viewing.  Careful  examination  of  the  calculated  y- 
shift  values  from  all  profiles,  the  traverse  system,  and  the  size  of  the  measurement 
volume,  leads  to  the  conclusion  that  the  uncertainty  is  on  the  order  of  the  measurement 
volume  diameter.  Consequently,  the  uncertainty  is  ±50  pm. 


C.2  Alignment  of  Probe  with  Element 

The  alignment  of  the  LDV  probe  with  the  center  of  each  element  was  done  in  an  iterative 
fashion,  and  the  procedure  is  discussed  in  Chapter  2,  Section  2.6.2.  By  watching  the 
beam  profiles  change  along  the  sides  of  the  element  and  the  wall  of  the  tunnel  (when  the 
measurement  volume  was  higher  than  the  element  height)  it  was  possible  to  have  an 
uncertainty  in  positioning  of  ±60  microns  in  the  x  and  z  directions. 


C.3  Velocities 

The  uncertainties  associated  with  this  LDV  system  have  been  evaluated  by  Kuhl  (2001) 
and  Tang  (2004).  Through  these  previous  investigations,  the  biases  associated  with  the 
alignment  of  the  probe  and  the  measurement  of  the  beam  angles  were  found  to  be 
negligible.  The  fringe  bias,  geometric  bias,  angular  bias,  and  velocity  gradient  bias  are 
negligible  due  to  the  fact  that  the  measurement  volume  of  the  LDV  is  small  and  nearly 
spherical.  These  prior  studies  also  show  that  the  velocity  bias  is  negligible  for  this 
system.  The  measured  signal  and  related  data  processing  of  that  signal  can  include 
different  uncertainty  effects.  Likewise,  Tang  (2004)  discovered  that  the  effect  of  finite 
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transit  time  broadening  on  the  LDV  was  negligible.  Olcmen  et  al.  (1998)  investigated  the 
instrument  broadening  effects  of  the  LDV  system.  The  same  Macrodyne  signal 
processors  that  were  used  in  this  research  were  analyzed  and  proven  to  have  negligible 
bandwidth  broadening  effects. 

Once  the  beam  angles  were  measured  for  the  current  setup,  a  two-dimensional  boundary 
layer  profile  was  taken  in  order  to  validate  the  LDV  system.  Near-wall  DNS  data  as  well 
as  the  law  of  the  wall  for  a  two-dimensional  turbulent  boundary  layer  were  used  to 
validate  the  LDV  measurements.  As  is  shown  in  Figure  2.1 1,  the  systematic  error  for  the 
LDV  system  has  been  taken  into  account  and  alleviated.  Consequently,  the  uncertainty 
associated  with  the  LDV  system  was  analyzed  by  looking  at  the  measurement  procedure 
as  a  whole.  The  uncertainties  in  the  velocities  obtained  from  the  LDV  system  are 
calculated  following  the  same  procedure  used  by  Olcmen  and  Simpson  (1995).  Two 
profiles  were  taken  on  two  consecutive  days  at  the  same  location  inside  the  wind  tunnel. 
The  location  was  116.72  cm  downstream  of  the  contraction  exit  and  at  the  spanwise 
center  of  the  wind  tunnel  test  section.  The  profiles  taken  were  in  a  smooth  wall  turbulent 
boundary  with  no  roughness  element  present.  The  uncertainties  of  the  velocity  values  are 
calculated  using  Chauvenet’s  criterion  to  calculate  the  standard  deviation.  Equation  C.l 
was  used  to  give  21:1  odds  uncertainties  in  the  data. 

—JSSL  =  1.15  (C.l) 

<7 

Here,  dmax  is  the  average  of  the  half  of  the  differences  between  two  data  values  and  a  is 
the  standard  deviation.  The  constant  1.15  is  chosen  to  give  21:1  odds  uncertainties 
calculated  as  ±2a.  The  calculated  uncertainties  can  be  seen  in  Table  C.l  below. 
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Table  C.l:  Uncertainties  in  measured  quantities 


u/u, 

±0.28 

u\!u) 

±0.19 

v/u. 

±0.07 

u2w/ul 

±0.19 

w/u, 

±0.17 

v2w/u 2 

±0.03 

u2/u; 

±0.51 

±0.07 

7/ul 

±0.06 

uw2  ju] 

±0.18 

w'/u; 

±0.11 

vw2  ju] 

±0.04 

mju) 

±0.05 

uvw/u ] 

±0.05 

uw/u2 

±0.07 

s/v; 

±0.79 

vw/U2 

±0.03 

s/u; 

±0.06 

w’/u; 

±0.64 

C.4  Turbulent  Kinetic  Energy 

The  turbulent  kinetic  energy  (TKE)  is  calculated  using  the  equation 

TKE  =  \^ll)  =  If7  +  v2  +  ^7  u  ,  where  u2 , v2 ,w2  are  the  streamwise,  normal-to-wall, 
and  spanwise  Reynolds  normal  stresses  respectively.  Through  the  application  of  the  root- 
mean-square  (RMS)  error  analysis  technique,  the  uncertainties  are  calculated  using 
Equation.  C.2, 


STKE  =  ± 


/  \2 
'dTKE' 


du2  j 


Su2  + 


/  \2 
dTKE ' 


dv 2 


Sv 2  + 


dTKE 

dw2 


\2 


Sw2 


(C.2) 


where  Su2  ,Sv2  ,Sw2  are  the  associated  uncertainties  of  the  respective  Reynolds  normal 
stresses.  These  uncertainty  values  can  be  found  in  Table  1  above.  Taking  the  derivatives 

and  substituting  the  necessary  uncertainties  gives  STKE  =  S^2  / 2]  =  ±0.265  m2/s2. 


C.5  Skin  Friction  Velocity 

The  error  in  the  skin  friction  velocity,  UT,  at  20:1  odds  was  obtained  using  the  following 
equation,  Equation  C.3, 
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SU,  = 


2.365o- 

■Tn 


(C.3) 


where  n  is  the  number  of  identical  profiles  taken  and  a  is  the  standard  deviation  between 
profiles.  The  total  number  of  two-dimensional  profiles  taken  that  were  used  to  evaluate 
the  uncertainty  was  7.  The  profiles  were  all  taken  at  the  same  location  on  different  days. 
The  value  of  2.365  was  obtained  using  the  student’s  t-distribution  table  for  a  95% 
confidence  interval.  Following  the  substitution  of  values,  the  uncertainty  in  UT  =  ±5%. 


C.6  Flow  Angle  Data 

The  mean  flow  angle  (FA)  is  calculated  via  FA  =  tan_l(lF/f/),  where  W  is  the  mean 
spanwise  velocity  and  U  is  the  mean  streamwise  velocity.  The  uncertainties  in  the  FA 
can  be  calculated  using  the  RMS  error  analysis  technique  seen  in  Equation  C.4, 


SFA  =  ±, 


dFA 

dU 


^  f  dFA'' 


SU1  + 


dW 


5W 2 


(C.4) 


where  8U  and  8W  are  the  uncertainties  in  the  mean  velocities.  In  this  case  8U=0.266  m/s 
and  8W=0.165  m/s.  The  U  and  W  values  are  taken  to  be  U=13.582  m/s  and  W=0.034 
m/s  and  were  obtained  by  taking  an  average  over  all  single  element  profiles.  Therefore, 
taking  the  necessary  derivatives  and  substituting  in  the  proper  values  the  uncertainty  in 
the  flow  angle  FA  is  8FA  =  ±0.695°.  Thus,  the  flow  angle  uncertainty  is  taken  to  be  8FA 
k  ±0.7°. 

The  shear  stress  angle  (SSA)  is  calculated  using  the  equation,  SSA  =  tan~1(-vw/-«v). 
Following  the  same  principal  as  above,  the  uncertainties  in  the  SSA  are  calculated  using 
Equation  C.5  below. 


dSSA  =  ± 


dSSA 
V  dvw  ) 


Svw  + 


dSSA 


v  duv  ) 


Suv 


(C.5) 


The  uncertainties  in  uv  and  vw  are  as  follows,  Jwv=0.050  (m/s)2  and  Jvvv  =0.030 


(m/s)2.  Similar  to  above,  the  wv  and  vw  values  are  taken  to  be  uv  =  -0.725  (m/s)2  and 
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vw=  -0.008  (m/s)  and  were  obtained  by  taking  an  average  over  all  single  element 
profiles.  After  performing  the  necessary  steps,  8SSA  =  ±2.37°  ~  ±2.4°. 


C.7  Drag 

The  uncertainty  of  the  drag  on  an  isolated  roughness  element  was  calculated  via  a 
perturbation  (jitter)  analysis  performed  on  the  data  reduction  code  used  to  calculate 
Equation  4. 1 .  The  standard  deviations  of  the  mean  streamwise  U  velocity  as  well  as  the 

streamwise  Reynolds  normal  stress  u 2  were  calculated  using,  £(  )/1.96,  where  )  is 

the  uncertainty  associated  with  either  U  or  u2  .  Vectors  of  random  numbers  were 
generated  from  a  normal  distribution  with  a  mean  equal  to  zero  and  a  standard  deviation 
equal  to  one.  These  random  vectors  had  lengths  that  were  equal  to  the  interpolated 
velocity  profiles.  These  vectors  were  then  multiplied  by  the  respective  standard 

deviations  of  U  and  u 2  and  added  to  the  original  U  and  u 2  velocity  vectors  in  order  to 
obtain  a  ‘perturbed’  vector.  The  uncertainty  in  the  measurement  location  was  also  taken 
into  account  via  the  same  process  as  above,  with  one  exception.  The  uncertainty 
associated  with  the  y  and  z  positions  is  only  relative  to  the  wall  and  centerline, 
respectively.  Thus,  all  y  and  z  position  values  were  perturbed  as  a  whole.  This  is  in 
contrast  to  the  above  mentioned  method  in  that  all  of  the  velocity  values  were  perturbed 
with  respect  to  each  other. 

The  same  procedure  that  was  used  to  calculate  Equation  4. 1  was  performed  using  the  new 
‘perturbed’  vectors.  This  process  was  repeated  ten  times  and  an  average  and  standard 
deviation  value  were  obtained  from  these  runs.  From  this  analysis,  the  uncertainty  in  the 
calculated  drag  was  determined  to  be  ±6%. 


Appendix  D  Gaussian  Fetch  Data 

D.l  Introduction 

This  appendix  contains  a  brief  analysis  and  discussion  of  the  data  obtained  on  the 
Gaussian  fetch  of  roughness.  More  in  depth  analysis  will  be  conducted  in  the  future. 
Figures  of  the  related  roughness  data  can  be  found  in  the  figures  at  the  end  of  the 
appendix. 


D.2  Brief  Analysis 

The  flow  conditions  within  the  tunnel  for  both  the  fetch  of  Gaussian  roughness  and  the 
isolated  roughness  element  cases  were  the  same,  i.e.  dp/dx=0,  Ue=27.5  m/s.  The 
boundary  layer  thickness,  8,  at  the  measurement  region  for  the  fetch  of  Gaussian 
roughness  was  52.2  mm.  The  friction  velocity,  UT>  was  found  by  assuming  a  constant 
shear  stress  region  and  calculated  using  Equation  D.  1 . 

u.-f^  (D.l) 

The  constant  Reynolds  shear  stress  values  in  the  log  region  are  averaged  to  give  uvavg  . 

The  actual  value  obtained  for  the  Gaussian  fetch  regime  was  UT  =1.43  m/s.  This  value  of 
UT  was  used  to  non-dimensionalize  the  measured  velocity  quantities  presented  for  the 
Gaussian  fetch  of  roughness.  The  non-dimensionalized  roughness  height  based  on  UT  for 
the  Gaussian  fetch  is  k+=217. 

As  can  be  seen  in  Figure  1,  the  U4  velocity  profiles  along  the  centerline  of  the  element 
show  a  drastic  difference  as  compared  to  the  two-dimensional  smooth  wall  case.  There  is 
not  only  a  AU4  shift  downward,  but  also  a  reduction  in  the  log  layer  as  well.  The  U+ 
profiles  tend  to  collapse  to  one  another  just  after  the  height  of  the  roughness.  A  pitot  tube 
boundary  layer  profile  was  also  taken  at  the  centerline  location  of  x/d=1.08  in  order  to 
obtain  the  boundary  layer  thickness  and  to  compare  with  the  LDV  data.  The  V+  profiles 
show  a  slight  peak  in  downwash  just  above  the  element  heights.  The  profiles  closest  to 
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the  element  show  that  there  is  a  positive  V+  velocity  directing  fluid  away  from  the  wall 
below  the  height  of  the  element.  The  Reynolds  normal  and  shearing  stresses, 

u2 ,  v2  ,w2  ,uv  show  peaks  in  magnitudes  at  or  just  above  the  height  of  the  element.  This 
is  due  to  the  fact  that  there  is  lower  momentum  fluid  near  the  wall  and  around  the 
roughness  elements.  Thus,  higher  Reynolds  stresses  will  be  produced  just  above  the 
height  of  the  element  where  the  slower  near  wall  fluid  meets  the  faster  out  region  fluid. 
This  interaction  will  create  a  shear  layer  at  the  interface  of  these  two  regions  that  will 
increase  turbulence  quantities.  The  related  contour  plots  of  the  Reynolds  normal  and 
shear  stresses  at  the  center  plane  between  elements  (x/d=1.08).  Figures  D.5  through  D.8, 
also  show  this  same  phenomenon.  These  plots  show  a  core  of  peak  values  near  the 
roughness  element  height  which  radiate  outward.  Although  the  magnitudes  are  smaller 
for  the  fetch  of  roughness,  this  trend  is  comparable  to  the  single  element  case. 


The  secondary  velocity  vectors  at  the  center  plane  location  between  elements,  x/d=1.08, 
are  presented  on  the  U+  contour  plot  in  Figure  D.2.  These  flow  vectors  show  fluid  being 
swept  in  towards  the  wall  and  around  the  element  at  this  center  plane  location.  Contours 


of  normalized  turbulent  kinetic  energy 


(TKE),  iF/2y< 


,  at  the  x/d=1.08  plane  can  be 


seen  in  Figure  D.21 .  There  is  a  high  region  of  TKE  near  the  element  height.  These  peak 
TKE  values  for  the  fetch  are  significantly  lower  than  the  values  for  the  single  element 

case.  The  normalized  diffusion  velocity  vectors  Vq  /Ur  =  (v<?2/ q2  ]/UT  and 

are  also  presented  in  Figure  D.21.  These  vectors  show  the 


transport  of  TKE  to  be  much  smaller  than  the  single  element  case.  There  is  still 
dissipation  of  TKE  away  from  the  core  value  near  the  height  of  the  element  but  the  vector 
magnitudes  and  directions  are  not  nearly  as  distinct  as  they  are  for  the  single  Gaussian 
element  case.  This  incident  is  not  surprising  due  to  the  fact  that  there  is  lower 
momentum  fluid  down  near  the  roughness  elements  in  the  fetch  regime.  Although  there 
are  significant  differences  between  the  single  roughness  element  case  and  the  fetch  of 
roughness,  there  are  also  distinct  similarities.  This  leads  to  the  conclusion  that  a  thorough 
knowledge  of  each  case,  single  element  and  fetch,  is  necessary  in  order  to  completely 
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understand  the  effects  that  roughness  has  on  the  flow  physics  in  a  turbulent  boundary 
layer. 
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D.3  Appendix  D  Figures 


(cont) 


u3/  Ur3  uw^Ut3  u2v/Ut 


Figure  D.l:  Centerline  profiles  of  Gaussian  fetch  (mean  velocities,  Reynolds  stresses,  and  triple 

products) 
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Figure  D.2:  Contours  of  U/UT  in  the  y-z  plane  at  a  location  of  x/d=1.08  with 
normalized  secondary  velocity  vectors  V/UT  and  W/UT 
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Figure  D.3:  Contours  of  V/Ur  in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.4:  Contours  of  WjUT  in  the  y-z  plane  at  a  location  of  x/d=l. 08 
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Figure  D.5:  Contours  of  u1  ju]  in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.6:  Contours  of  v2/i7r2 


in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.7:  Contours  of  w2/u 2 


in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.8:  Contours  of  uv/u2 


in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.9:  Contours  of  uw/lf 2 


in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.10:  Contours  of  vw/U r2  in  the  y-z  plane  at  a  location  of  x/d=1.08 
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FigureD.il:  Contours  of  u 2 v/lf*  in  the  y-z  plane  at  a  location  of  x/d=l. 08 
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Figure  D.12:  Contours  of  u2w/ul  in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.13:  Contours  of  v2w/u ,  in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.14:  Contours  of  uv2 fu] 
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Figure  D.15:  Contours  of  uw2/u] 
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Figure  D.16:  Contours  of  vw2 


in  the  y-z  plane  at  a  location  of  x/d=1.08 
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in  the  y-z  plane  at  a  location  of  x/d=1.08 


Figure  D.17:  Contours  of  uvw/u 3 
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Figure  D.18:  Contours  of  «3 /u] 


in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.19:  Contours  of  v3/t/3  in  the  y-z  plane  at  a  location  of  x/d=1.08 


Appendix  D:  Gaussian  Fetch  Data 


-267- 


Plane  at  x/d=1.08 


Gaussian 


0.5  0 

Fetch-k=2.54  mm 


-0.5 


4.00 

3.00 

2.00 

1.00 

0.00 

-1.00 

-2.00 

-3.00 

-4.00 

-5.00 


Figure  D.20:  Contours  of  wi/u3r 


in  the  y-z  plane  at  a  location  of  x/d=1.08 
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Figure  D.21:  Contours  of  TKE 

V,/U,  and  WjU, 


with  diffusion  velocity  vectors, 
=  {wq*  /q2]/UT ,  in  the  y-z  plane  at  x/d=1.08 


Appendix  E  Flow  Angle  Data 

E.l  Mean  Flow  Angle 

The  mean  flow  angle  (FA)  is  calculated  via  Equation  E.l  below,  where  W  is  the  mean 
spanwise  velocity  and  U  is  the  mean  streamwise  velocity. 

(  O'} 

FA  =  tan-1  —  (E.l) 

The  flow  angle  enables  one  to  determine  the  degree  of  cross-flow  (in  the  spanwise 
direction)  due  to  three-dimensional  effects.  If  the  flow  angle  is  larger  than  45°  at  any 
point  within  the  flow,  the  cross-flow  is  deemed  to  be  locally  larger  than  the  flow  in  the 
streamwise  direction.  The  calculated  uncertainty  in  the  flow  angle  is  5FA  ~  ±0.7°;  see 
Appendix  C,  Section  C.6.  Figure  E.l  shows  the  mean  flow  angle  variation  along  the 
centerline  of  the  elements.  The  mean  flow  angle  is  expected  to  be  zero  along  the 
centerline  due  to  symmetry  reasons  within  the  flow.  Within  the  bounds  of  the 
uncertainty,  this  is  generally  the  case  at  all  streamwise  centerline  measurement  locations 
except  for  the  first  location  at  x/d=1.36.  Due  to  the  proximity  of  this  location  with 
respect  to  the  element  and  the  highly  turbulent  flow  that  is  present;  even  the  slightest 
error  in  probe  alignment  will  produce  results  that  deviate  from  the  expected  zero  flow 
angle  value.  As  is  shown  at  the  subsequent  downstream  measurement  locations,  these 
small  alignment  errors  do  not  propagate  very  far  downstream.  Figures  E.2  and  E.3  show 
contour  plots  of  the  flow  angles  at  the  first  measurement  plane  downstream  of  the 
elements.  All  axisymmetric  elements  show  FA  values  that  are  very  comparable  with 
respect  to  each  other.  At  the  measurement  location,  x/d=2.75,  mean  U  is  positive 
throughout  all  boundary  layer  profiles.  Consequently,  negative  flow  angles  indicate  that 
the  flow  is  being  directed  in  the  (-)z-direction,  (-)spanwise  direction,  where  the  positive 
flow  angle  values  indicate  the  opposite  flow  behavior.  The  flow  angles  are  related  to  the 
strength  of  the  secondary  flow  and  thus  directly  related  to  the  strength  of  the  horseshoe 
vortex  structures  produced  by  each  element.  Therefore,  as  is  discussed  in  the  previous 
chapters,  the  cube  at  45°  has  the  largest  magnitude  in  FA  followed  by  the  cube  at  90°  and 
the  axisymmetric  elements.  The  flow  angle  created  by  the  cube  at  45°  indicates  that  the 
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local  cross-flow  is  nearly  50%  of  the  local  streamwise  velocity  value.  The  range  of  the 
flow  angles  across  the  measurement  plane  is  a  good  indicator  of  the  three-dimensionality 
of  the  flow.  Thus,  there  is  highly  three-dimensional  flow  near  the  wall  (y+<10)  in  the 
vicinity  of  the  horseshoe  vortex  structures.  This  three-dimensionality  in  the  flow  tends  to 
dissipate  with  distance  from  the  vortex  structure. 

Contours  of  the  FA  at  a  downstream  measurement  location  of  x/d=10.0  can  be  seen  in 
Figures  E.4  and  E.5.  As  the  flow  progresses  downstream  the  strength  of  the  horseshoe 
vortex  structure  dissipates  and  causes  a  decrease  in  flow  angle  magnitude.  The  cube 
elements  still  produce  the  most  significant  cross-flow  at  this  location  with  the  cube  at  45° 
being  the  largest.  The  axisymmetric  elements  are  still  very  comparable  to  one  another 
but  the  large  grooved  element  and  the  hemisphere  tend  to  have  the  lowest  values  of  cross- 
flow  at  this  downstream  location. 


E.2  Flow  Gradient  Angle 

The  flow  gradient  angle  (FGA)  indicates  the  flow  strain  change  associated  with  the  flow 
and  is  calculated  using  Equation  E.2  below. 


FGA  =  tan"' 


dW/dy 

dU/dy 


(E.2) 


A  magnitude  in  the  FGA  larger  than  45°  indicates  that  the  flow  in  the  spanwise  direction 
changes  faster  than  the  flow  in  the  streamwise  direction.  Contour  plots  of  the  flow 
gradient  angle  at  the  x/d=2.75  plane  (x/d=1.36  for  fine  grooved  element)  can  be  seen  in 
Figures  E.6  and  E.7.  The  fine  grooved  element  at  x/d=1.36  shows  very  large  values  of 
FGA  near  the  peak  of  the  element  which  indicate  that  the  spanwise  flow  changes 
significantly  faster  than  the  streamwise  flow.  Again,  the  cube  elements  have  the  most 
effect  on  the  FGA  followed  by  the  axisymmetric  elements  which  have  comparable 
magnitudes.  In  these  contour  plots  it  is  possible  to  see  the  regions  where  there  is  the 
fastest  change  in  flow  direction.  As  the  distance  from  the  wall  increases,  a  region  of 
large  FGA  is  seen  around  15<y+<40.  In  this  region  dW/dy  decreases  rapidly  as  the 
distance  to  the  core  of  the  vortex  gets  less  and  less.  Likewise,  dW/dy  increases  rapidly 
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as  you  move  away  from  the  center  of  the  vortex  in  the  positive  y-direction.  It  is  in  this 
region  that  negative  FGA’s  exist,  in  the  neighborhood  of  100<y+<300.  Contour  plots  of 
the  FGA  in  the  y-z  plane  at  a  location  of  x/d=10.0  can  be  seen  in  Figures  E.8  and  E.9.  As 
is  expected  the  magnitudes  of  the  FGA’s  have  decreased  to  values  close  to  50%  of  the 
upstream,  x/d=2.75,  magnitudes.  The  higher  magnitude  regions  have  moved  farther 
away  from  the  wall  and  centerline  as  compared  to  the  upstream  measurement  plane.  This 
result  is  quite  obvious  due  to  the  fact  that  the  FGA’s  are  directly  related  to  the  horseshoe 
vortex  structures.  These  vortices  have  previously  been  shown  to  move  away  from  the 
wall  and  centerline  as  they  propagate  downstream. 


E.3  Shear  Stress  Angle 

The  shear  stress  angle  (SSA)  is  used  to  describe  the  direction  of  the  effective  turbulent 
stress  in  the  flow  field  {Tang,  2004).  The  definition  of  the  SSA  can  be  found  in  Equation 
E.3. 


SSA  = tan  1 


\  —  vw | 


(E.3) 


The  calculated  uncertainty  in  the  shear  stress  angle  is  5SSA  ~  ±2.4°;  see  Appendix  C, 
Section  C.6.  For  similar  reasons  expressed  in  Section  E.l,  the  SSA  should  be  zero  along 
the  centerline.  Likewise,  due  to  the  same  proximity  issues  previously  discussed  the  SSA 
is  not  exactly  zero  at  the  near  element  locations  as  is  evidenced  by  Figure  E.10.  At  the 
downstream  locations  the  SSA  is  seen  to  collapse  to  zero  except  at  very  near-wall 
locations.  Contour  plots  of  the  SSA  in  the  y-z  plane  at  x/d=2.75  (x/d=1.36  for  fine 
grooved)  and  x/d=10.0  for  all  elements  can  be  found  in  Figures  E.l  1  through  E.l 4.  There 
are  higher  values  of  SSA  located  near  the  peaks  and  edges  of  each  element  indicating 
regions  of  large  directional  changes  in  the  turbulent  stresses  produced  by  the  roughness 
elements’  presence  in  the  boundary  layer.  These  areas  of  significant  direction  change  in 
the  effective  turbulent  stress  have  diminished  by  the  x/d=10.0  measurement  plane. 
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Figure  E.l:  Mean  flow  angle  variation  along  the  centerline,  where  FA  =  tai\-'(W/U) 
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Figure  E.2:  Contours  of  mean  flow  angle  (FA)  in  the  y-z  plane  at  x/d=2.75  (Fine  Grooved  x/d=1.36) 
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Figure  E.4:  Contours  of  mean  flow  angle  (FA)  in  the  y-z  plane  at  x/d=10.0 
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Figure  E.6:  Contours  of  flow  gradient  angle  (FGA)  in  the  y-z  plane  at  x/d=2.75  (Fine  Grooved  x/d=1.36) 
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Figure  E.7:  Contours  of  flow  gradient  angle  (FGA)  in  the  y-z  plane  at  x/d=2.75 
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Figure  E.8:  Contours  of  flow  gradient  angle  (FGA)  in  the  y-z  plane  at  x/d=10.0 
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Figure  E.9:  Contours  of  flow  gradient  angle  (FGA)  in  the  y-z  plane  at  x/d=10.0 
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Figure  E.10:  Shear  stress  angle  (SSA)  variation  along  the  centerline,  where 

SSA  =  tan"1  (-  vw/-  mv) 
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Figure  E.12:  Contours  of  shear  stress  angle  (SSA)  in  the  y-z  plane  at  x/d=2.75 
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Figure  E.13:  Contours  of  shear  stress  angle  (SSA)  in  the  y-z  plane  at  x/d=10.0 
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Figure  E.14:  Contours  of  shear  stress  angle  (SSA)  in  the  y-z  plane  at  x/d=10.0 
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Table  F.l  gives  the  y-z  coordinates  that  define  the  shape  of  both  the  large  and  fine 
grooved  elements.  The  shape  of  the  profiles  for  all  other  roughness  elements  can  be 
defined  via  their  respective  equations. 

Table  F.l:  Coordinate  locations  for  profile  shapes  of  grooved  elements 


Fine  Grooved  Element  1 

mm 

0 

0.0013462 

0.00002032 

0.00133604 

0.00004064 

0.00132588 

0.00006096 

0.00131572 

0.00008128 

0.00130556 

0.0001016 

0.0012954 

0.00012192 

0.00128524 

0.00014224 

0.00127508 

0.00016256 

0.00126492 

0.00018288 

0.00125476 

0.0002032 

0.0012446 

0.0002032 

0.00102616 

0.00022352 

0.001037844 

0.00024384 

0.001049528 

0.00026416 

0.001061212 

0.00028448 

0.001072896 

0.0003048 

0.00108458 

0.00032512 

0.001096265 

0.00034544 

0.001107949 

0.00036576 

0.001119633 

0.00038608 

0.001131317 

0.0004064 

0.001143 

0.0004064 

0.000924561 

0.00042672 

0.000936245 

0.00044704 

0.000947929 

0.00046736 

0.000959613 

0.00048768 

0.000971297 

0.000508 

0.000982981 

0.00052832 

0.000994665 

0.00054864 

0.001006349 

0.00056896 

0.001018033 

0.00058928 

0.001029717 

0.0006096 

0.001041401 

0.0006096 

0.000822961 

0.00062992 

0.000834645 

Large  Grooved  Element 

HUH 

0 

0.0009398 

0.0000508 

0.00097028 

0.0001016 

0.00100076 

0.0001524 

0.00103124 

0.0002032 

0.00106172 

0.000254 

0.0010922 

0.0003048 

0.00112268 

0.0003556 

0.00115316 

0.0004064 

0.00118364 

0.0004572 

0.001214121 

0.000508 

0.001244601 

0.000508 

0.000685801 

0.0005588 

0.000716281 

0.0006096 

0.000746761 

0.0006604 

0.000777241 

0.0007112 

0.000807721 

0.000762 

0.000838201 

0.0008128 

0.000868681 

0.0008636 

0.000899161 

0.0009144 

0.000929641 

0.0009652 

0.000960121 

0.001016 

0.000990601 

0.001016 

0.000431801 

0.0010668 

0.000462281 

0.0011176 

0.000492761 

0.0011684 

0.000523241 

0.0012192 

0.000553722 

0.00127 

0.000584202 

0.0013208 

0.000614682 

0.0013716 

0.000645162 

0.0014224 

0.000675642 

0.0014732 

0.000706122 

0.001524 

0.000736602 

0.001524 

0.000177802 

0.0015748 

0.000208282 

Appendix  F:  Element  Profile  Shape  Data 


-286- 


0.00065024 

0.000846329 

0.00067056 

0.000858013 

0.00069088 

0.000869697 

0.0007112 

0.000881381 

0.00073152 

0.000893065 

0.00075184 

0.000904749 

0.00077216 

0.000916433 

0.00079248 

0.000928117 

0.0008128 

0.000939801 

0.0008128 

0.000721361 

0.00083312 

0.000733045 

0.00085344 

0.000744729 

0.00087376 

0.000756413 

0.00089408 

0.000768097 

0.0009144 

0.000779781 

0.00093472 

0.000791465 

0.00095504 

0.000803149 

0.00097536 

0.000814833 

0.00099568 

0.000826517 

0.001016 

0.000838201 

0.001016003 

0.000619763 

0.001036323 

0.000631447 

0.001056643 

0.000643131 

0.001076963 

0.000654815 

0.001097283 

0.000666499 

0.001117603 

0.000678183 

0.001137923 

0.000689867 

0.001158243 

0.000701551 

0.001178563 

0.000713235 

0.001198883 

0.000724919 

0.0012192 

0.000736601 

0.0012192 

0.000518162 

0.00123952 

0.000529846 

0.00125984 

0.00054153 

0.00128016 

0.000553214 

0.00130048 

0.000564898 

0.0013208 

0.000576582 

0.00134112 

0.000588266 

0.00136144 

0.00059995 

0.00138176 

0.000611634 

0.00140208 

0.000623318 

0.0014224 

0.000635002 

0.0014224 

0.000416562 

0.00144272 

0.000428246 

0.00146304 

0.00043993 

0.00148336 

0.000451614 

0.00150368 

0.000463298 

0.001524 

0.000474982 

0.0016256 

0.000238762 

0.0016764 

0.000269242 

0.0017272 

0.000299722 

0.001778 

0.000330202 

0.0018288 

0.000360682 

0.0018796 

0.000391162 

0.0019304 

0.000421642 

0.0019812 

0.000452122 

0.002032 

0.000482602 

0.0020828 

0.000457202 

0.0021336 

0.000431802 

0.0021844 

0.000406402 

0.0022352 

0.000381002 

0.002286 

0.000355602 

0.0023368 

0.000330202 

0.00235712 

0.000320042 

0.00235966 

0.000315277 

0.00238506 

0.00029889 

0.00241046 

0.00027925 

0.00243586 

0.000255608 

0.00246126 

0.000226715 

0.00248666 

0.000190195 

0.00251206 

0.000140207 

0.00253746 

4.30303E-05 

0.00254 

0 

Appendix  F:  Element  Profile  Shape  Data 


-287- 


0.00154432 

0.000486666 

0.00156464 

0.00049835 

0.00158496 

0.000510034 

0.00160528 

0.000521718 

0.0016256 

0.000533402 

0.0016256 

0.000314962 

0.00164592 

0.000326646 

0.00166624 

0.00033833 

0.00168656 

0.000350014 

0.00170688 

0.000361698 

0.0017272 

0.000373382 

0.00174752 

0.000385066 

0.00176784 

0.00039675 

0.00178816 

0.000408434 

0.00180848 

0.000420118 

0.0018288 

0.000431802 

0.0018288 

0.000213362 

0.00184912 

0.000225046 

0.00186944 

0.00023673 

0.00188976 

0.000248414 

0.00191008 

0.000260098 

0.0019304 

0.000271782 

0.00195072 

0.000283466 

0.00197104 

0.00029515 

0.00199136 

0.000306834 

0.00201168 

0.000318518 

0.002032 

0.000330202 

0.002032003 

0.000111764 

0.002052323 

0.000123448 

0.002072643 

0.000135132 

0.002092963 

0.000146816 

0.002113283 

0.0001585 

0.002133603 

0.000170184 

0.002153923 

0.000181868 

0.002174243 

0.000193552 

0.002194563 

0.000205236 

0.002214883 

0.00021692 

0.0022352 

0.000228602 

0.00225552 

0.000218442 

0.00227584 

0.000208282 

0.00229616 

0.000198122 

0.00231648 

0.000187962 

0.0023368 

0.000177802 

0.00235712 

0.000167642 

0.00237744 

0.000157482 

0.00239776 

0.000147322 

0.00241808 

0.000137162 

0.0024384 

0.000127002 

Appendix  F;  Element  Profile  Shape  Data 


-288- 


Vita 


Jeremy  Bennington  started  his  undergraduate  career  at  the  University  of  Northern 
Colorado  (UNC)  in  Greeley,  CO  the  fall  of  1997.  While  at  UNC  he  studied  Pre- 
Engineering  and  played  baseball  at  the  varsity  and  semi-professional  levels.  In  the  fall  of 
1 999  he  transferred  to  Kansas  State  University  in  Manhattan,  KS  to  pursue  a  degree  in 
Mechanical  Engineering.  In  the  spring  of  2000  he  transferred  to  the  University  of 
Wyoming  located  in  Laramie,  WY.  He  graduated  cum  laude  in  the  spring  of  2002  with  a 
BS  in  Mechanical  Engineering  while  majoring  in  Mechanical  Engineering  and 
Mathematics.  In  pursuit  of  higher  education,  he  entered  the  graduate  program  in  the 
Aerospace  and  Ocean  Engineering  Department  at  Virginia  Polytechnic  Institute  and  State 
University  in  Blacksburg,  VA.  He  completed  the  requirements  for  a  Master  of  Science 
degree  in  Aerospace  Engineering  in  the  summer  of  2004. 


